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Abstract 



We present a strategy capable of describing basic features of the dynamics of crowds. The 
behaviour of the crowd is considered from a twofold perspective: both macroscopically and 
microscopically. We examine the large scale behaviour of the crowd (considering it as a 
continuum), simultaneously being able to capture phenomena happening at the individual 
pedestrian's level. We unify the micro and macro approaches in a single model, by working 
with general mass measures and their transport. 

We improve existing modelling by coupling a measure-theoretical framework with basic ideas 
of mixture theory formulated in terms of measures. This strategy allows us to define several 
constituents (subpopulations) of the large crowd, each having its own partial velocity. We 
thus have the possibility to examine the interactive behaviour between subpopulations that 
have distinct characteristics. We give special features to those pedestrians that are repre- 
sented by the microscopic (discrete) part. In real life situations they would play the role of 
firemen, tourist guides, leaders, terrorists, predators etc. Since we are interested in the global 
behaviour of the rest of the crowd, we model this part as a continuum. 

By identifying a suitable concept of entropy, we derive an entropy inequality and show that 
our model agrees with a Clausius-Duhem-like inequality. From this inequality natural re- 
strictions on the proposed velocity fields follow; obeying these restrictions makes our model 
compatible with thermodynamics. 

We prove existence and uniqueness of a solution to a time-discrete transport problem for 
general mass measures. Moreover, we show properties like positivity of the solution and con- 
servation of mass. Although afterwards we opt for a particular form, our results are valid for 
mass measures in their most general appearance. 

We give a robust scheme to approximate the solution and illustrate numerically two-scale 
micro-macro behaviour. We experiment with a number of scenarios, in order to capture the 
emergent qualitative behaviour. 

Finally, we formulate open problems and basic research questions, inspired by our modelling, 
analysis and simulation results. 

Keywords: Crowd dynamics; Social and behavioral sciences; Conservation laws; Micro- 
macro models; Mass measures; Thermodynamics; Mixture theory; Social networks; Initial 
value problems; Simulation 
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I want to suggest that even with our woeful ignorance of why 
humans behave the way they do, it is possible to make some 
predictions about how they behave collectively. 

Philip Ball, Critical Mas^ 



Quotation taken from 0, p. 6. 
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Chapter 1 

Introduction 



I would like to use this first section to introduce the central challenge treated in this thesis: 
the modelling of the behaviour of crowds. This section is also used to place this subject in a 
broader context. 

1.1 Process leading to this thesis 

In 2010 I was offered the opportunity to take part in the Honors Program Industrial and Ap- 
plied Mathematics: a new initiative of the Department of Mathematics and Computer Science 
(Eindhoven University of Technology) to challenge "the best students of lAM" by working 
in an academic, scientific setting. During these months I analyzed and extended a number 
of currently existing approaches to the modelling of crowd dynamics, under supervision of 
Dr. Adrian Muntean and Prof.dr. Mark Peletier. Since this research project took place in 
a relatively short period of time, we ultimately had to conclude that there were many more 
things to be explored in this field. As a natural consequence, we thus decided to dedicate my 
master's thesis to the same topic. 

The next part of this introduction is intended to emphasize that we have not been dealing 
with an irrelevant problem. Indeed, a series of everyday-life events have shown the impor- 
tance of being able to predict the dynamics of a crowd. Several attempts to capture human 
behaviour in mathematical models have found their way into useful applications. In the 
following sections these statements are discussed in more detail. 

1.2 Societal background 

In his 1895 book La psychologic des foules, the French social psychologist Gustave Le Bon 
wrote: "L'age ou nous entrons sera veritablement I'ere des foules" Although Le Bon could 
probably at the time (the fin de siecle) not quite foresee what the twentieth century would 
look like, his words have turned out to be prophetic. In their most literal sense - disregarding 
any political or metaphorical meaning that Le Bon might have intended to attach to them - 
they describe what our world has become. Indeed, during the twentieth century the global 

^"The age we are about to enter, will truly be the era of crowds". The French quotation is taken from [39] , 
p. 12, a republication of the original work published in 1895. 
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population has been ever-increasing. Occasions involving large numbers of people in crowded 
areas have become familiar sights. We experience such situations both under everyday urban 
conditions, and during large-scale manifestations with huge audiences (which take place oc- 
casionally) . 

Most of the "everyday situations" mentioned above happen in a 'normal' setting, that is, 
without the people being in a state of panic. Although this tends to sound reassuring, it 
does by no means imply that no special attention is to be paid. In order to design safe and 
comfortable public space, the presence of pedestrians cannot be disregarded. The behaviour 
of individual pedestrians, possibly clustering together to form larger crowds, is an important 
factor in urban design and traffic management. Indeed, unsafe and uncomfortable situations 
are often related to congestion and high densities. To assess the quality of the pedestrian 
environment, the following questions may be taken as a guideline (see |40] ) : 

• Are routes direct, leading the people where they actually want to go? 

• Is it easy to find and follow a certain route? 

• Are crossings easy to use; how long do pedestrians have to wait before they can cross a 
road? 

• Are footways well-lit, and sufficiently wide; what obstructions are there? 

• If pedestrians, cyclists and vehicles are mixed, does the typical speed of road users allow 
for this? 

These questions also indicate which measures can prevent the aforementioned cases of un- 
safety /discomfort. The necessity of these measures depends on the typical scale of pedestrian 
traffic (i.e. the 'crowdedness') that the area has to deal with. The effectiveness of each specific 
measure is also determined by its interplay with the others. 



The comments made here, should be extended to the semi-public domain. To be more con- 
crete, we do not only need to consider the behaviour of people walking in the street, on 
sidewalks and in parks, but also in railway stations (see Figure 1.1), sports stadiums and 
shopping malls. In fact, the investigation of crowd dynamics is essential in any setting in 
which a person's desired motion is impeded by the presence of other people, to prevent that 
inconvenience escalates into danger. 



If proper understanding of, or response to a crowd's behaviour is lacking, events in the 
(recent) past have shown us what the consequences are. These are dramatic situations during 
which people's lives are at stake. An example is the pilgrimage (Hajj) to Mecca and other holy 
places in Saudi Arabia, that Muslims perform every year during the Dhu al-Hijjah month. 
Since all adult Muslims are required to perform such Hajj at least once during their lifetimes 
(unless physical or financial circumstances make this impossible), huge numbers of people 
(several millions) gather every year during a relatively short period of time. Near the town 
of Mina, pilgrims as a ritual throw pebbles at three stone pillars (jamarahs), representing 
the devilj^ A special bridge (the Jamarat Bridge) was built is the 1960s to allow pilgrims to 

God/ Allah had commanded Abraham/Ibrahim to sacrifice his son Ishmael, but Satan urged Abraham to 
disobey God. The 'stoning of the devil' ritual is to commemorate Abraham rejecting the temptation of the 
devil. [2] 
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Figure 1.1: Passengers rushing towards the exit of Eindhoven Railway Station, after alighting 
from two trains that have just arrived. At the transition from the platform to the staircase a 
bottleneck occurs. (Photo: Joep Julicher) 



take part in the ritual at two floor levels. However, during the 1424 AH (2004) and 1426 AH 
(2006) Hajjs several hundreds of people were trampled to death due to overcrowding at the 
bridge. As this urged authorities to take action, crowd experts were asked to investigate the 
situation, and give recommendations for improvements; see [6l I30j. 

The site was completely reconstructed afterwards. The three pillars were replaced by wider 
and taller oblong concrete structures, to allow more pilgrims at the same time, and pre- 
vent them from accidentally throwing pebbles at each other. Crowds flow around these new 
'pillars' more easily, by which congestion is reduced. A new multi-storey bridge was built 
with better entrance facilities and more emergency exits. Bottlenecks were removed. More- 
over, the protocols for crowd management by stewards and other officials were reconsidered. [1] 

Unfortunately, crowd disasters do not always happen far away from us. Recently, in July 
2010, 21 people died at the Loveparade in the German city of Duisburg, at a distance of less 
than 100 kilometers from Eindhoven. This festival took place in a closed-off area, which could 
be entered via a small number of tunnels. Each of these ended at a common ramp/staircase, 
eventually leading to the festival ground. Due to overcrowding at the bottom end of the 
staircase, a stampede occurred. All lethal victims died because of suffocation. 
After these disastrous events, there was a lot of criticism on the safety precautions that had, 
or had not, been taken. The organization of the Loveparade was blamed for having provided 
too few emergency exits to the Loveparade area. Moreover, the number of security agents 
was smaller than promised. On the other hand, local authorities are said to have ignored the 
police's objections, because they were so keen on having this prestigious event in their town. 
For more information, see e.g. [H El [3]. 

The above emphasizes the importance and urgency of being able to properly describe the 
dynamics of a crowd. This is a prerequisite for predicting crowd behaviour, which in its turn 
is needed to anticipate life-threatening situations. 
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1.3 Modelling approaches 

Ever since the Renaissance scientists have wondered how human behaviour can be captured 
in mathematical formalism (cf. e.g. [7]). We focus now in particular on the developments of 
the last decades, during which the description and analysis of a crowd's motion has gained 
the attention of the scientific world. This is mainly due to the gradual increase of the number 
of large-audience events being organized, and the accidents happening at such events. An 
illustration of these developments has been given above. Several models have been developed 
and explored to catch the crowd-related phenomena we experience in real life, in a scientific 
(mathematical) framework. These models were treated both analytically and numerically, 
and were based on a number of distinct perspectives. 

Focussing on what happens at the level of an individual pedestrian, one obtains discrete, 
agent-based or microscopic models. In this approach the dynamics of each single person in 
the crowd are modelled and traced. This perspective is adopted by Dirk Helbing et al. within 
the framework of a so-called social force model, see e.g. [31]. Roughly speaking, each pedes- 
trian is driven via Newton's Second Law, by means of a social force that depends on the 
presence and social behaviour of other people. Simulation is their main tool for obtaining 
qualitative and quantitative information. 

On the other hand, one can also 'zoom out', considering a crowd on a more global scale. 
Working on this macroscopic level, one uses densities, rather than individual pedestrians. 
Via this approach one ends up in what is essentially a fluid-dynamic setting. A specific 
property of these models is that we cannot capture local interactions any more. Macroscopic 
modelling is thus only useful if we are interested in average characteristics of the crowd and 
its motion. This way of modelling is for instance adopted by Bertrand Maury et al. in [4lj . 
where initially a gradient flow structure is proposed. Their results are derived analytically. 

Up to the initial conditions (which might be generated by means of a random sampling), 
the models described above are fully deterministic. 

In this thesis we focus an approach that differs in nature from the aforementioned ones. 
Our motivation is that we do not want to be forced to choose either of the two perspectives; 
we aim at 'marrying' the two perspectives in one single model. This strategy was described 
by Benedetto Piccoli, Andrea Tosin et al. in |I451 I17j . In most of the fluid-dynamic mod- 
els of pedestrians, non-linear hyperbolic conservation laws appear. Certain analytical and 
numerical problems inherently arise when treating these PDEs, especially in more than one 
space dimension. For instance the solution might not be unique, shock waves can occur, it 
is important whether the corresponding flux is convex or not, boundary conditions are diffi- 
cult to impose, etc. Numerically, positivity of the density might not always be guaranteed. 
In order to circumvent these problems, Piccoli et al.'s work takes place in a (time-discrete) 
measure-theoretical framework. This thesis is mainly built upon the fundaments Piccoli cum 
suis provides. 

However, we want to be able to cover a much wider range of situations, for instance a crowd 
consisting of several distinct subpopulations (rather than only one population) , each of which 
is willing to move with its own desired velocity. Moreover, we draw parallels with mixture 
theory and thermodynamics, and extract inspiration from those fields. 
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The microscopic, macroscopic and micro-macro approaches described above, are not the only 
existing strategies. See e.g. [H] (pp. 212-214) for a systematic classification of models. Par- 
tially based on [H], we describe a number of different perspectives. 

At an intermediate level between micro and macro are so-called mesoscopic models. These 
models do not distinguish between individuals and hence cannot trace their trajectories. 
However, individual behaviour can be specified. The description is in terms of probabil- 
ity densities: mostly, the probability to find an individual with specified speed in a certain 
location at a certain time. To give an example. Dirk Helbing describes vehicular traffic in 
such framework (also indicated by 'Boltzmann-like' or 'gas-kinetic' models) in [29j. Moreover, 
in [28] he derives hydrodynamic equations for pedestrians, based on the Boltzmann-approach. 

Up to now, all models have mainly been based on considerations related to physics. This 
mostly matches our point of view. Other views are (of course) also possible. 

The fact that we do not fully understand the underlying mechanisms of human behaviour, is 
most naturally incorporated in a model by adding stochasticity. A small amount of external 
noise can be added to avoid undesired situations. One should think here for example of a 
situation in which two individuals are positioned 'head-on' and both want to move straight 
ahead. In a model it might happen that a deadlock occurs, although in reality the two most 
likely just move aside slightly, and get passed one another. A bit of stochastic fluctuation 
forces a breakthrough if the described deadlock happens. Away from such configurations, the 
influence of the noise is only small. 

Fully stochastic models are different in the sense that they do not just add random perturba- 
tion to intrinsically deterministic dynamics. Here, the underlying decision-making processes 
are influenced directly by random effects. If you take the deterministic limit in these models 
(i.e. vanishing stochastic effects) the overall phenomena are completely different from those 
in the stochastic 'normal' case. This makes fully stochastic models intrinsically different from 
models that contain external noise only. 

In a cellular automata model all variables are discrete. The crowd is described at the in- 
dividual's level. The spatial domain is subdivided into a number of cells, where each cell 
possesses a state variable (e.g. = 'empty', 1 = 'occupied'). At specifled (discrete) points 
in time, the model decides which cells are occupied in the subsequent generation; the state 
variables are then updated accordingly. For reality-mimicking models the total number of 
occupied cells is constant. Typically, the update is rule-based, i.e. each occupied cell (read: 
particle/pedestrian) makes a decision where to go, based on the current situation and its 
goals. These rules are often supported by psychological arguments. Possibly, some stochastic 
effects are included. The update can either be executed in parallel (for all cells simultane- 
ously), or for randomly selected individuals only. 

Serge Hoogendoorn and Piet Bovy model individuals as players in a differential game. This 
approach was flrst applied to driving behaviour in traffic fiows in |34| . The ideas presented 
therein were subsequently extended to pedestrian dynamics; see [33]. Vehicles and/or pedes- 
trians are assumed to maximize their expected success or profit, or follow a Zipfian principle 
of least effort. In such a game-theoretical approach, individuals are allowed to modify their 



8 Modelling Crowd Dynamics 



Technische Universiteit Eindhoven University of Technology 



control decisions based on their observations, and on predictions of the behaviour of other 
players in their neighbourhood. 



1.4 A few people in the field 

Throughout 2010 and 2011 Adrian and I have been in contact with many people, whose area 
of expertise was related to our research. When meeting them, we were mainly interested in 
what kind of questions these experts would like to be answered by crowd models. Moreover, 
we obtained insight in the wide range of possible applications. 

The person who provided us with the actual idea for studying the topic of crowd dynamics 
is Prof. Chris Budd (University of Bath, UK). Mark Peletier invited me to meet Prof. Budd 
during his stay in Eindhoven, Spring 2010. At that time, I needed to decide what I wanted 
to do in my Honors Program, and these two people suggested the idea (that has eventually 
led to this thesis). In January 2011, Adrian Muntean and I visited Mark Peletier during his 
sabbatical stay in Bath, and we met Prof. Budd again. He was the one drawing our attention 
to two-scale phenomena in nature (birds and fish), by showing us a couple of fascinating BBC 
movies 13 



Prof.dr.ir. Serge Hoogendoorn, Dr.ir. Winnie Daamen and Mario Campanella MSc (Delft 
University of Technology) work on modelling/simulation of pedestrian/traffic dynamics, and 
on the calibration of their models by real-life experimentsj^ Winnie Daamen gave a talk in 
Eindhoven (during the CASA colloquium) in September 2010. Adrian Muntean and I paid a 
visit to their Transport & Planning Department in January 2011. 

They have developed a software package called Nomad, which is a microscopic simulation 
tool that can e.g. be used to assess the geometry of infrastructure. The underlying model is 
the one described in [33] . The behaviour of the individuals is comprised in their acceleration, 
which consists of an uncontrollable and a controllable part. The uncontrollable part is due 
to physical interactions with their surroundings, that is, with other individuals and obstacles 
(part of the geometry of the domain). The controllable part contains the tendency to mini- 



mize the walking cost (cf. Section 1.3). 



Calibration of the simulation is done by observing real-life traffic and pedestrian flows (as far 
as privacy regulations permit to do so) and laboratory experiments; see e.g. [HI |35j. These 
laboratory experiments take place in standardized environments, and participants' positions 
are obtained by video image analysis. An example of such standard setting (which is very 
well-known in the field) is the counterflow or bidirectional flow: a narrow corridor in which 
two groups of people want to move in opposite directions. In this experiment, one expects 
to observe a specific phenomenon of self-organization: lane-formation. That is, people tend 
to align in such a way that they just follow someone going in the same directionj^ In Figure 



1.2 a schematic impression of a bidirectional flow is given. Another standard scenario is a 



bottleneck. It occurs for instance if a corridor suddenly gets narrower. At this transition, in- 

of this thesis. 



4.8 



''These phenomena will be addressed in Section 
^See www. pedestrians. tudelft .nl for more information on either of these two areas of expertise. 

video of a bidirectional flow experiment is available at www . youtube . com/watch?v=J4J__100V2E. Lanes 
are formed without the participants being instructed to do so (note that these lanes are unsteady). The 
experiment was performed in connection with the German Hermes project; in Delft similar experiments are 
done. 
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dividuals get clogged up in circular structures (arches) which are hard to break; this scenario 
is indicated in Figure 1.3 See also [48j, pp. 418-419, for a description of both settings and of 



the emerging phenomena. 



Figure 1.2: Schematic drawing of a counterflow scenario with lane- formation. At both ends 
of the corridor individuals are supplied; they intend to walk to the other side. Halfway, four 
lanes emerge. 




Figure 1.3: Schematic drawing of arches at a bottleneck, the flow is obstructed there. 

Ing. Gilian Brouwer works for Adviesburo Nieman: a company specialized in consultancy 
on quality, safety and building physicsj^ He is specialized in fire safety engineering. At Nie- 
man, several software tools are applied to simulate the evacuation of people in case of fire in 
a building. In case a building does not comply with Dutch regulations ('Bouwbesluit'), simu- 
lations are needed to make feasible that nevertheless an equivalent level of safety is provided. 
In |38| a comparison is made between two of the software packages that are used by Nieman 
(FDS + Evac and Simulex). Both packages are based on a social force-like model, where 
FDS + Evac also includes the effect of smoke. An individual's desired route is based on the 
geometry of the domain and the configuration of other pedestrians. The route is chosen such 
that it is expected to take a minimal amount of time to reach the destination. In Simulex 
however, in each spatial position the shortest route to the exit is precalculated (considering 
the distance). The evacuation is assessed by recording the times at which pedestrians leave 
the building, and consequently calculating exit flow rates or average exit times. 

Prof.dr.ir. Bauke de Vries is the chairman of the Urban Management &: Design Systems 
group at the Department of Architecture, Building and Planning (Eindhoven University of 
Technology). We met him in February 2011 and spoke about his working interests and area 
of expertise. He mainly focuses on interior design of buildings, where for instance routing 
and the positions of obstacles are taken into account. In the past, he also performed real-life 
experiments in his department, to obtain data about the walking behaviour of his coworkers 

''www . nieman . nl 
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within the office space. During our conversation, he also explained the work of Prof.dr. Harry 
Timmermans, who is a member of the same group. He is more specialized in the design of 
exterior urban space, such as shopping malls. An example can be found in |19j . in which the 
two experts cooperated. Their simulation models work with agents (read: pedestrians) that 
have an agenda and an updating mechanism for their 'to do'-list; the individuals behave and 
move accordingly. 

We also talked with people about crowds in a broader context. An example is Dr. Petru 
Curseu who is an expert in group processes and decision-making at Tilburg University. To 
get an impression of his area of interest, see [18]. When he visited us in Eindhoven (November 
2010), he showed his interest in microscopic pedestrian models. He proposed to use them as 
a metaphor for more abstract societal phenomena. For instance, if there is strongly repul- 
sive interaction between a large group and a small group, the question arises whether this 
enables the large group to leave a room relatively faster than the small group. According to 
Dr. Curseu, a parallel can be drawn between this pedestrian simulation and a majority in a 
country that discriminates a minority. The majority group turns out to have better access to 
resources. 

I took part in the Study Group Mathematics with Industry 2011 (organized at the Vrije Uni- 
versiteit Amsterdam, 24-28 January 2011). The problem was provided by Chess An ad hoc 
wireless network is considered, in which each node can broadcast and receive messages. How- 
ever, it is uncertain whether sent messages reach a receiver, and if yes, how many/which nodes 
are reached. The aim was to let each node send, receive and process 'intelligent' messages such 
that it can estimate the number of nodes in its neighbourhood and in the total network. The 
application possibilities of such unreliable wireless networks are closely linked to the subject 
of this thesis. This was also pointed out to us at an earlier stage by Prof.dr. Fabian Wirth 
(University of Wiirzburg, Germany), an expert in logistic networks. Typically, each node is 
a simple microprocessor, with limited computational power, that measures a property of its 
environment, like temperature or position. A large number of such sensors can be used to 
detect forest fires, but also to obtain information about e.g. bird flocking, social behaviour 
and crowd dynamics. The proceedings of the study group are to be expected. 
In March 2011 Adrian Muntean and I had a conversation with Dr. David Fedson, a medical 
expert in the area of epidemics. We talked about Malcolm Gladwell's The Tipping Point, 
which treats phase transitions that take place everywhere around us. Certain events suddenly 
take off in social behaviour (fashion), epidemics of infectious diseases, and in animal groups. 
For the latter see [36]. Another topic he brought to our attention is synchronization. A 
connection between synchronous behaviour in nature and in the financial world is made in 
|47j . Prof. Steven Strogatz also explains this phenomenon in a fascinating wayF] 



1.5 Content and structure of this thesis 

In this thesis we focus on the mathematical background of crowd dynamics models. This 
means that the emphasis will not be on direct applications. We are mainly interested in iden- 
tifying underlying structure and mechanisms, and in understanding the coupling between the 
micro- and macro-scale. We work in a fully deterministic setting. For more information in 

^www. chess.nl 

*www. youtube . com/watch?v=aSNrKS-sCEO 
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the direction of stochastic modelhng approaches (such as stochasticahy interacting particles), 
the reader is referred e.g. to [18] . 

We first describe the theoretical fundaments on which we build our model. In Section [2] 
basic measure-theoretical concepts are introduced. The most important parts are the (re- 
fined) Lebesgue decomposition and the Radon-Nikodym Theorem. 

Section [3] is dedicated to mixture theory, a branch of continuum mechanics that turns out 
to be suitable for describing a crowd consisting of multiple subpopulations. The correspond- 
ing theory is presented in measure-theoretical language. It provides us with conservation of 
mass equations, that are the basis for our model. In Section [3] we also cover aspects of ther- 
modynamics that are used later to derive an entropy inequality for a particular crowds setting. 

In Section [4j we apply the modelling ideas of the preceding sections to obtain a model for 
the dynamical behaviour of a crowd. After a number of (mainly formal) calculations we first 
present a continuous-in-time model and the accompanying entropy inequality. We also pro- 
pose a velocity field, governing the dynamics. From this model, we derive a discrete- in-time 
version. All our mathematical analysis concerns this time-discrete model. The main result 
is a proof of global existence and uniqueness of a time-discrete solution. We also derive a 
discrete-in-time equivalent of the entropy inequality. 

Section [5] contains the full description of the numerical scheme and parameter setting, which 
we use to simulate our two-scale crowd setting. The aim of the simulation part is to inves- 
tigate the basic patterns produced by the interplay between one or two individuals and a 
macroscopic crowd. Section [6] contains our simulation results. 

In Section [7] we review the work done and the obtained results. Moreover, we give sug- 
gestions for future work and pose a few basic questions that are still open. 

Our first attempt [23] of modelling pedestrians using the approach described in this the- 
sis, was published in Nonlinear Phenomena in Complex Systems. Its content can be found in 
Appendix |F] 
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Chapter 2 

Measure theory 



2.1 Measures and their Lebesgue decomposition 

Let [Q,B{Q)) be a measurable space, where 7^ O C M*^. Here, B{ft) denotes the cr-algebra 
of Borel subsets of $7. Suppose that fi and A are positive, finite measures defined on B{Q). 

Definition 2.1.1 (Absolutely continuous measures). The measure is said to be absolutely 
continuous with respect to the measure A if for any fi' G i3{^l), A(O') = implies //(O') = 0. 
Notation: 

/u < A. 

Definition 2.1.2 (Singular measures). The measures /i and A are said to be mutually singular 
if there exists a B £ B{0,) such that 

/x(J7\B) = X{B) = 0. 

Notation: 

/i _L A. 

Although the relation n ± X is symmetric, it is also often said that 11 is singular with respect 
to A. 



Definition 2.1.2 is taken from [22], p. 40. An alternative definition is given by [16], p. 120, 



which we will now show to be equivalent. 
Lemma 2.1.3. The following two statements are equivalent: 
(i) There exists a B £ B{0,) such that \ B) = X{B) = 0. 

(a) There are disjoint A^.A^ G B{Vl) such that = ^{Q.' n Ai) and X{Vl') = X{Q.' n A2) 

for all n' G B{n). 

Proof. 1. Assume that holds. Define Ai := B and A2 := ^\B (these sets are obviously 
disjoint). For all O' G B{Cl) we have (using property ^) 

fi{n' \ Ai) ^ \ Ai) = \b) = o, 

and 

x{^}' \ A2) ^ x{n \ A2) = x{n \{n\ b)) = x{b) = o. 
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It follows that 

n{n') = n Ai) + \ Ai) ^ fi{n' n ^i), 

and 

\{n') = x{n' n A2) + a(!^' \ A2) ^ x{n' n A2). 

Since ft' D Ai and 0,' D A2 are subsets of 0,' , it is clear that 

fjL{n') > ^{Q'nAi), 

and 

x{n') ^ x{n'r\A2), 

thus we have that 

H{n') = /x(Q'n^i), 

and 

A(O') = A(S7'n A2). 
2. Assume that (|ii]) holds. Define B := Ai, then 

\ B) = \ B) n Ai) = \ Ai) n Ai) = = o, 

and 

X{B) = X{B n A2) = X{Ai n A2) = A(0) = 0. 

□ 

In the following theorem we relate any positive, finite measure to absolutely continuous 
and singular measures. 

Theorem 2.1.4 (Lebesgue decomposition). If fj, and X are positive, finite measures defined 
on B{^}), then there exists a unique pair of positive, finite measures ^ac (ind defined on 
B{Q), such that: 

(i) fJ. = fJ.ac + fJ-s, 
(a) fiac < A, 

(Hi) fis -L A. 

We call fj,ac the absolutely continuous part and fig the singular part of fj, w.r.t. A. The pair 
{iXacilJ^s) is called the Lebesgue decomposition of fi w.r.t. X. 

Proof. The proof of Theorem |2.1.4| is given in Appendix □ 

Definition 2.1.5 (Discrete measures). The measure fi is said to be a discrete measure with 
respect to the measure X if there exists a countable set A := {xi, X2, . . .} C 1^ such that 

^i{n\A) = X{A) = 0. 

Lemma 2.1.6. Let fi be a positive, finite measure on B{0,). Then the following two statements 
are equivalent: 

(i) is a discrete measure with respect to the Lebesgue measure X'^ . 
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(a) There is a countable set {xi, X2, ■ ■ ■} C and a set of corresponding nonnegative co- 
efficients {ai,a2, ■ ■ ■} C M, such that fi = Yl^i'^i^Xi- Here Sx^ is the Dirac measure 
centered at Xi . 

Proof. 1. Assume that ^ holds. Let A := {yi,y2, . . .} C O be the collection of points such 
that \ ^) = X{A) = 0. Without loss of generality, assume that all elements of A 
are distinct. (In case not all elements of A are distinct, just delete yj from A if there is 
ayi £ A satisfying i < j and yi = yj.) 
Since \ ^) = 0, we have for any fi' G B{^}): 

nA)(^ fi{n') = fi{Q' r\A) + fi{n' \ A) ^ fi{n' r\A) + fi{n \a) = n A), 

thus = n A). For fixed fi', let J' be the index set, such that i £ J implies 

G Jl' n A. Since ^ A = is a disjoint union, it follows that 

/i(i7') = /i(J^'nA) = ^My.)- 

If we define oti := ^(y?) ^ 0, and write 1 for the indicator function, then the above is 
equivalent to 

oo 
i=l 

This is exactly the definition of the Dirac measure, so we can also write 

oo 
i=l 

2. Assume that <^ holds, and let the set {xi,X2, . . .} be called B. It follows by definition 
of the Dirac measure that \ B) = 0. i? is a countable collection of points, which 
(obviously) all have Lebesgue measure zero. Thus X'^{B) = Yli^i ^'^{^i) = 0- We 
conclude that /i is discrete with respect to A'^. 

□ 

Definition 2.1.7 (Singular continuous measures). The measure fi is said to be singular con- 
tinuous with respect to the measure A if fj,{x) = for all x £ ^1, and there is a B £ B{il.) such 
that 

/x(0\B) = X{B) = 0. 

Remark 2.1.8. We only consider measures defined on B{^}), i.e. we can only apply them to 
sets. If we write /i(x), we therefore actually mean //({x}). 

Theorem 2.1.9 (Decomposition of singular measures). If Hs is a positive, finite measure 
defined on B{0,), which is singular w.r.t. X, then there exists a unique pair of positive, finite 
measures fj,d and ^sc defined on B{Vt), such that: 

(i) fig = fid. + fJ-sc, 

(a) fid is a discrete measure w.r.t. X, 
(Hi) fisc is singular continuous w.r.t. X. 
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We call the discrete part and figc the singular continuous part of w.r.t. A. The singular 
continuous part is also called the Cantor part of the measure. 



Proof. The proof of Theorem |2.1.9 is given in Appendix [B] □ 



Corollary 2.1.10 (Refined Lebesgue decomposition). Assume the hypotheses of Theorems 



2.1.4 and 2.1.9. Then for any positive, finite measure /d, there are unique measures ^ac, l^d 
and fisc such that 

fJ- = lJ'ac + ll'd + fJ-sc, 

where fiac ^ A*) fJ-d is a discrete measure, and ^sc is singular continuous w.r.t. A. 



Remark 2.1.11. The statement \ B) = X{B) = (cf. Definitions 2.1.2 2.1.5 and 



2.1.7), is equivalent to: A(-B) = and fJ-{B) = fi{^}). This is due to the identity fJ-{0,) = 
H{n nB) + \B) = ii{B) + /^(O \B). if S is such that ^(5) = we cah B a set of 

full measure. 

2.2 Radon-Nikodym Theorem 

The Radon-Nikodym Theorem is of vital importance in this context, especially for the appli- 
cability of measure theory to mixture theory as arising in real-life situations mimicking the 
dynamics of crowds. The following theorem gives sufficient conditions for a measure to be 
expressed in terms of a density, with respect to another measure: 

Theorem 2.2.1 (Radon-Nikodym for finite measures). Suppose fi and A are positive measures 
on a measurable space (^1,^3(17)) such that < /^(ri) < 00, < < 00, and let fj, 

be absolutely continuous with respect to A. Then there exists a real, nonnegative, B{^)- 
measurable function h on Q such that 



= j hdX, for all n' G B{n). 



Proof. A comprehensible, well-structured proof of Theorem 2.2.1 can be found in □ 



Remark 2.2.2. The density h is often called Radon-Nikodym derivative and is denoted by 

Lemma 2.2.3. Assume the hypothesis of Theorem \2.2.1\ Then the Radon-Nikodym derivative 
is unique in the following sense: if both hi and /12 satisfy /i(f^') = Jq, h^dX, for all Vt' G B{Q), 
where k £ {1,2}, then hi = /i2 almost everywhere (w.r.t. X) infl. 

Proof. This statement can easily be verified. For any 17' G B{i}) we have that 
{hi - h2)dX = j hidX - j h2dX = /i(O') - ^(0') = 0. 

Thus, /ii — /i2 = almost everywhere w.r.t. A in Q, i.e. hi = /i2 A-a.e. in Q. □ 
Remark 2.2.4. Theorem 2.2.1|also holds for more general situations, for instance, including 



the case of signed cr-finitc^ measures (see [2T], e.g.). 

^The measure A is called (T-finite if there exists a sequence of sets {Ek}'kLi C B(Q) with — IJfet^i ^fe ^-^d 
|A(£fc)| < 00 (for all k e N); see [21], p. 269. 
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2.3 Properties of Radon-Nikodym derivatives 

In the fohowing lemma we state a number of useful properties of Radon-Nikodym derivatives: 



Lemma 2.3.1 (Basic properties). Assume the hypothesis of Theorem 2.2.1 on the measures 
V, [I, X, ui, V2, fJ-i and ^2- The Radon-Nikodym derivatives satisfy the following general 
properties: 

^ rr ^ .1 dv dfj, 

1. If u <^ fi <^ X, then -rr = -r-fv- 

dX dfi dX 

2. If X and g G L^(il), then J^, gdfj, = J^, 9^dX for all W G B{il.). 

dfi fdX\-^ 

3. If fi <^ X and A ^ //, then ^ = 

dX ydfiJ 

4. If (fij, S(r2j), /ij), i £ {1,2}, are two measure spaces with Vi <^ fii, i £ {1,2}, then 

, d(i/i (g) U2) dvi dv2 

vi®V2'^IJ'i® and — r = — . 

(g) 1^2) dfii dfi2 

Proof. The proof of Lemma |2.3.1| is given in Appendix [Cj □ 

2.4 Mass measure fi^ 

Here, we already give an indication of the particular measures we intend to use in the rest of 
this thesis. Let 17 C M*^ be a domain (read: object, body) with mass. For physically relevant 
situations, we consider d £ {1, 2, 3}. Let fim{^') be defined as the mass contained in Q' C ^l. 
Note that we use the concept of a mass measure very much in the spirit of [S] . 

Remark 2.4.1. As a rule, whenever we write Q' C 17, we actually mean that is such that 
Q' G B{0,). We assume fim to be defined on all elements of B{il). 



In Sections 2.4.1 and 2.4.2 we consider two specific interpretations of this mass measure 



depending on the localization of the information we wish to capture. 
2.4.1 Microscopic mass measure 

Suppose that contains a collection of point masses (each of them of mass 1), and denote 
their positions by {pk}k=i ^ ^) for A G N. We want the mass measure mm to be a counting 
measure with respect to these point masses, i.e. for all il' G B{Q) 

mmi^') = #{pk G n'}. (2.4.1) 



By (2.4.1 ) we mean that mm counts the number of individuals located in 17' (with correspond- 
ing, known positions pk G 17'). This can be achieved by representing mm as the sum of Dirac 
masses, with their singularities located at the points pk, G {1, 2, . . . , A^}: 



N 

mm — 

k=l 
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If rrim is defined as in (2.4.1) or (2.4.2), we call it a microscopic mass measure. Note that 



such measure satisfies the conditions in Definition 12.1.51 and is thus a discrete measure. The 



characteristics of such measure are illustrated by Figure 2.1 



n 




Figure 2.1: A microscopic mass measure counts the number of points pk that are contained 
in r^', a subset of il. In this case rrimi^') = 2 since p2 and lie in 



2.4.2 Macroscopic mass measure Mm 

Let us now regard a different mass measure M^- Assume that the following postulate applies 
to Mm-. 



Postulate 2.4.2 (Properties of Mm). 1- Mm > 0. 



2. Mm is a-additive. 



3. Mm is finite. 



4- Mm ^ A"' (with X"^ the Lebesgue measure in M.'^). 



By Parts [H [2] and [3] of Postulate [2A2l we have that Mm is a positive, finite measure on 
Q, whereas Part |4] implies that there is no mass present in a set that has no volume (w.r.t. 
A*^). We refer to a mass measure satisfying Postulate 2.4.2 as a macroscopic mass measure. 
Now Theorem 2.2.1 guarantees the existence of a real, nonnegative density p G L^^^{Q) such 
that 



Mmin') 



p{x)dX'^{x), 



for ah n' £ B{n). 



An illustration of such measure is given in Figure 2.2 
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n 




Figure 2.2: A macroscopic mass measure is used to obtain the mass in Q' C fi. Mm{^') 
depends on the density p of the area in W. The density is here indicated by grey shading. 



ModeUing Crowd Dynamics 19 



Chapter 3 

Mixture theory and 
thermodynamics 



In this section, we present some ideas from [IH H2| [50] regarding the theory of mixtures. We 
cast their description in a measm'e-theoretical framework. UnHke [111 I42j . here we mainly 



consider the Eulerian point of view. In Sections 3.1-3.5 we describe the mixture as a contin- 



uum; this is the classical and most common way of doing so (cf. |1H I42| e.g.). In Section 3.4 
we extend the ideas of the preceding sections to more general mass measures. 



3.1 Description of a mixture 

At this stage, a mixture is defined as a continuum consisting of a certain number of con- 
stituents (also called: components). Constitutive relations typically differ from component to 
component. 

Let Q C M*^ be the domain in which the mixture is located, and consider the measurable 
space (0,;S(f2)), where B{Q) is the cr-algebra of Borel subsets of Q. Suppose the mixture 
consists of v constituents, with index a S {1,2,... For each constituent a we define its 
volume and mass present in Q' € B{^1) at time t G (0,T) by r°(t,0') and ;u"(t,i7'), respec- 
tively. Note that in fact t is the time variable and can be understood as a parameter. The 
objects r"(t, •) and /^"(i, •) on B{Q) are for each t € (0, T) fixed measures. 
We assume that the following postulate applies: 

Postulate 3.1.1 (Properties of and /i°). For all a G {1,2, ... ,z^} and all t G (0,T) we 
postulate: 

1. ^ 0. 



2. r" is a -additive. 



3. r" is finite. 

4- ^ A*^ (with X'^ the Lebesgue measure in M.'^). 
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5. ^ 0. 



6. /i" is a -additive. 



7. is finite. 



H'^ <. r° 



Note that the concepts of Section 2.4.2 are incorporated in this postulate. 



Remark 3.1.2. It follows in a straightforward way from Parts [4] and [8] of Postulate 3.1.1 
that < A"* for ah a G {1, 2, . . . , i'} and ah t G (0, T). The absolute continuity of both r" 
and complies with the fact that we define mixtures at the continuum level. 

Remark 3.1.3. We characterize the presence of constituent a at a certain x by the presence 
of mass and volume here. Mathematically, this means that if at time t G (0, T) a subset 
ri' G B{Q) contains some fraction of constituent a, then both //"(i, O') > and r"(t, Q') > 
must hold. On the other hand, if a constituent a is not present in Q' , then //"(t, f2') and 
T°(t, $7') are both zero. Following this intuitive characterization, we assume r" <C /x" for all 
aG {l,2,...,u} and for all t G (0,r). 

Let r and // be the (time-dependent) total volume and mass measures, given by 



T{t,n') 



Y.r^{t,n'), 



a=l 



for ah n' G I3{n) and aU t G (0, T). 



(3.1.1) 



(3.1.2) 



a=l 



By the definitions of r and /i given in (3.1.1) and (3.1.2) and as a consequence of Postu- 
late 3.1.1, we have <^ t <^ X'^. Note that since we also assumed <C /x" for each a, it 



follows by definition that also t <^ fj. holds. 

For any t G (0,T), the Radon-Nikodym Theorem provides the existence of the unique non- 
negative densities (Radon-Nikodym derivatives) Q{t,-), Q"{t,-), 6°'{t,-) G L^rf(r2) (for each 
aG{l,2,...,z.}): 

e := 4^, (3.1.3) 



dr" 



(3.1.4) 
(3.1.5) 



The Radon-Nikodym derivative d" arising in (3.1.5 ) is called the volume fraction of component 



a at time t. As a result of (3.1.1): 

V 

^6*° = 1, a.e. w.r.t. r, and for ah t G (0,r). 
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Note also that due to Part [T] of Lemma 2.3.1 we have that for all a G {1, 2, . . . , z^}, and for 
all t E (0, T) 

9" = re, a.e. w.r.t. A'^. 

Also for the mass measures unique Radon-Nikodym derivatives exist. For all a G {1, 2, . . . , z^} 
and for all t G (0, T), we define 



P ■-- 



dT 
d^i°' 



(3.1.6) 
(3.1.7) 
(3.1.8) 



We call p the density of the mixture, the partial density of component a, and p" the 
intrinsic density of that component. By (3.1.2) we have that 



^p", a.e. w.r.t. r, and for all t G (0,r). 



(3.1.9) 



a=l 



From (3.1.3 )-( 3. 1.5) and (3.1.6)- (3. 1.8) a number of identities can be derived (again via Lemma 



2.3. 1| Part[l|. For ah t G (0,T) we have 

-a ^ ^aga^ a.e. W.r.t. T, for aU a G {1, 2, . . . , I/}, 
p»Q = p^Q'', a.e. w.r.t. A"', for ah a G {1,2,...,!/}. 



(3.1.10) 
(3.1.11) 



Remark 3.1.4. Note that the latter p°'Q equals the Radon-Nikodym derivative dpf^ /dX^. In 
practice it will be much more convenient to work with A'* than to work with r. Therefore we 
define 

P":=^=P"0, (3.1.12) 
and henceforth refer to (rather than to p") as the partial density. 

Furthermore pG is the Radon-Nikodym derivative of the total mass measure p with respect 
to A'^. Analogously to p", we define 

and refer to this p as the density, from now on. 



Now, (3.1.9) implies 

V 

p="^p", a.e. w.r.t. A"*, and for ah t G (0, T). 

a=l 

The mass concentration of constituent a at time t is defined as 

0° 

c" := 

P 



(3.1.13) 



(3.1.14) 



and it follows naturally from (3.1.13) and (3.1.14) that 



a=l 



a.e. w.r.t. p, and for all t G (0, T). 
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Remark 3.1.5. The mass concentration c° is only defined if p > 0. For any t E (0,r), the 
subset 

:= {x G n\p{t,x) = 0} cn 

is a null set w.r.t. the measure p. This can easily be seen: 

H{t,n"t) = J p{t,x)dX^{x) = J OdX'^ix) = 0. 

This implies that is defined almost everywhere w.r.t. /x. We will see later, that defining 
c° actually only is useful in regions where p > is satisfied (i.e. where mass is present), and 
that there is thus no serious problem here. 

Remark 3.1.6. It is obvious that — = — — - = — holds /u-almost everywhere. We thus 

p pe p 

could have defined c° as /5"/ p, and would have obtained the same properties almost everywhere 
w.r.t. p. 

3.2 Kinematics 

The kinematics of component a of the mixture is dictated by a motion mapping x", such 
that 

denotes the position at time t of an a-particle, which was initially (i.e. at t = 0) situated in 

By assumption x" is smooth, and the inverse mapping exists, such that 

X" = H"(t,a;). 

Also by assumption, H° is smooth. Such smooth motion mapping x") with smooth inverse, 
is called a dijfeomorphism. 

For each component a its velocity at position x G il. at time t G (0, T) is denoted by v"{t, x) 
and can be derived from the motion mapping x": 

v-{t,x):=^^:=-X^{t,X^). (3.2.1) 



Remark 3.2.1. The right-hand side in (3.2.1) is indeed a function of {t,x), since we have 



to read X°' as X°' = E"{t,x). The derivative d/dt is taken however only with respect to the 
first variable of x"- 

Definition 3.2.2 (Time derivatives). Let g denote some (physical) scalar quantity associated 
with component a and depending on space and time. We write g = g{t,X°') if we consider 
the quantity from a Langrangian point of view, and g = g{t, x) from a Eulerian point of view. 
The time derivative of g is now defined as 
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In the above definition, note that g[t,X'^) = g(t,x"{t,X°')^. 

Now, we define the barycentric velocity (i.e. the velocity of the centre of mass of a volume 
element) of the mixture by 



= ^c"^;". (3.2.2) 



V := 

P a=l a=\ 

The barycentric velocity is also called mean velocity^ or just velocity of the mixture. 



Remark 3.2.3. A similar expression as in Definition 3.2.2 holds for time derivatives of scalar 
functions g associated to the total mixture: 

Dg dg{t,X) dg{t,x) 

3.3 Balance of mass 

We assume that there is no mass exchange between the components of the mixture. The 
conservation of mass implies that for any VL' G B{Q), for all t £ (0,T) and for all a G 
{1,2,. 

j^fi''{t,x''{t,n'))=0. (3.3.1) 
Note that x"(^) ^') is the configuration at time t of all particles of the component a that were 



initially located in 0'. For ;u"(t, x°'{'t, ^')) in (3.3.1 ) to be well-defined, we need that x"(^) ^') 



is an element of B{^}) for each t € (0,T) and for all fl' G B{n). 

Lemma 3.3.1. If for each t £ (0, T), the mapping x"(^i ■) : ^ ^ is a diffeomorphism, then 
X"(t,f^') e B{n), for each t G (0, T) and for all Q' e B{n). 

Proof. Given that for any t G (0, T) fixed, x"(^) ') is a diffeomorphism from Q to Q, we have 

in particular that its inverse H"(t, •) is a continuous mapping from Q to 0,. By definition of 

continuous functions (cf. [l6], p. 8), (H")~^(t, Jl') is an open set, for any ^7' C O open. 

The Borel cr-algebra is defined as the smallest cr-algebra containing all open subsets of $7. It 

follows that we have that (H")~-'^(i, J7') G B{Q) (i.e. it is a measurable set). By definition of 

measurable functions (cf. [l6], p. 8), we now have that H"(t, •) is measurable. 

Theorem 1.12 (b) ([S], p. 13) provides that {E°')-^{t,n') e B{Q) not only for all W C n 

open, but also for all 0,' G B{Q). That is, is a Borel function (cf. [l6], p. 12). 

H"(i, •) is the inverse of x"(^) ■)> so the above statement also states that for all t G (0, T) fixed 

X"{t, n') E B{n) for ah fi' e B{n). □ 

We now present a few calculations involving Reynolds' transport theorem, that are based 
on [39] (pp. 18-21). We use the notation det(Vx") for the determinant of the Jacobian matrix 
of the motion mapping. The following identity is derived e.g. in [39] (p. 20): 

i- det(Vx°) = V ■v'^ det(Vx"). 
dt 



24 Modelling Crowd Dynamics 



Technische Universiteit Eindhoven University of Technology 



From the conservation of mass statement, formulated in (3.3.1), we derive 
= |/i-(t,x"(t,f^')) = Jt J P''it,x)d\\x) 

d 

~ di 



p"{t,x''{t,X))det{Vxn{t,X)dX\X) 
l^detiVxldX'' 

+ 1 Vp'^-^detiVx^dX" 
n' 

+ 1 p"|det(Vx")dA'^ 

y + Vp" • z;" + /9"V • v"^ det(Vx°)dA' 



n 



^+V.{p-v-)'jdet{VxndX'' 

^+V.{p»v")yX^ (3.3.2) 

Since t and fi' can be chosen arbitrarily, we conclude from the expression in the seventh line 
of ( [3X2I ) that for all a G {1, 2, . . . , i^} and for ah t € (0, T) 

+ V • (p"-^")) det(Vx") = 0, A'^-a.e. in (3.3.3) 

Since the inverse of x" is assumed to be differentiable, we know that det(Vx") 7^ (cf. 
remark on p. 4 of [39]). Thus, we deduce from (3.3.3) that 

^ + V • = 0, A'^-a.e. in Q. (3.3.4) 

ot ' 



The weak formulation of (3.3.4) leads to 



^ / V"p"dA'^ = I v'^- Vip^p'^dX'^, (3.3.5) 



n n 

for ah a G {1, 2, ... , u}, for ah t G (0, T), and for ah V" G C^{^)- 



Remark 3.3.2. From (3.3.4) we can derive the local mass balance for the mixture as a whole. 
To achieve this, we sum (3.3.4) over all a and take (3.1.13) and (3.2.2) into consideration. 
This yields that for all t G (0, T) 

dp 



dt 



+ V-{pv) =0, 



A -a.e. in 0. 



(3.3.6) 
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The weak formulation of (3.3.6) can be either deduced directly from (3.3.6), or derived from 
(3.3.5). For the latter way, we take the same test function ip°' = ip £ Cq for each a, again 
sum over all a and use (3.1.13) and (3.2.2). This procedure results in 



dt 



(3.3.7) 



In Remark 3.1.5 we indicated that defining c"^ is actually only useful in regions where p > 0. 



This is made clear by (3.3.6) and (3.3.7). The mass concentrations c"' are incorporated in the 
definition of v. However v only appears in combination with p as the product pv. If p = the 
product is zero any way, and thus (at least physically) it does not matter how (or whether) 
is defined. 



3.4 Generalization 



We would like to extend the ideas presented in Sections 3.1-3.3, such that they hold not only 
for those measures that are absolutely continuous w.r.t. A*^. Indeed we have seen in Section 



2.1 that in general, a measure might also contain a singular part. The mixture-theoretical 



concepts presented so far do however not hold for singular measures. 

Let p"(t,r2'), like before, denote the mass of constituent a present in Q' E B{Q,) at time 
t G (0,r). To obtain a more general framework, we need to revoke the assumption that 
^ £qj. a E {1, 2, ... , u}. In this setting we thus postulate the following: 

Postulate 3.4.1 (Properties of //"). For all a G {1, 2, . . . , i/} and all t G (0, T) we postulate: 

1. /i° ^ 0. 

2. /x° is a-additive. 



Furthermore, as in (3.1.2), we define the total mass by 

V 

p{t,n') :=^p"(f,f7'), 



(3.4.1) 



a=l 



for all n' G B{n) and all t G (0,r). 

This new setting, without the absolute continuity demand, implies that in most cases of 



Section 3.1 we cannot apply the Radon-Nikodym Theorem anymore. However, note that due 
to ( |3.4.1D and Part [T] of Postulate [sXT we have that /i° <C p for all a G {1, ... , v}. Thus for 
each /i" a unique Radon-Nikodym derivative exists with respect to p. Let us denote 



dp ' 



(3.4.2) 



which is also called the mass concentration of constituent a at time t. 
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Remark 3.4.2. If we still would have been in the absolute continuous case of Section [3?H both 
jj T and T fi were true. Due to Part [3] of Lemma 2.3.1 this implies that 1/p = dr/dfi, 
which is defined ^-almost everywhere. Part [T] of the same lemma then implies that 



dfj, 



dfi'' dr 
dr dfj, 



-a 



Also: 



p 



p pQ 

and since all expressions above are defined p-a.e., we have that 



P 



;U- almost everywhere. 



Hence, in the absolutely continuous case (3.4.2) is equivalent to the definition of c" given in 
Section 13.11 



It follows from (3.4.1) that: 



a=l 



a.e. w.r.t. p, and for all t G (0, T). 



(3.4.3) 



Up to the definition of the barycentric velocity, in Section 3.2 we have not used the absolute 
continuity of any of the mass measures. The ideas presented there are also applicable in the 
generalized setting. We only have to redefine (in a very natural way) the barycentric velocity 
of the mixture by 



(3.4.4) 



a=l 



From a notational point of view it is not difficult to generalize the weak formulation of the 
balance of mass concept, as given in (3.3.5). Since we know that />" = dp'^/dX'^, the weak 
form can be written as 

d 



ft I 



V^l^'^dp" 



(3.4.5) 



for aU a G {1, 2, . . . , i/}, for ah t G (0,r), and for ah V'" e C'd(^)- This 'trick' is merely 



a matter of notation, and the above is equivalent to (3.3.5) for absolutely continuous mass 



measures. However, we still need to make sure that this expression also makes sense for more 
general mass measures. 



We are running ahead by mentioning this now, but we will see later (see e.g. Section 4.8) 
that we are mainly interested in measures that are either absolutely continuous, or discrete 
(i.e. sum of Dirac distributions), or possibly a combination of the two. This means that we 



explicitly exclude the singular continuous part from the measure, and only show that (3.4.5) 



makes sense for discrete measures. Note that the measures we do allow, are (combinations 



of) exactly those measures introduced in Section 2.4 



We consider /i" to be a single Dirac measure centered at a time-dependent position x[t) G 
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for each t G (0,T), and assume its motion to be described by dx{t)/dt = v°^{t,x{t)). We 
choose an arbitrary function -0" E C'o(^) and take the inner product with this function (eval- 
uated in x{t)) on both sides of the equality. In the resulting left-hand side, we recognize the 
chain rule (recall that tp"^ is once continuously differentiable) , so 



Now, we have obtained 



--r{x{t)) = Vi:^{x{t))-v^{t,x{t)). 
at 



(3.4.6) 



We apply the definition of the Dirac measure, and since /Li" = 5a,(t), we find that (3.4.6) can 
also be written as 

d 



dt 



iP"dfj,° 



To extend this approach, now let /x" be a linear combination of Dirac deltas centered at Xi{t), 
and let each of these move according to 



d 
dt 



Xi{t)=v"{t,Xi{t)). 



(3.4.7) 



Here the index i is taken from a countable, possibly infinite, index set J . If J is infinite, the 
coefficients of the Dirac deltas must have finite sum. Similar arguments as before yield that, 
for each i £ J', we have 



dt 



^"(x,(t)) = VV^"(xi(t))-?;"(t,x,(t)). 



If fi" = X^jgj- ai6x^(t) with nonnegative coefficients such that Xliej' ^ then we have 



^ a,-r{xi{t)) = aiVrixiit)) ■ v'^it, Xi{t)), 



iej 



which is, by the definition of the Dirac measure, again equal to (3.4.5). 



Note that the weak formulation (3.4.5) also makes sense for linear combinations of abso- 
lutely continuous and discrete measures. For the absolutely continuous part we have the 
balance of mass-interpretation; for the discrete part we interpret it via (3.4.7). 



Remark 3.4.3. We have seen in Lemma 2.1.6 that any discrete measure can be written as 
a linear combination of Dirac distributions. However, if (for every t) ^"{t, •) is discrete, not 
only the centres of the Dirac deltas, but also the coefficients might be time-dependent. It 
will be shown later that the time-dependent discrete measures we work with have constant 
coefficients (see Corollary 4.5.8). This means that the time-dependence is only present in the 
location of the Dirac deltas. It will also be shown that these positions Xi{t) satisfy (3.4.7); cf. 
Lemma 14.5.91 

In this spirit, (3.4.5) makes sense for the types of measures that are relevant for us. 
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Remark 3.4.4. Analogously to Section 3.3, we can deduce from (3.4.5) the structure of a 



weak formulation that applies to the total mass measure /Lt. We take the same test function 
V'" = tp ^ C'o(^) for each a. Each integral w.r.t. /x" is transformed into an integral w.r.t 
using (3.4.2). Consequently, summing over all a, and using (3.4.3) and (3.4.4), yields 



d 
dt 



V ■ VTpdfi. 



(3.4. 



3.5 Entropy inequality 



In this section, we place the mixture concepts discussed in Sections 3J^ 3.4| into a thermody- 
namical context. Our aim is to derive an entropy inequality. This inequality is strongly related 
to concepts like the second axiom of thermodynamics, or the Clausius-Duhem Inequality. 
For more information on thermodynamics and on the role of entropy, the reader is referred 
to e.g. [52] or p7]. 



To each constituent of the mixture we assign a function ry" : (0, T) x — )■ M, which is 
called the entropy density (per unit mass) of component a. The entropy density for the total 
mixture is defined as 



(3.5.1) 



a=l 



Note that this definition has the same structure as the definition of the barycentric velocity 



in (3.2.2). The entropy Sn'{t) assigned to fi' C at time t is defined by integration of r] with 
respect to the density p, that is 

Sn'it) := [ ripdX". 



(3.5.2) 

To each component we also assign a temperature T° : (0,T) x — t- M^. 

Following the ideas of [42j (p. 7), [20j (p. 347), and p4j (p. 28), we now postulate the following: 
Postulate 3.5.1 (Local entropy inequality). Locally the following inequality holds: 



Dt 



(3.5.3) 



a=l 



Here jr) is the entropy flux vector, and r° is the volume supply of external heat to constituent 



a. 



Remark 3.5.2. In the term Yla=i P^^f^/T", which is the volume supply of external entropy 
to constituent a, we recognize a similar structure as in (3.2.2) and (3.5.1). It can therefore 
be considered as some "barycentric quantity". 



Remark 3.5.3. Sometimes (3.5.3) is written as 

d 



dt 



(prj) + V • {p7]v + jr,) -J2p"^^^^ 



Q = l 
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cf. e.g. [4.2] , p. 7. Using the local balance of mass (3.3.6) and Definition 3.2.2, we can indeed 
show that these two formulations are equivalent: 



dt 



{pr]) + V • {pit]v) 



dt 
D 



(prj) + V{pri) ■ V + prjS/ ■ v 
(prj) + pr]V ■ V 



Dt 

Dp Drj , , ^ X 

:V + Pj^ + ?/(V • [pv] - Vp-v) 



Dt 



Drj ,Dp ^ ^ , , 



Dt] 
Dr] 



dp 
,dp 



P-F^ + + »?V • (pv) 



V • {pv)) 



p 



Dt] 



(3.5.4) 



For an arbitrarily fixed i}' C (such that dQ' is sufficiently regular) we now calculate the 
time derivative (at fixed time t E (0, T)) of the total entropy contained in 0,': 



d_ 
dt 



rjpdX" 



d_ 
dt 



{7^p)dX^ 



E 

a=l 



/ 



a=l 



{prjv + jri) ■ ndX 



d-l 



(3.5.5) 



an' 



where we have used the local entropy inequality in the second line of (3.5.5). 



We call (3.5.5) the global entropy inequality. 



Remark 3.5.4. It would have been possible to derive the entropy inequality for the whole 
mixture from the partial entropy inequalities formulated individually, per constituent; cf. e.g. 
|25j (pp. 475-476). However, nowadays the general consensus is that this approach gives a 
too restrictive result: the motion of the mixture is overconstraint. For more details on this 
fundamental issue, the reader is referred to |32] (pp. 11.12-13), and [8j (p. 865). 

Remark 3.5.5. In general we are also interested in mass measures that are not exclusively 
absolutely continuous. This naturally leads to the question whether we can generalize the 
concept of entropy and the accompanying entropy inequality. However, at this point, we 
should ask ourselves what the physical meaning of such generalization is. More understand- 
ing is needed in order to judge the physical relevance of an extension to a broader class of 
measures. There are also practical objections. For example, caution is needed when general- 



izing boundary terms. It is all but obvious in what way the integral over dVL' in (3.5.5) should 
be defined for a general mass measure, that may also contain a discrete part. 



In Section 4.3 we introduce an explicit entropy density, and derive the corresponding entropy 



inequality. Despite of the aforementioned difficulties, we show afterwards (see Section 4.3.3) 
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that structurahy the same inequahty can be obtained for discrete measures (compare the 
statement of Theorem 4.3.3 with ( 4.3.30| )). 
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Chapter 4 

Application to crowd dynamics 



In this section we explain and extend the measure-theoretical approach, developed by Bene- 
detto Piccoli et al. (cf. e.g. [441 1451 [T7] ). We intend to fit some of their ideas to our framework 
presented in Section [3} 

We consider a population located in a given domain 17 C M*^. To capture physically real- 
istic situations we take d £ {1,2,3}. Let T G (0,c«) be the fixed final time of the process. 
We define a measure n : (0,T) x 13{Q) — )• M"*", such that fi{t,Q') represents the mass of the 
part of the population present in a region il' G B{i^) at time t E (0,r). 

We assume that the population consists of a fixed number of subpopulations (these were called 
constituents in Section |3]), indexed a G {l,2,...,i^}. The mass of each subpopulation present 
in Q' E B{Q,) at time t E (0, T) is given by the time-dependent mass measure n"{t,Q'). 
The mass measures /i and /i" are related via fi = Yla=i explicitly assume that fJ-°'{t, •) 

is a finite measure for all t E (0, T) and a E {1, 2, . . . , u}. As a result, fj,{t, •) becomes a finite 
measure for all t E (0,T). 



4.1 Weak formulation 



Each subpopulation moves according to its own motion mapping x", from which a velocity 
field v°'{t, x) follows. Note that for each a, the dependence of on t indicates the functional 



dependence on all time-dependent measures /i", as we will see in Section 4.2 



Remark 4.1.1. We henceforth disregard the assumptions with respect to the motion map- 
pings and velocity fields v°', which were done in the derivation of the balance of mass 



equations presented in Section |3j The result, Equation (3.4.5), is taken here as a starting 
point for further modelling, without considering the underlying assumptions. 

Recapitulating: the fact that here we deal with u time-dependent mass measures /i", 
transported with corresponding velocities v", translates into 



9^° 
~dt 



+ V • ifi^v' 



0, 



for ah a E {1, 2, ... , i^}. 



(4.1.1) 



Equations (4.1.1) are accompanied by the following initial conditions: 



/^"(Or 



Mo, 



for each a E {1, 2, ... , u}, 



(4.1.2) 
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for given ^q. 

This partial differential equation in terms of measures is a shorthand notation for the weak 



formulation presented in Sections 



3.3 



and 



3.4 



Namely, for all test functions V" ^ Cq{CI), 
where a £ {1, 2, . . . , z/} and for almost every t £ (0, T) the following holds: 



d_ 
dt 



V'"(x)(i/x°(t,x) 



a £ {1,2,. 



(4.1.3) 



Remark 4.1.2. If, for some a G {1,2, . . . the mapping t i— )• J^il:"^{x)d^"{t,x) is abso- 
lutely continuous on (0, T), for any choice of if)" £ Cq{Q.), then it is differentiable at almost 
every t £ (0,T). See e.g. Theorem 7.20 in |46j . p. 148. This means that the left-hand side of 
(4.1.3) exists for almost every t. 

Remark 4.1.3. I f we de mand that £ ((0, T); L^c (fi)) for ah a G {1,2, . . . then the 
right-hand side of (4.1.3) is well-defined. Indeed, since ip" £ Cg(il), we have that ||VV'"||lcx)(q-) 
is finite for each a £ {1, 2, . . . , i'}. We thus have for each a £ {1, 2, . . . , i/} that 



T 



n 



T 



dt 



n 

T 



llVniL-(Q) / \v''{t,x)W^{t,x)dt 



l|VniLoo(c)ll^"(t,-)|lL^(Q)rfi 



< oo. 



lL°°(n) 



dt 



(4.1.4) 



In particular, it follows that J^v"'{t,x) ■ Vxp°'{x)dn°'{t,x) is finite for almost every t £ (0,T) 
and for all a £ {1, 2, . . . , z^}, and thus the right-hand side of (4.1.3) is well-defined. 



Definition 4.1.4 (Weak solution of (4.1.1)). The vector of time- dependent measures 



where each ^u" satisfies the initial condition (4- 1-2), is called a weak solution of (4-1 -D, if for 
all a £ {1,2, ... the following properties hold: 

(i) for all -0" £ Cq{^), the mappings 1 1— )• J^ilj"'{x)d^°'{t,x) are absolutely continuous, 

(ii) £L\{Q,T);Ll^{n)), 



(Hi) (4- 1.3) is satisfied. 
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We deduced a weak formulation with respect to the total mass measure fj, in Section 3^ 
We use the shorthand notation 

|^+V.(H = 0, (4.1.5) 

accompanied, for given ^q, by the initial condition 

/i(0,-) = ^o. (4.1.6) 

This shorthand notation should, again, be interpreted as (cf. ( |3.4.8 )) 

d f f 

i;{x)d^i{t,x) = v{t,x) ■Vip{x)dfi{t,x), (4.1.7) 



dt 

for all ijj G Cq{Ci) and almost every t G (0,r). 



A weak solution to (4.1.5) will be understood in this section always in the sense of Defi- 
nition I4.1.5t 



Definition 4.1.5 (Weak solution of (4.1.5)). The time- dependent measure {^)t^o, satisfying 
the initial condition is called a weak solution of {4-1-5), if the following statements 

hold: 

(i) for all ip G Cq{Q), the mapping 1 1— )• f^Tp{x)dfi{t, x) is absolutely continuous, 

(ii) vGLmo,Ty,Llin)), 

(Hi) (x^. j. ?| ) is fulfilled. 

The main difference between our approach here and the one presented in ^7\, is that we 
take the freedom to allow each subpopulation to have its own velocity field f". Note that 
the crucial modelling steps take place when we decide what v" actually looks like. We then 
characterize the motion of the corresponding subpopulation. In particular, we define the way 
in which this motion is influenced by the crowd surrounding it (belonging to any of the com- 
ponents). 

In [T7j merely the setting of transporting the total mass measure is considered (see Defi- 
nition 4.1.5). This means that only the barycentric velocity v of the crowd seen as a whole 
can be prescribed. However, we thus lose the ability of modelling the interaction between 
subpopulations of different types. Distinct behaviour of subpopulations is namely driven by 
underlying partial velocities of distinct nature. 

To be more specific, |17j considers a combination of an absolutely continuous measure (a 
'cloud' of people, in which individuals are indistinguishable) and a sum of Dirac measures 
(point masses). Only v is prescribed in the framework of that paper, partial velocities being 
not included. This implies that on the macroscopic scale (the crowd) the absolutely continu- 
ous and discrete parts essentially move according to the same velocity field. Correspondingly, 
the Dirac masses cannot evolve independently from the cloud; they are in fact part of the 
cloud with some "pointer" attached to them. 

Note that this situation is incorporated in our model as a special case, where we take one /x° to 
be a combination of Dirac measures and an absolutely continuous part. However, we feel that 
it is of no use drawing attention to these point masses, if they cannot behave independently. 
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4.2 Further specification of the velocity fields 

Until now, we have not explicitly defined the velocity fields v" (a G {1, 2, . . . , i^}). In Section 



4.1 we have already remarked that the crucial modelling step takes place at the moment when 
we decide on the explicit form of v". By the choices we make here, we define the character- 
istics of a subpopulation. That is, we assign, so to say, a personality to the members of that 
subpopulation. 

In the definition of v" we can incorporate the way in which an individual is infiuenced by the 
people around him. This individual might be shy, wanting to keep distance from others. The 
individual's character might also be quite the opposite, driving him to come close to other 
people. We even have the freedom to distinguish between the way the individual responds 
to the presence of others, based on the subpopulation that other individual belongs to. For 
example, one subpopulation might consist of acquaintances (triggering attractive behaviour), 
while a second subpopulation consists of 'enemies' (from which the individual is repelled). 

In this section, we show the way we want to model the velocity fields. Very much inspired by 
the social force model by Dirk Helbing et al. [31], the velocity of a pedestrian is modelled as 
a desired velocity v^^^ perturbed by a component v"^^. The latter component, which is called 
social velocity, is due to the presence of other individuals, both from the pedestrian's own 
subpopulation and from the other subpopulation. This effect is modelled by functional de- 
pendence of v^oc on the set of measures { /i^ , , . . . , /u"^ } . The desired velocity is independent 
of these mass measures, and represents the velocity a pedestrian would have had in absence 
of other people; this implies that v^^g is independent of t, provided the environment of an 



individual does not change in time. Cf. Remark 4.2.1 



The velocity is thus for a G {1, 2, . . . , i/}} defined by superposing v^^^ and its pertur- 
bation , 



^SOC" 



v°'{t,x) ■.= v2e^{x)+v^^^{t,x), for ah t G (0, T) and X G 17. (4.2.1) 

The component v"^^ models the effect of interactions with other pedestrians on the current 
velocity. This is essentially a nonlocal contribution. Since the effects (in general) differ from 
one subpopulation to the other, we assume that v"^^ has the form: 



(t,x):=^ I /^(|y-x|)5(0°^)^— ^d/(t,2/), for ah a G{1, 2,..., z.}. (4.2.2) 

^=^n\{x} 

Remark 4.2.1. Note thus that the t-dependence of v" and v"^^ is not an explicit time- 



dependence. As (4.2.2) shows, the social velocity depends functionally on the time-dependent 
mass measures fJF. Although we have not done so, we are also allowed to choose v^es = 
Vdesit,x). This ought to be an explicit dependence on t due to the interpretation given to 
Udes- the velocity a pedestrian would have if he was the only person in the room. 
Obviously, this term should depend only on the geometry of the domain (and not on the 
crowd it contains). If v^es depends explicitly on the time, this refiects that the domain is 
non-static. For example, such situation occurs when certain areas are not always accessible 
(opening hours), or are temporarily obstructed. For the moment this is too farfetched. 



In (4.2.2), we have used the following modelling ansatze: 
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is a function from M4. to M, describing the effect of the mutual distance between 
individuals on their interaction. The subscript and superscript denote that this specific 
function incorporates the infiuence of the members of subpopulation /? on subpopulation 
a. In principle we distinguish between two kinds of effects: 

— Attraction-repulsion: this is the interaction between acquaintances (and similar 
kinds of interaction). In this case, fp is a composition of two effects: at short 
distance individuals are repelled, since they want to avoid collisions and congestion, 
but if their mutual distance increases they are attracted to other group mates, in 
order not to get separated from the group. 

— Repulsion: i.e. 'shy' behaviour. Here, contains only the repulsive part, since 
the individual in subpopulation a does not like to come close to subpopulation /3. 



Graphical representations of these two choices of are depicted in Figure |4.1 



9^y denotes the angle between y — x and i^2^es(^)- angle under which x sees y if it 
were moving in the direction of v^^^^x). 

g is a function from [— vr, vr] to [0, 1] that encodes the fact that an individual's perception 
is not equal in all directions. We choose: 

g(e) :=a + {l- ^)1±^^, for 9 G [-7r,7r]. (4.2.3) 

This definition ensures that an individual experiences the strongest infiuence from some- 
one straight ahead, since g{0) = 1 for any a G [0, 1]. The constant a is a tuning param- 
eter, called potential of anisotropy. It determines how strongly a pedestrian is focussed 
on what happens in front of him, and how large the infiuence is of people at his sides 



or behind him. The effect of the choice of this parameter can be seen in Figure 4.2 
Note that cos{9"y) can be calculated in a convenient way: 



cos(6': 



^ (y - x) ■ v^csi^) 
"^y' \y - x\\v2^^[x)\' 



Remark 4.2.2. Other choices for v^oc s-re also possible. Based on the current velocities we 
can make an individual anticipate the distance he expects to be from another pedestrian in 
the (near) future. Although this is probably more realistic, it increases the complexity of 
the model dramatically. If these changes are made, namely depends on the velocities 
(/3 E {1, 2, . . . , z^}); in particular on itself. The definition of v'^ becomes thus implicit, and 
is much harder to work with. 



4.3 Derivation of an entropy inequality 



In this section, we derive an entropy inequality concept in the spirit of (3.5.5) in Section 3.5 



Here, we also specify under which conditions this inequality holds. For simplicity and clarity, 
calculations are performed for a crowd without distinct subpopulations first; see Section 4.3.1[ 
A multi-component crowd is considered afterwards in Section 4.3.2 The inspiration for the 



one-population case was provided by pU] (but we slightly deviate from their definitions). 
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Rr 



Figure 4.1: Graphical representation of typical examples of the function attraction- 
repulsion (left) and repulsion only (right). Here, Rr is the radius of repulsion: an individual 
is repelled if its distance to someone else is smaller than Rr, while Ra is the radius of attraction: 
an individual is attracted to an 'acquaintance' if their mutual distance is between Rr and Ra- 
Note that Rr might be chosen differently in either of the two choices. Typically Rr will be 
larger in the repulsion case than in the attraction-repulsion case. 








Figure 4.2: Graphical representation of the function g. We have ^(0) = 1 always, and 
g{—TT) = 5(vr) = a. Furthermore g is symmetric around 9 = and increasing on (— vr, 0) (thus 
decreasing on (0,7r)). The plot has been made for a = 0, 0.3, 0.5, 0.7, 0.9, 1 (bottom to top). 



We work with a mass measure that is absolutely continuous w.r.t. the Lebesgue measure, 
that is: we do have a mass density. 

If one wants to deviate from the setting of [TB], an alternative approach is suggested in 
Appendix [D] There, the entropy density follows, if the so-called free energy F is explicitly 
provided. The crucial decision is thus the choice of this free energy. In Appendix |D| this 
is done for the example of the ideal gas. If one would find a good choice of F, the same 
arguments could be used for crowds. 



4.3.1 One population 

We consider the presence of a single population in the domain Q C M'^. The time-dependent 
density is p : (0, T) x — )• M'^. We recall that it satisfies the balance of mass equation (3.3.6): 

^ + V-{pv)=0, a.e. in 0. (4.3.1) 

Assumption 4.3.1. We assume that the velocity field u : (0, T) x — )- is of the form 

V :=VV + VWi^p, (4.3.2) 
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where y : O — )• M is called confinement potential, and W : M*^ 
potential. We also assume that W is symmetric, that is 



is called interaction 



wio = w{-o, 



for all ^ € 



Note that in |16j a minus-sign is added to the definition in (4.3.2) 
The convolution VW * p is defined as 



{VW*p){t,x) := J VW{x-y)p{t,y)dX'{y), 
n 



for all t G (0, T), and x G 0, (4.3.3) 



and thus this part is time-dependent via p. For brevity we often write p{x) instead of p{t, x) 
in the sequel. The time-dependency is to be understood implicitly. Analogously, if we write 
(Vl^*p)(j;), we actually mean {VW -k p){t,x). 

Remark 4.3.2. The component Vy is the desired velocity fdes of Section [4. 2[ Similarly, the 
component VVF* p is the social velocity v^oc- 



We already proposed a velocity field of the following form (cf. Section 4.2): 
v{t,x) ■.= VAes{x)+ / f{\y - x\)g{e^y)^ — —p{t,y)dX^{y). 



(4.3.4) 



We can take O as our domain of integration here, instead of J7 \ {x} (cf. (4.2.2)). This is be- 
cause {x} is a A'^-negligible set; the integral has the same value with or without the exclusion 
of this set from the domain. This holds for absolutely continuous measures, but not for any 
measure in general. 



The velocity as defined in (4.3.4) fits to the structure of Assumption 4.3.1, if: 



• 'I'des can be written as W for some potential V . Note that in [441 a comparable situation 
is covered. The potential is found by solving the Laplace equation Ay = in the domain 
$7. For us this would mean that ?;des is divergence-free. However, the difference with 
our situation is that HI] normalizes and rescales W afterwards. 



we disregard here the factor g{6xy), that is, we set g = 1 (or a = 1). This is necessary at 
this stage, because the inclusion of g makes it impossible to find a symmetric interaction 
potential 



there is a function F 

y -X 



^ such that / 



\y - x\ 



-f{\x-y\) 



x-y 
\x - y\ 



F'. Then: 

x-y 



F'{\x-y\) 



\x - y\ 



VFi\x-y\), (4.3.5) 



where the last step of the calculations is due to the chain rule, and the fact that V|x — 
X — y 

y\ = r. Now write W(x — y) := Fi\x — y|); note that this indeed implies that 

\x — y\ 

W{i) = W{-i) for all i £ R'^. 



^This is an inevitable, but disappointing choice. We thus recommend further research in this direction. 
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It is, for example, possible to have v^es = a G M'^; that is, the desired velocity has constant 
magnitude and direction. To comply with the assumption that Vdes can be written as W, 
we have to take V{x) = a ■ x. 

It follows from the aforementioned assumptions on g and / (and the subsequent calcula- 
tions) that Vsoc = VVF -k p. In Figure 4.3 an impression is given of the functions F that 
correspond (via the relation / = —F') to those functions / plotted in Figure 4.3 These 
functions are unique up to additional constants. 




Figure 4.3: Graphical impression of the function F corresponding, via the relation / = —F', 



to the functions plotted in Figure 4.1 attraction-repulsion (left) and repulsion only (right). 
As in [16], we define the entropy density rj as 



r]{t,x) ■.= Vix) + -{W*p){t,x). 



(4.3.6) 



The corresponding entropy of the system in Q, is then, according to (3.5.2), given by: 



S{t) = J r]{t,x)p{t,x)dX'{x). 
n 



(4.3.7) 



The central question here is: Does the choice (4.3.6) - (4.3.7) satisfy the Clausius-Duhem 



Inequality? We answer this question in Theorem 4.3.3 



Theorem 4.3.3 (Entropy inequality). Assume that v satisfies Assumption \4.3A . that the 
entropy S is given by (4-3.6) - (4.3.1), and also that the system is isolated^ 
Then the following inequality holds: 

dS 

^ ^ 0. 

dt 

Proof. We investigate directly the time-derivative of the entropy S", i.e. we have: 

Tt= j + ^^{W^p){x)p{x)d\\x) + \ j iW ^ p){x)^-^^{x)d\\x) 

n n n 

V{x)^-f^{x)d\\x) + \j{wJ^^){x)p{x)d\\x) + \ j{W^p){xf-^^{x)d\\x). 
n n 

(4.3.8) 



■^The statement that the system is isolated means that there is no flux through the boundary of fi, i.e. 
pv ■ n — at d^l. 
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Note that 



(W^*^)(x)p(x)dA^(x) 



W{x-y)^^{y)dX''{y)p{x)dX''{x) 



n n 



n n 



n n 



n n 



W{x-y)^{y)p(x)d\\y)d\\x) 



W{y-x)^^{x)p{y)dX\x)dX%) 
W{x-y)p{y)d\%)^^{x)d\''{x) 



{W*p){x)^{x)dX^{x). 



(4.3.9) 



We have replaced x by y, and vice versa, to obtain the third equality (this is solely a matter 
of notation). To obtain the fourth one, we used that for all x,y £ Q, W{y — x) = W{x — y) 
and interchanged the order of integration. 



From (4.3.9), we conclude that (4.3.8) can be written as 



d^ 



V{x)^{x)dX''{x)+ I {W^p){x)^{x)dX\x) 



^V{x) + {W*p){x))^{x)dX''{x). 



(4.3.10) 



Substitution of the balance of mass (4.3.1) in (4.3.10) yields 



d5 



{V{x) + (VF*p)(x))V • {p{x)v{x))dX'^{x) 



{V{x) + {Wi^ p){x))p{x)v{x) • ndX'^-^ix) 



dn 



+ J V{V{x) + {W*p){x)) ■ p{x)v{x)dX'^{x). 
n 



(4.3.11) 



By the hypothesis that pu-n = at the boundary dft, the boundary term in (4.3.11 ) vanishes. 
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We proceed as follows: 
dS 



V{V{x) + {Wi.p){x)) ■ p{x)v{x)dX'^{x) 

{VV{x) + V{W * p){x)) ■ p{x)v{x)dX'^{x) 

{VV{x) + (VVF*/9)(x)) • p{x)v{x)d\'^{x) 

v{x) ■ p{x)v{x)dX'^{x) 

\v{x)fp{x)dX'^{x) 
^ 0. (4.3.12) 



n 



Note that the relation 



ViW*p){x) =V,J W{x-y)p{y)dX'{y) 
V^W{x - y)p{y)d\\y) 
: j{VW){x-y)p{y)d\\y) 



n 

--{VWi.p){x) 



was used in the third step. 



In (4.3.12), we recognize the entropy inequality we were looking for: 

^ ^ 0. (4.3.13) 

dt ^ ' 

□ 

Remark 4.3.4. In Theorem |4.3.3 and in its proof we have implicitly assumed a sufficient 
amount of regularity of the boundary dQ. This is also important for the steps we are about 
to take. 

Let us comment on the situation in which pv ■ n = does not necessarily hold at 50. 
This means that we allow mass (and consequently also entropy) to escape from or enter the 
domain of our focus. Instead of (4.3.12), this would yield 

^ ^ - y {V{x) + {W* p){x))p{x)v{x) ■ ndX'^-^x) 
an 

rj{x)p{x)v{x) ■ ndX'^-^{x) - j ^{W p)ix)p{x)v{x) ■ ndX'^-^{x). (4.3.14) 
an an 
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After defining the entropy flux jn as 



we are thus in the setting of (3.5.5) if there is no external volume supply of heat: 

{vpv+jr,) ■ndX'^-\ 



dS 



dt 



an 



4.3.2 Multi-component crowd 



In this section, we extend the results of Section 4.3.1 to the situation in which the crowd 
consists of a given number of subpopulations. Let the crowd in the domain Q C M*^ consist 
of v subpopulations, indexed by a G {1, 2, ... , u}. The time-dependent density of component 
a is denoted by : (0, T) x J7 M+. Cf. the definition ( |3.1.12[ ) of as a Radon-Nikodym 
derivative in Section [3Tl 



Assumption 4.3.5. 1. The velocity field of component a, denoted by u " : (0, T) x — ?■ M*^ 

is assumed to be of the form 



:= + VWj^ ★ /, for all a G {1, 2, ... , i/}, 



(4.3.15) 



/3=1 



where V"' : — )• M is called the confinement potential of component a, and : — t- M 
is called the interaction potential of component f3 (affecting a). 

2. For each a, /3 G {1, 2, . . . , i^} we assume that W^{0 = W^(-^) holds for all ^ G R'^. 

3. We assume symmetric interactions, that is 

W^ = W^, for alia, P e {1,2,..., u]. (4.3.16) 



Remark 4.3.6. Part [3] of Assumption 4.3.5| implies for example that we do not allow a 
"predator-prey relation" between two subpopulations. In such relation a predator should be 
attracted to the prey-population, but a prey should be repelled from the predators. These 
interactions are asymmetric. 



The following balance of mass equation is satisfied for each aG{l,2,...,i/}: 



dt 



+ V • = 0, a.e. in 0. 



(4.3.17) 



We recall from Sections 3.1-3.2 the definitions of the total density p and bary centric velocity 



v. 



■.=yp-v'^. 
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On the macroscopic scale, the following balance of mass equation is satisfied: 

dp 



at 



+ V • [pv) = 0, a.e. in Q. 



(4.3.18) 



For all t S (0, T) and x €z Q we define the partial entropy density of component a as 

7?-(t, x) := V"{x) + ^ * 

/3=1 

The entropy density of the whole crowd is defined as 

Consequently, the entropy of the system at time t is given by 
S{t)= I ri{t,x)p{t,x)dX'^{x) 



= E / {v^i^) + lj2^W^^p^){t,x))p^{t,x)dX''{x). (4.3.20) 
In the spirit of Theorem |4.3.3 we can formulate the following theorem: 

Theorem 4.3.7 (Entropy inequality). Assume that for each a we have that satisfies 
Assumption 4-3.5. Moreover, assume that the entropy is given by (4-3.20) and that the system 
is isolated. 

Then the following inequality holds: 



dS 
~dt 



^0. 



Proof. The approach here is very much in the spirit of the proof of Theorem 4.3.3 We thus 
consider the time-derivative of S: 



dS_ 
~dt 



E / V^{xf-^{x)d\\x) 
"=in 

+ ^E/ E|(W^^"*P^)(^K(^)dA'^ 



a=l^ 0=1 



1 



V ,. u 



dp" 
dt 



dt 



E / V'^{x)'^{x)dX''{x) 

+ ^EE/(^^"*^')(-)^(-)^^'(^ 

"=1/3=1 n 



(4.3.21) 
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Following the idea of (4.3.9), we derive 



EE j{W^*^^){x)p-{x)dX 



[X) 



"=1/3=1 n 



iziZl I W^{x-y)^-^{y)d\%)p'^{x)d\\x) 

EE / / W^{x-y)^{y)p^x)dX%)dX%) 
"=1/3=1 n n 

jzt. f f W^{y-x)^{x)pHy)dX'{x)dX''{y) 
/^=i"=i6 l 

t^iZj I Wi{x-y)p^{y)dX%)^{x)dX''{x) 



a=l 0=1 



n n 



a=l 13=1 



EE Uw!*p^)ix)^ix)dX'{x). (4.3.22) 



dt 



We have replaced x by y, a by /3, and vice versa, to obtain the third equality (this is again 
solely a matter of notation). To obtain the fourth one, we used the fact that Wa{y — x) = 
Wa{x — y) for all a, /3 G {1, 2, . . .}. Moreover, we interchanged the order of summation and 
integration. 



We now combine (4.3.21) and (4.3.22), and conclude that 



dS_ 
~dt 



a=l p=l 



) -^{x)dX\x). (4.3.23) 



We substitute (4.3.17), the balance of mass per component, in (4.3.23), by which we obtain 



Itt 



-jZl [^"(^) + E(2(^^" + ^a^) */)(^)] V • [p'^{x)v^{x))dX\x) 

a=l^ /3=1 

-E / [v'^{x) + jZ{\{W^ + W^)^p^){x)y{x)v%x)-ndX''-\x) 

a=l/^ /3=1 

+ E / y\v''{x) + jZ{\{W^ + W^)^p^){x)]- p%x)v'^{x)dX\x). (4.3.24) 



a=l 



/3=1 



Under the assumption that there is no flux of mass through the boundary of for any of the 
constituents a, we have p^^v" • n = at for all a € {1, 2, . . . , z^}. This means that the 



boundary term in (4.3.24) vanishes. Taking also Part ^ of Assumption 4.3.5| (i.e. symmetric 
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interactions) into consideration, (4.3.24) reads 
dS _ 

a=l 



X] 1 |^"(x)|V"(x)ciA'^(x) 



>0. 



(4.3.25) 



By (4.3.25) we thus have the following entropy inequality: 

dS 



dt 



> 0. 



(4.3.26) 

□ 

If we allow mass (and consequently also entropy) to escape or enter through the boundary 
of J^, we should consider the boundary term in ( |4.3.24| . Let us define the entropy flux jjj as 



a=l (3=1 



Still assuming symmetric interactions, and using (4.3.19), (4.3.24) would now read 

){x) p"(a;)v°(x) •ndA'^~^(x) 



dS^ 
~dt 



^r?"(x)/>°(x)i'"(x) • ndX'^'\x) - I j^{x) ■ ndX'^'\x). 



an 



a=l 



(4.3.27) 



dn 



Remark 4.3.8. Note that this does not bring us to the setting of (3.5.5) if there is no 
external volume supply of heat. This is because, in general 

V 

T] p V ^ TjpV. 

a=l 



Compare this to Remark 3.5.4 Using Yla=i ff^p'^v'^ in the entropy inequality fits the approach 
of [25j. This is however rejected by [8l|32j, who go for using T]pv^ as it is mentioned in (3.5.5). 



Remark 4.3.9. The symmetry of the interactions, as imposed by Part|3]of Assumption 4.3.5 



is crucial in the proof of Theorem 4.3.7 However, it is a quite restrictive assumption, that 
disallows many interesting settings and thus deserves further research. At a later stage we 
hope to include a drift in the interactions, which is such that we can still formulate an entropy 
inequality. 
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4.3.3 Generalization to discrete measures 



hi this section, we make feasible that it is possible to derive (at least from a mathematical 
point of view) an analogon for discrete measures of the entropy inequalities mentioned in 



Sections 4.3. 1| and 4.3.2 We already announced that we would do so in Remark 3.5.5 



We consider a single population with corresponding discrete mass measure of the form 



where J7 C N is some index set. The evolution in time of the centres Xi is governed by the 
velocity field v via 

-J-Xi{t) = v{t,Xi{t)). 



dt 



This velocity field is assumed to be of the form (cf. ( 4.3.2| )) 



(4.3.28) 



The convolution VW -k fi is a generalized form of (4.3.3), defined as 



{VW * fi){t,x) := j VW{x-y)dfj,{t,y), for ah t G (0, T), and a; G f]. (4.3.29) 

n\{x} 

The point x itself has been excluded from the domain of integration to avoid interaction of a 
point mass/pedestrian with itself. We assume again that W is symmetric: 

W{C) = W{-0, for ah C e R'^. 



In the spirit of (4.3.6) we define the entropy density rj as 



1 

r]{t,x) := V{x) + -{W -k fi){t,x). 



The corresponding entropy of the system in Q is 



S{t) = J ri{t,x)dfj,{t,x). 
n 

We explicitly restrict ourselves to the situation that all point masses remain in This 
restriction implies that an integral with respect to the measure n can be represented by a 
sum, in which all centres Xi contribute. We have that 



where all positions Xi are time-dependent. The indicator function Ix^en is 1 if Xj G O is true, 
and otherwise. 



''Allowing point masses to leave (or enter) the domain, means that a boundary measure, or a trace of the 
mass measure on the boundary, needs to be defined properly. It is all but trivial to do so for general measures. 



46 Modelling Crowd Dynamics 



Technische Universiteit Eindhoven University of Technology 



The requirement that all centres Xj remain in ft, makes that Ixi&n never becomes 0. For all 
i ^ J thus lx,er2 = 1 holds, for all t G [0,T]. As a result, we lose time-dependence in these 
indicator functions, and we can just write 



Remark 4.3.10. To derive an entropy inequality we examine the time derivative of S. Con- 
sider the 'forbidden' situation that Xi{t) leaves (or enters) the domain, say at time t = t*. 
Then 'i-Xi{t)en is discontinuous in t = t*. We disallow this situation, because the time deriva- 
tive of S does not exist in such point (not even in a weak sense). 



dS 

~dt 



We take the time derivative of S{t): 
dxj 1 



dt 2 



dxj 1 



^ VV{Xi) ■ V{t, + ^ Yj "^^i^i - ^i) • ^(*' - ^ E - ^j) ■ 

Xj^X^ Xj^Xi 

^V{xi) ■ v{t, Xi) + ^ VVF(xi - Xj) ■ v{t, Xi) + l Y ^W{xj - Xi) ■ v{t, xj) 

Y^ViXi)-vit,Xi)+ Y ^WiXi-Xj)-v{t,Xi) 

Yv{t,Xi)-{vV{xi)+ Y ^W{xi-xj)} 
Y^i^^i) ■ {"^Vixi) + (VVF*/x)(x,)} 



\v{t, x)p(i/u(t, x) 



n 

> 0. 



(4.3.30) 



The seventh equality is due to (4.3.29). By (4.3.30) we have derived an equivalent statement 
as in Theorem 14.3.31 



Following similar lines of argument, an entropy inequality as in Theorem 4.3.7 can be derived 
also if we allow subpopulations with both discrete and absolutely continuous mass measures. 
We will, however, not go into further details in this direction. 
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4.4 Derivation of the time-discrete model 

hispired by \n\ H5] , we derive in this section a discrete-in-time counterpart of the weak for- 



mulation presented in (4.1.3) and Definition 4.1.4 



For this aim, we introduce a strictly increasing sequence of discrete points in time {tn}nej\f C 
[0, T]. Here, := {0, 1, . . . , Nt} is the index set, with Nt £ N. The set of points in time is 
chosen such that to = and tiy^, = T. 

We define At„ := tn+i — t n- Since {ininsA^ is a strictly increasing sequence, Atji > holds 
for all n. 

To emphasize that we are working in a time-discrete setting, we write as the time-discrete 
equivalent of •) : B{Q,) — )• M"'' for all n G J\f and all a S {1, 2, ... , z^}. Analogously, we 

write <(•) for ■) : Q ^ R'^. 



For arbitrary n € M, integration in time of (4.1.3) over the interval yields for 

ah ae {l,2,...,u} and ah V'" G C^{^)- 



tn + l 



'0"(x)d;u"(t„+i,x) -y V'"(a;)d;u"(t„,x) = J J v''{t,x) ■V^''{x)dfi''{t,x)dt. (4.4.1) 
n n t„ n 

Assuming all necessary regularity, we can expand the right-hand term in a Taylor series 
around t„. For the sake of brevity we define A{t) := v'^{i, x) ■ V'il^'^{x)dfi'^{i, x)dt. Note 

that A{tn+i) is the right-hand side of (4.4.1). Now 



A{tn+l) = A(t„) + At„ 



d 



dA 
~dt 



t=t, 



Note that A{tn) = and ^^(*) = /n^'"!*'^) ' V?/;"(x)(i/x°(t, x). By writing <(•) instead 



of /i"(t„,-), and w"(-) instead of v"(t„,-) (above we already announced to do so), (4.4.1) 
transforms into 

r(3^)rf/x"+i(x) - j rix)dii'^ix) = Atn J v'^{x)-Vrix)dfi'^^ix) + 0{Atl). 
n n n 

This can also be written as 

J r (x)d/i^+i(x) = J (^"(x) + At„<(x) • VV°(x))d<(x) + 0{Atl). (4.4.2) 
n n 

If v" is 'well-behaved' we expand 

V'"(x + At„t;^^(x)) =V"(x) + AtX(^)-VV°(x) + 0(A4), (4.4.3) 



or 



r (x) + AtMx) ■ (x) = r + At„<(x)) + 0{Atl). 



(4.4.4) 



By 'well-behaved' we mean that is (at least) //"-uniformly bounded. That is, for fixed 
n G M there exists a non-negative constant M„ such that |f"(x)| < M„ for //"-almost every 
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X, by which //"-almost everywhere: At„t>"(2;) = 0{Atn)- 
We substitute this in (4.4.2) to obtain 

r (x)d<+i(x) = I [r{x + At„<(x)) + 0{Atl)y^,-{x) + 0{Atl) 

r{x + At„<(x))d<(x) + 0{Atl), 



(4.4.5) 



where we have taken the ©(At^) -terms outside the integral, and used that /u°(i7) is finite, 
since n"{tn, •) is a finite measure for all choices of n G M. 



We neglect the C'(At^)-part and obtain 

'0"(x)d/i^+i(x) 



(4.4.6) 



Although (4.4.6) is an approximation, we treat it as an equality from now on. Whenever we 
refer to (4.4.6), we thus mean the equality rather than the approximation. In (4.4.6) we have 
moreover used the definition: 

Definition 4.4.1 (One-step motion mapping). The one-step motion mapping Xn defined 
by 

X^(x) :=x + Ai„<(x), (4.4.7) 

for all n £ {0,1, ... , Nt — !}• For simplicity, it will henceforth just be called motion mapping. 
It provides the position at time step n + 1 of the point located in x at time step n. 

Assumption 4.4.2 (Properties of the motion mappings). For all n G {0, 1, . . . , Nt — 1} and 

each a G {1, 2, . . . , z^} we assume that : — ci homeomorphism. This means that: 

(i) Xn ^-^ invertible, 
(a) Xn ^-^ continuous, 
(Hi) {Xn)~^ "^^ continuous. 



Remark 4.4.3. From the proof of Lemma 3.3.1 we extract the statement that a continuous 
mapping is Borel. Regarding Assumption 4.4.2, this implies that Xn {Xn)~^ Borel. 
Respectively, that is: 

(i) for ah n' G ^(!^) we have that {x^y^fl') G ^(!^); 

(h) for ah n' G B{n) we have that Xni^') e S(J1). 

Note that Assumption |4.4.2 is the somewhat weaker counterpart of the assumptions on the 



motion mappings in Section 3.2 



Remark 4.4.4. The expression in (4.4.6) makes sense even if we do not restrict ourselves to 
taking only ■0" G Cg(il). In the sequel we allow ■0" to be any function that is integrable on 
Q with respect to the measure /i^+i, that is: 0° G Lj^a^^{Q,). 
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Fohowing Remark AAA, we can take ■0" = IfJ' the characteristic function for any $7' G 



B{rt). As a resuh, (4.4.6) reduces to 

<+i(J7') = M^^((x:^)-'(f^'))- (4.4.8) 

Definition 4.4.5 (Push forward). The measure ??n+i is called the push forward of the measure 
Tjn via the motion mapping xZ> notation: 

Vn+l = Xn#r]n, (4.4.9) 

if 

is satisfied for all fi' £ B{Q). 



We have now derived the time-discrete version of the problem formulated in Section 4.1 
Definition 4.4.6 (Time-discrete solution). The vector of time- discrete measures: 

is called time-discrete solution, if: 

(i) /i^=o ~ /^o i^ satisfied for each a G {1, 2, . . . , i/}, for some given set of initial measures 
fj,Q that are positive and finite, 

(a) the evolution of fi" is for each a determined by the push forward fJ-^+i = XnH^l^n' 

(Hi) for each a, /i" is positive and finite for all n. 

We refer to finding a time-discrete solution as Problem {P). 

4.5 Solvability of Problem {V) and properties of the solution 

In this section, we prove that there exists a unique time-discrete solution in the sense of 
Definition gAl 

4.5.1 Solvability 

Theorem 4.5.1 (Global existence of time-discrete solutions). Suppose that for each n £ N" 
there exist constants c„ > and Cn > 0, such that for each a E {1, 2, . . . , i/} 

CnX'^i^') ^ X^iiXnY^i^')) ^ CnX'^in'), for each Q' G ]3{n). (4.5.1) 

Suppose furthermore that for each a the initial measure fiQ is given in its refined Lebesgue 



decomposition (cf. Corollary 2.1.10) 



..a _ a I a , a 



where ^aco ^ ^'^> l^do i^ discrete w.r.t. X"^ and fi'^^Q is singular continuous w.r.t. X'^. Assume 
that fJ-acfl! l^do '^"'^ ^sc,o '^'^^ positive and finite measures. 
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Then a time-discrete solution as defined in Definition 4-4-6 exists and it is of the form 

(4.5.2) 



,,o _ ,,0! _|_ a I a 
H-n ~ H-ac,n "r r-d.n "r H-sc.ni 



where fJ-acn ^ f^dn ^'^ discrete w.r.t. , is singular continuous w.r.t. X'^, for all 

n G N , and each component is positive and finite. 



Proof. The proof of Theorem 4.5.1 partly follows the lines of arguments of |17j . 
Let a S {1, 2, . . . , z/} be fixed but arbitrary. The proof goes by induction. The statement in 
(4.5.2) is true for n = due to the given initial condition. Moreover, each of the components 
A*ac,0' ^do ^sc,o is positive and finite. 

The induction hypothesis is that for some (arbitrary, fixed) n E {0, 1, . . . , Nt — 1} the time- 
discrete solution exists, that it is of the form (4.5.2), and that each of the three measures 
A'ac ri) n ^'^'^ l^sc n ™ the Corresponding refined Lebesgue decomposition is positive and finite. 

We now prove that if the induction hypothesis holds for this n, then it also holds for n + 1. 



1. By the definition of the push forward operator formulated in (4.4.8), for any VL' G B{Vt) 
the following holds: 

<+i(i^') = <((x;:)-'(i^')) 

= /^^c,n((x:?)-'(J^')) +/^2,n((xJ^)-'(J^')) +^s%((X^^)~'(J^'))- 

We can thus write 

2. We now define 

..a ,_ ^oiiL a 

We wish to show that /i"^, is absolutely continuous w.r.t. A*^. Let Q,' G B{i^) be such 
that X'^in') = 0. By ( |4.5.l[ ): 

x'iUnrH^')) ^ CnX'in'), 

and thus A'^(O') = implies X'^[{Xn)~^{^')) = 0. It is part of the induction hypothesis 
that < A'', and thus it follows from X'^{{Xn)'^{^')) = that ^ac,n((x")~n^')) = 
0. Thus' 

by which we have proven that /i"^ ^+1 ^ 

Since n positive by the induction hypothesis, we also have that 
and thus A^ac,n+i positive. 

Similarly, finiteness of A'acn+i follows from finiteness of /^ac n- 

where (x")~^(f]) = J7 due to the assumed invertibility of the motion mapping. 
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3. Similarly, define 

Md,n+1 XnifH'd,n- 



By the induction hypothesis //^ ^ is a positive, finite and discrete measure. Lemma 2.1.6 
provides that we can thus write A^d n ~ SieJ' '^i^^i^ some countable index set J'cW, 
a set {xi}ii^j C 17 and a set of corresponding nonnegative coefficients {ajjigj- C M, such 
that Yliej < 

By definition of have that for any $7' G 

= X^^^'^X^C^O^^'^- 
This proves that fi'^^^i is a discrete measure with respect to A'^. 

Positivity of /Uj follows from positivity of the Dirac measure and the fact that ^ 
for all i ^ J . Since Xn maps homeomorphically from to fi, obviously {xi\iizj C 
implies {Xn(^«)}ieJ' C As a result 

i&J iej 

and thus fJ'2n+i ^^^^ finite. 
4. Finally, we define 

..a ._ -.^01 41 ..a 

"sCjTl+l • X.nTTH'sC,'!!' 

d 



We want to prove that this measure is singular continuous w.r.t. A . 
For any x € Q the definition of the push forward implies that n+i(^) — 1^% n{(.Xn)~^ i^)) ■ 
As 

/^sc n is by the induction hypothesis singular continuous, 

/^s%((x°)"'(x)) = 0, for any {x?irHx) G Q. 

Thus A^°c,n+i(^) = foi' a; G 0. 

Let the set -B„ G be such that ^"j. \ = X'^{Bn) = 0, which exists by 

definition of singular continuous measures (see Definition 2.1.7). Because 



the following identity is true: 



(Xn \ Bn 

/^"c,„+i(f^\x:?(i?n)). 
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In the last step we used that Xni^) — ^ {due to invertibility of Xn)- The bottom line 
is that 

The last equality follows from the way we have chosen Bn- 



By hypothesis of the theorem, see (4.5.1), we have that 



A'(x^K)) ^ -X'Ux^y'ix^iiBn))) = -X'iBn) = 0, (4.5.3) 

where it is important that c„ > 0, and X'^{Bn) = by definition of B^- Consequently, 
(4.5.3) proves that A'^(^Xn(^n)) = 0. 



Define Bn+i ■= Xn{Bn)- From the fact that /U^cn+i(^) ~ ^ x G 17, and 

^J■?c,n+l{^ \ Bn+l^ = X'^(^Bn+l^ = 0, 

it follows that fJ-f^n+i ^ singular continuous measure w.r.t. A'^. 
Since fif^. „ is positive by the induction hypothesis, we obtain 

for all W G 13{Q). Thus, fJ-fcn+i positive. Similarly, fJ-scn+i finite because figcn is 
finite: 



5. We have now proven that if a time-discrete solution of the form ( 4.5.2 ) exists for n, it also 



exists for n + 1. Positivity and finiteness of /J^^+i follow from positivity and finiteness of 
its three components, which we have proven above. The induction argument guarantees 



that (4.5.2) holds for all neJ\f. 



□ 

Theorem 4.5.2 (Uniqueness of global time-discrete solutions). Assume the hypotheses of 



Theorem 4.5.1. Then the global time-discrete solution: 



..a _ ..a _|_ a , a 

r-n ~ H'ac,n ^ r-d.n ^ H'sc.n^ 



for all n ^ N 



IS unique. 



Proof. Uniqueness of the time-discrete solution follows from the fact that = Xn#/^n 

defines the push forward unambiguously. To see this, assume that a unique //^ exists, and 
there are two measures and ^2,n+i such that ^ifj^_^i = Xni^l^n each i G {1, 2}. Then 

we have that for any 0,' G B{Q) the following holds: 

Mt„+i(0') = ^°((x^)-i(i7')) = 

thus = ^2,n+l- 

It now follows by induction that if the initial measure /Iq is unique, consequently is defined 
uniquely for all n G M. Note that the decomposition of /x" in its three parts is unique by 
Coronary[2Tl0l □ 
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Remark 4.5.3 (Connections to Reference [17]). The resuhs of [T7j are a special case of our 



setting; we already said so in our comments at the end of Section 4.1 



Let mo be a discrete measure: mo := "^jLi^Pj, and let Mq be an absolutely continuous 
measure (with density p{t, x)). There are two ways to recover their situation. One way is to 
set v = 1 and consider 

/xj = 6'mo + (1 - e)Mo, 
where € [0, 1] is a tuning parameter. 

The second way to recover the results of [T7] is by setting v = 2, and then by taking 



Brno, 

{l-e)Mo, 

henceforth only considering the total mass measure and the barycentric velocity. 



4.5.2 Basic properties of the time-discrete solution 

In this section we will formulate and prove a number of properties of the time-discrete solution 
provided by Theorems |4.5.1 and 4.5.2 These properties follow mainly from the subsequent 
steps done in the proofs of those theorems. 



Corollary 4.5.4 (Conservation of mass). Assume the hypotheses of Theorem 4-5.1 For each 
n £ M the initial mass is conserved by each of the three components of the time-discrete solu- 



tion provided by Theorems 4-5.1 and 4-5.2, That is, for all n £ M and each a G {1, 2, . . . ,1^}: 
(z) <,,„(17) = M°,,o(^^)' 

(iU) /i-,„(l^)=/i%(^^)- 



As a result, also fJ-ni^) = /^o(^)- 

Proof. For each n G {0, 1, . . . , Nt — 1} and a G {1, 2, ... , v}, consider the measure p^ 



4.5.1 



where 
Parts 



00 G {ac, d, sc}. By definition of /x^ „+i (see the constructive proof of Theorem 
§[3] and [4]): 

..a . ^au a 

H-ui,n+l ■~ XnTTH'ojjn- 

For each n we thus have 

and {xn)~^{^) = ^ due to the invertibility of Xn- This implies that fJ-Z n+ii^) = /"3,n(^) for 
all n G {0, 1, . . . , Nt}, and by an inductive argument: /i^ n(^) = o(^) each n. 
Since = n + /^d n ^ /^sc n ^11 n G N, the above implies trivially that 



/.^(0) = M^(1^), 



for all n G J\f. 



Corollary 4.5.5. Assume the hypotheses of Theorem 4-5-1 If ^jl'^cO = 0' /^do = ^' 



□ 



or 



/^sc = respectively, then for each n the corresponding component of /i" vanishes. That is, 
for each a, the following statements are true: 



(i) fJ-acO = implies Hac,n = for all n G M , 
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(ii) /i^o = implies = for all n G J\f, 
(Hi) fj."^^ = implies fi^cn = /^^^ ^ 

Proof. For any a G {1, 2, . . . , z^}, let w G {ac, d, sc} be arbitrary. Then the corresponding 
part in the refined Lebesgue decomposition of n^+i follows uniquely from ^ by the 

push forward n^^n+i — Xn#^3,n- See Parts [2| |3] and |4] of the proof of Theorem 



4.5.1 



Assume that /i^ n = 0- Then for each W G B{Q) 
We thus conclude that /x^ n+i = 0- 

It follows by an inductive argument that /^S n = ^ if ^3 o = ^ given. □ 

Remark 4.5.6. If /i"^ o = 0' then the assumption that there is a Cn > such that 



can be removed from the theorem. This is because, even without this assumption, the push 
forward of /x^c „ = is an absolutely continuous measure (of course, more specifically fj,"^ = 
0). 

Due to similar arguments, fi'^^ o = ^ allows us to remove the assumption that there exists a 
c„ > for which 



Remark 4.5.7. If we set fJ-^co = ^^'^ /"sc o = then Theorem 4.5.1 provides that /i" is 
discrete w.r.t. A'^ for all n G M. We have shown in Lemma 2.1.6 that any discrete measure 
w.r.t. A'^ can be written as a linear combination of Dirac measures. Note however, that in 
both the centres of the Dirac masses and the coefficients might in principle depend on 
n. (Since the index set J7 is necessarily countable, is suffices to allow only the positions and 
coefficients to be n-dependent.) 

However, due to the specific structure of the push forward, we have the following corollary 
of Theorem 14.5.11 

Corollary 4.5.8. If the initial measure is a linear combination of Dirac measures: ;Uq = 
Yliej diSxi for some countable index set J , then for each n the measure is of the form 
= ^i^j o^i^x^^, with the same (constant) coefficients Oj. 

Proof We extract from Partjsjof the proof of Theorem 4.5.1 that for any n G {0, 1, . . . , Nt — 
1}: < = T^iej^i^xi implies //"+i(0') = ^^^j ai6^..(^^){n'). Assume that = EieJ«»'^^xo 
has been given. By inductive reasoning we can now deduce that = '^^iej o-i^x^ „ for all 
n £ M, where 

Xi^n ■■= {Xn-l ° Xn-2 ° • • • ° xt){Xi^<d)- 

□ 

Lemma 4.5.9. // = E fflif^Xi „ for each n £ M then, for any i £ J , the position 
Xi^n satisfies a discretized version of dxi[t)/dt = v'^{t,Xi{t)), evaluated at t = tn for all 
nG{0,l,...,A^T-l}. 
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Proof. For any i E ^7, a Taylor series expansion of Xi around t = tn yields (evaluated in 

t = tn+l)'- 

Xi{tn+l) = Xi{tn) + Mnj^Xiitn) + C(At^), 



by which we find the expression 



dt 



Xi{tn+l) - Xi{tn) 



(4.5.4) 



The discretized version of dxi{t)/dt = v"'{t,Xi{t)) is now obtained by substitution of (4.5.4), 
and after neglecting the C'(At^)-part. Write Xi^n as the discrete-in-time equivalent of Xi{tn), 
and furthermore use v!^{xi^n) = V^it 

njXi,n) like before. The discretized version becomes 



or 



(4.5.5) 



We have seen in the proofs of Theorem 4.5.1 and Corollary 4.5.8 that Xi^n+i and Xi^n are 
related via 



and by definition of the motion mapping (see Definition |4.4.1 ) this is exactly (4.5.5). 



□ 



Remark 4.5.10. In Remark 3.4.3 we have already anticipated the statements of Corollary 
14.5.81 and Lemma I4.5.9[ 

Lemma 4.5.11. Assume that, for each n £ M and a G {1, 2, . . . , ly}, the velocity field is 
Lipschitz continuous, with Lipschitz constant strictly smaller than 1/At„. That is, there is a 
constant ^ Kn < 1/Atn, such that 

- v"{y)\ ^ Kn\x - y\, for all x,y en. 

For each n S A^, the measure jj!^ „ is discrete, and thus it is the linear combination of Dirac 



masses (Lemma 2.1.6). Let {xi^n]i^j be the corresponding set of centres (at time-slice n). 



Now, if {xifi}i(zj consists of distinct elements, then also {xi^n}iej consists of distinct ele- 
ments, for each n £ M. 

Proof. The proof goes by contradiction. 

Assume that n £ {1, 2, . . . , N^} is such that not all elements of {xi^n}iej are distinct, and, 
more specifically, let j,k £ J satisfy j ^ k and Xj^n = 3;fc,n- 

We know that = {Xn-i ° Xn~2 ° ■ • ■ ° Xo)i^i,o) foi' all i £ J', as this was shown in the 
proof of Corollary 4.5.8[ 

By assumption of the lemma Xj^ ^ Xkfi- For xj^n = Xk,n to hold, there must therefore be an 
h £ {0, 1, . . . , n — 1} such that 

Xj,h 7^ Xk^fi, 

but 

By definition of the push forward, the last equivalence can be written as 

Xj,h ~\~ ^l^nVfiixj^n) = Xk^fi + Atjjf^ (xfc^n) ; 
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or 

Atn{v^{Xj^n) - V^{xk^fi)) = Xk,h " Xj^fi- (4.5.6) 

By assumption of the lemma, v"^ is Lipschitz continuous, and there is a Kn < 1/Atn such 
that 

- v^{y)\ ^ Kn\x-y\. 

In particular, this implies that 



However, as a result of (4.5.6) we already have that 



Atn 



and thus we have a contradiction. This finishes the proof. □ 

4.5.3 Relaxing the conditions on the motion mapping and velocity fields 

Up to now we have made a number of assumptions with respect to the one-step motion map- 
ping the (discretized) velocity field v". We first give a short recapitulation of these 
assumptions here. Afterwards we indicate which of these assumptions can be relaxed, and 
in what way this can be done. The results of this section have a preliminary character: we 
expect that more is to be discovered as soon as more analysis effort is invested in this direction. 

The following restrictions on the motion mapping and velocity field were imposed so far: 



1. In Section 4.4, we assumed to be //^-uniformly bounded. We needed this to be able 



to write the Taylor expansion in (4.4.3). This assumption means that, for each n G M 



and a £ {1,2, ... ,1^}, there exists a constant Mn for which 

< Mn, for /i"-almost every x G Q. 



2. By Assumption 4.4.2, we postulated that Xn • ^ ~^ ^ is a homeomorphism, for all 
n £ {0,1, ... ,Nt - 1} and a £ {1,2, ... , u}. 

3. In case /i"^. q ^ and /i^^, q ^ 0, we need for all n G {0, 1, . . . , Nt — 1} strictly positive 
constants c„ and C„,, such that 



for each Q' £ B{n), 



to prove existence and uniqueness of the time-discrete solution (see Theorems 4.5.1 and 



4.5.2). 



4. The Lipschitz continuity of the velocity field is demanded in the hypothesis of Lemma 
4.5. 11[ More specifically, we assume that for each n G {0, 1, . . . , Nt — 1} a constant 



Kn < 1/Atn exists, such that 

IKix) - v^{y)\ ^ Kn\x-y\, 



for all x,y £ Q. 



This restriction is only needed fi^o ^ 0> to make sure that no two centres of Dirac 
masses can 'merge'. 
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Before we can relax these conditions, we introduce the following concept: 

Definition 4.5.12 (Support of a measure). The support of a (positive, finite) measure n on 
B(i}) is defined as the set 

supp fi := {x €z i} \ > for all ^l' G B{^}) open such that x € 0,'}. 

One can show (see [l3], pp. 27-28) that the support is a closed set of full measure, that is 

/i(supp fl) = fi{fl). 

As a result 

^(il \ supp fj.) = /u(r2) — fi{il n supp /i) = /i(l^) — /x(supp fj,) = 0. (4.5.7) 

For any measurable set Q' C 0,\ supp /x it follows that 

H{n') lj.{n\suppfi) = 0. (4.5.8) 

We have introduced the concept of a support, because it is in fact only important to consider 
the properties of the motion mapping Xn supp /i". As (4.5.8) shows, there is no mass con- 
tained in any region outside the support. For our purposes, it is irrelevant whether the image 
of such Q' C 0,\ supp /i" is again a subset of 0, or not. In order to have mass conservation 
within the only thing that matters is whether Xni^') ^ ^' supp^". 

We do not change the restriction in Part [T] of the list above. However, we remark that if 
we demand that |7;"(2;)| < M„ holds for all x £ supp^", then it also holds for ^^-a.e. X. 
Indeed, if |w"(2;)| < M„ for all x E supp/U°, then 

■.= {xen\Mni^ Kix)\} C n \ SUpp/i^, 



thus, cf. ( |4.5.7D , 

^J^(17M)^<(f^\suppO = 0. 



In Part [2] of the list of assumptions, we pay special attention to the restriction of Xn 

SUpp<, that is, Xn ■= XnLupp/.- 

mapping is a homeomorphism by: 



supp/i", that is, Xn '■— ^n|supp/i"' Moreover, we replace the assumption that the motion 



Assumption 4.5.13. For all n G {0, 1, . . . , Nt — 1} and each a G {1, 2, . . . , i/} we assume 
that the motion mapping : — I^'^ is such that the following properties are satisfied: 

(i) the range of xZ ■= Xn\snpp,i'^ contained in Q. That is: Xn ■ supp^u^ ^; 
(a) for all G B{^) we have that {Xn)~^{^') is measurable; 
(Hi) for all Q' G B{Q) we have that x"(i7') is measurable. 

Here, {Xn)~^ should not be understood as the inverse mapping. The set (x")"^^') is the 
pre-image of Vt' . We do not assume invertibility of the motion mapping. 
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Note that throughout this thesis, anywhere we used the inverse mapping (Xn) '^^ should 
now read it in the pre-image sense. This is e.g. also the case in the definition of the push 



forward, Definition 4.4.5 



Under the new set of restrictions on the motion mapping, we still want the Theorems and 



Lemmas of Section 4.5 to hold. We want this especially for Theorems 4.5.1|and|4.5.2 (global 



existence and uniqueness of the time-discrete solution). To achieve this. Part [S] of the list of 
assumptions needs reconsideration. We make it slightly stronger: 

Assumption 4.5.14. Suppose that for each n £ {0,1, Nt — 1} there exist constants 
On > and Cn > 0, such that for each a G {1, 2, . . . , i/} 



hold for all Q' G B{Q). Again, {xn)~^{^') is the pre-image of O! . 



(4.5.9) 
(4.5.10) 



Note that (4.5.9) implies the left inequality of 

CnX\^') ^ A^((xJ^)-^(^^')) ^ CnX\n'), 

by substituting {xZr\^') for ^' , and by using xZ[{x'^r\^')) = ^' ■ 

We still demand that the velocity field is a Lipschitz continuous function if /^^q ^ (Re- 
striction [4]). We relax this restriction in this sense, that it is no longer required for all 
X, y G i7. It suffices to demand this property only for all x, y G supp /^^ Note that supp fjf^ ^ 
consists exactly of those points in which the Dirac measures of /i^ ^ are centered. In the proof 



of Lemma 4.5.11 t;" is only evaluated in these points. The proof is therefore still valid under 
the proposed weaker assumption. 

The results of Section |4.5| are still valid under the new assumptions proposed in this Section. 



We only need to modify some minor details in the proofs of Theorem 4.5.1 and Corollary 



4.5.4 These modifications are addressed in Appendix [E} 

Remark 4.5.15. In Lemma 14.5.111 we have shown that the centres of the Dirac measures 
cannot merge if the velocity field is a Lipschitz continuous function. We acknowledge that it is 
equally important that two Dirac masses that belong to distinct subpopulations cannot merge. 



Moreover, the integrand in the social velocity terms, see (4.2.2), typically has a singularity 
in 0. This means that at timeslice n, the velocity field v°(x) is unbounded if x is arbitrarily 
close to a Dirac mass. We thus propose as an extra demand that for all a G {1, 2, . . . , z^} and 
for each n G N 



supp 11^ n supp ^5 , 



for all (5 G {1, 2, . . . , I/} such that f5 ^ a, 



and 



supp/x^,„ n (supp/i^^c_„ n supp/i^c_„) = 0. 

If this condition is obeyed, no Dirac measures from distinct populations can be centred in 
the same point. Since the support of a measure is a closed set, we furthermore prevent f"(x) 
from being unbounded for all x G supp//^. What happens outside supp/u" is unimportant. 
We come back to this issue in Section [H 
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4.6 Reformulation of the proposed velocity fields in the time- 
discrete setting 



hi Sections 4.4 and 4.5 we have used the velocity fields without specifying their actual 



form. In this section we want to do so. It is quite self-evident that we just take the proposed 
form of Section [4. 2| and determine its time-discrete counterpart. 

For a G {1, 2, . . . , i/}, we thus take as the superposition of the desired and the social 



velocity like in (4.2.1): 



<(x) := v^esix) + <c,n(a;), 
For Vg^^ „ we have the following expression: 



for all X £ Q. 



^n\{x} 



y-x 
y- x\ 



for all a E {1, 2, . . . , I/}. 



4.7 Entropy inequality for the time-discrete problem 

In Section [4. 4| we have derived a time-discrete version of our problem. We derived an entropy 



inequality for the continuous-in-time problem in Section 4.3 Its discrete-in-time counterpart 



is presented in this section. We again work with absolutely continuous mass measures only, 
and treat (for clarity) the single-population case first. 



4.7.1 One population 



In Section 4.4 we derived, omitting C'(Ai^)-terms, (4.4.6): an equation relating the situa- 
tion at time t„ to the one at time tn+i- If the time-discrete mass measures are absolutely 



continuous, (4.4.6) reads 



'4>{x)pn+l{x)d\'^{x) = J ij{xn{x))pn{x)dX'^{x), 

n n 



(4.7.1) 



where we can take any ^p S ^li n+i (^ ) (^^- Remark 4.4.4). We deduce our discrete-in-time 
entropy inequality here, taking ( 4.7.1| ) as a starting point. 



The subscript n in the sequel denotes that we consider the time-discrete version of a function 



at time t„. Recall that Xn denotes the one-step push forward (see Definition 4.4.1), defined 

by 

Xnix) := X + AtnVn{x), 

Let Sn be the time-discrete equivalent of the entropy S{tn). The following theorem holds: 

Theorem 4.7.1 (Discrete-in-time entropy inequality). Assume that the time- discrete velocity 
Vn is of the form 

Vn:=VV + VW*pn, (4.7.2) 
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where W is such that W{^) = W{—^) for all ^ G M'^. Assume moreover that the evolution of 
the system is governed by ^J^TA^), and that the entropy is defined by 



Sn--= hvix) + -{W * pn)ix))pnix)dX'^{x), for all Tl & M . 



(4.7.3) 



Then for each n G {0, 1, . . . , Nt — 1} the following inequality holds, up to 0{At'^) -terms: 



Proof Fix n € {0,1, ... ,Nt - 1} and take 



By (4.7.1), we derive 

Sn+l 



(Vix) + -{W * Pn+l)ix))pn+lix)dX\x) 
{V{Xn{x)) + \{W * Pn+l){Xn{x)))pn{x)dX''ix). 



W{x -y- AtnVn{y))Pn{y)dX'^{y) 



W{x-y)pn{y)dX\y) 



n 



(4.7.4) 



(4.7.5) 



Note that, cf. - (|4A4]): 

V{Xn{x)) =V{x) + AtnVn{x) ■ VV {x) + 0{Atl), (4.7.6) 

{W ★ Pn+l){Xn{x)) ={W * Pn+l){x) + At„W„(x) • V {W ★ Pn+l){x) + ©(Ai^) 

= {W * Pn+l){x) + AtnVn{x) ■ {VW * Pn+l){x) + 0{Atl), (4.7.7) 

{W^Pn+i){x) = [ W{x-y)pn+i{y)dX''{y) 



W{x-Xn{y))pn{y)dX\y) 



- Atn J vn{y) ■ ^W{x - y)pn{y)dX\y) + 0{Atl) 
n 

--{W*Pn){x) - At„ J vn{y) ■ VW{x - y)pn{y)dX''{y) + 0{Atl). 



(4.7.8) 



In the second equahty of (4.7.8) we have used (4.7.1 ) with ^/'(y) = W{x — y) (in this scope, x 
is regarded as a parameter). Note that (4.7.8) can also be written as 

{W^Pn+l){x) = {W^Pn){x) + 0{Atn), (4.7.9) 
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which we use now. Combination of (4.7.7) - (4.7.9) yields 

{W * Pn+l){Xn{x)) ={W * Pn+l){x) + Af„t.„(x) • V {W ★ p„+l)(x) + 0{Atl) 

=iW ^ Pn){x) - Atn I Vn{y)-VW{x-y)pniy)dX'^{y) 



--{W ^ Pn){x) - Atn [ Vniy)-VW{x-y)pniy)dX\y) 



+ AtnVnix) • J VW{x - y)pn{y)d\''{y) + 0{Atl). (4.7.10) 



Substituting ( |4.7.6D and ( |4.7.10D into ( |4.7.5D , we find 

5„+i = / {V{Xn{x)) + ^(VF*p„+i)(Xn(:r)))/,,(x)dA'^(x) 

+ AtnVnix) ■ VVix))pn{x)dX\x) 
{W*Pn){x)pn{x)dX''{x) 

vn{y) ■ yW{x - y)pn{y)dX\y)pn{x)dX\x) 



n n 



+ ^Atn I vn{x) ■ I VW{x - y)pn{y)dX''{y)pn{x)dX\x) + 0{Atl). (4.7.11) 

Since W{i) = W{-i) for ah ^ e M^, we have that VW{i) = -VW{-i). Using some 
elementary calculus, we derive 



1 



vr,{y) ■ VWix - y)pn{y)dX%y)pn{x)dX\x) 



1 



At„ / vn{y)- / VW{x-y)pn{x)dX''{x)pn{y)dX\y) 



Atn / vn{y)- / VW{y - x)pn{x)dX'^{x)pn{y)dX'^iy) 



=^Atn J vn{x) ■ J VW{x - y)pn{y)dX\y)pn{x)dX\x). (4.7.12) 

In the last step we changed notation, replacing x by y and vice versa. In this final expression 
we recognize the last term of (4.7.11), which can in short be written as 

\Atn I Vn{x) ■ {VW * Pn){x)pn{x)dX'^{x). 
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The expression for Sn+i in (4.7.11) now transforms into 



Sn+l = J {V{X) + AtnVnix) ■VV{x))pn{x)dX'^{x) 

n 



1 



{W*Pn){x)pn{x)dX''{x) 



n 



+ Atn J Vn{x) ■ J VW{x - y)pn{y)d\\y)pn{x)d\\x) + 0{Atl) 
n n 

dr 



{V{X) + AtnVn{x) ■ VV{x))pn{x)dX''{x) 



n 



+ o iW^Pn){x)pn{x)dX''{x) 



Atn j Vn{x) ■ {VW * pn){x)pn{x)d\'^{x) + 0{Atl) 



{V{X) + -{W^pn){x))pn{x)d\\x) 



+ Atn j Vn{x) ■ {VV{X) + {VW i. pn){x))pn{x)d\\x) + 0{Atl). (4.7.13) 



Now we recognize the definitions of the entropy Sn, (4.7.3), and the velocity Vn, (4.7.2), in 



the last equality of (4.7.13). We thus proceed: 



Sn+l =Sn + Atn j \Vn{x)Y Pn{x)dX\x) + O(At^) 
+2\ 



^Sn + 0{Ati 

Neglecting C(Ai^) -terms, we thus have that 



(4.7.14) 



Sn+i-^Sn, forallnG{0,l,...,iVT-l}. 



□ 



Remark 4.7.2. Note that (4.7.4) is the discrete-in-time counterpart of (4.3.13), up to 
C'(At^)-approximation. 

Remark 4.7.3. Unlike the continuous-in-time case, no boundary terms appear in this time- 
discrete context. In Section [4.3.1[ indeed these terms did appear, and we either assumed them 
to vanish (isolated system), or we incorporated them in the entropy inequality. However, in 



the proof of Theorem 4.7.1| no boundary terms are encountered. This is because, originally, 
the test functions ij) were taken such that they vanish on dO,. This implicitly assumes that 
the time-discrete model is more suitable for describing an isolated system. 
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4.7.2 Multi-component crowd 

For a crowd that consists of multiple subpopulations, for each constituent of the mixture a 



governing equation like (4.7.1) holds: 



(4.7.15) 



for any S L]a (17). The one-step push forward Xn) is defined by 



Xnix) := x + AtnV^ix), 



where we take the time-discrete velocities of the form 



■= + * for ah a G {1, 2, ... , u}. 

13=1 



(4.7.16) 



We assume that for ah a,^ E {1, 2, . . . , z/}: W^{i) = W^{-$) holds for ah ^ e W^, and 

= wL 

Very much as in Section 4.7.1 we define the time-discrete entropy for the multi-component 
crowd to be 

^-■■=Y1 j[v''ix) + -Y,iW^^P?,)ix)]p^ix)d\''ix), forahneAA. (4.7.17) 



Theorem 4.7.4 (Discrete-in-time entropy inequality). Assume that satisfies {4-7. 16) and 

the accompanying conditions on W^. 



Assume furthernnore that the evolution of the system is governed by {4-7. 15), and that the 
entropy is defined by \4-7-17 ). 

Then for each n G {0, 1, . . . , Nt — 1} the following inequality holds, up to 0{/S.t^) -terms: 



(4.7.18) 



Proof. Fix n G {0, 1, . . . , iVy - 1} and take V'" := V + - Yl%i{Wp ★ Pn) for each a G 



{1,2,..., v}. Using (4.7.15), we derive that 

S„+i = j2j [V%x) + \ Y^iW^^ ^ pi^^){x)\pl^^{x)d\\x) 



a=i; 



/3=1 



Y [ [v^"(x^(^)) + 2E(H^;?*/'f+i)(x:?(x))]p^^(x)dA'^(x). (4.7.19) 



For fixed a, /3 G {1, 2, . . . , i^} similar expressions as in (4.7.6), (4.7.7), (4.7.8) hold. They can 
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be derived in an identical manner. Therefore, we list here only the results: 



(4.7.20) 



=(^^; + At„<(x) • {VW^ * pf+i)(:r) + 0{Atl), (4.7.21) 

(H^^" * pf+i)(x) =iW^ * y9^)(x) - At„ J vi{y) ■ VW^{x - y)pi{y)d\\y) + 0{Atl). 

n 

(4.7.22) 



We combine (4.7.21) and (4.7.22), and (omitting the details) obtain the analogon of (4.7.10): 



(Vf^;? *pf+i)(x^(^)) =iW^*p?,)ix) - AtnJ v^iy) ■ VW^ix - y)pi{y)d\''{y) 

+ At„<(x) • j VW^ix - y)pi{y)dX\y) + 0{Atl). (4.7.23) 



By substitution of (4.7.20) and (4.7.23) into (4.7.19) we get 



^-+1 = E / + I * P^+iKx'^ix))] P^^(x)dA^(x) 

= [ (v^"(2;) + At„<(x) • vy"(x))/)°(x)dA^(x) 

+ / Ew*p^)(^K(^)dA'^(x) 

a=l^ 15=1 

-^^*"E / E / ^^?(2/)-VW^;(x-y)/.^(2/)dA'^(y)pJ^(x)dA^(a;) 

+ ^^*-E / E<(^)- / yW^ix-y)p^niy)dX''{y)p'^ix)dX'ix) + OiAtl). 

a=l^ 13=1 I 



(4.7.24) 



We now use that W^{i) = W^(-0 for aU ^ G M"^. This implies that VW^(0 = -VW^{-i). 
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We obtain 



-lAtnj^ fjzl v^{y)-VW;^{x-y)pi{y)d\%)p':,{x)dX\x) 

■lAtnj^fl I ^^niy)- [ VW^iy-x)p'^ix)dX'^ix)p^MdX'iy) 

■lAtnj^fl J / '^Wiix-y)p^^iy)dX\y)p:{x)dX''ix) 

/5=i"=in n 

■l^injZ I jZ<^'^) ■ {'^W^^pi){x)p':,{x)dX\x). (4.7.25) 



In the third step we changed notation: we replaced x by y, a by /?, and vice versa. 



The expression for Sn+i in (4.7.24) transforms into 



Sn+i =j2[ (^"(^) + • Vy"(x))p^(x)dA'^(x) 

"=1q 

+ / j2iW^*p?,){x)p?^{x)dX%) 

+ At„X^ / f]<(x).(^(VH^; + VW^,^)*pf,)(xK(x)dA'^(x) + 0(A4) 
= E / [l^"(x) + 2E(W^;*p^)(x)]p^(x)dA'^(x) 

+ Atnj2 [ t'nl^^)- [vy°(x) + ^X^((VTy^" + VT^„^)*p^)(x)]p^(x)(iA'^(x) + 0(At2 

=5, + At„ / K(x)|Vn(^)dA^(^) + O(At^), (4.7.26) 

-=in 

where the last step is only valid under the assumption that the interactions are symmetric 
(this is a hypothesis of the theorem) . 



If we omit the C'(At^)-terms, we deduce from (4.7.26) that 

Sn+i ^ Sn, for all n £ {0,1, ... ,Nt - 1}, 
which finishes the proof. 



□ 



Remark 4.7.5. Note that (4.7.18) is the discrete-in-time counterpart of (4.3.26), again up 
to C'(At^)-approximation. 
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4.7.3 Generalization to discrete measures 



Following the lines of Section 4.3.3, we generalize the entropy inequality for the time-discrete 



model (cf. Theorem 4.7.1) to the situation in which the discrete-in-time mass measure is of 
a discrete type. 

We consider a single population with discrete-in-time mass measure 



The push forward of the centres is given by (4.4.7) such that 



Xi^n+l '■— Xn{Xi^n) — Xi^n + ^tnVn{Xi^n) i for all 2 G ^7^. 



(4.7.27) 



Here, the velocity field is the time-discrete version of (4.3.28): 



Vn := Vy + VW*/in, 

where the convolution is given by 

{VW-k^in){x) := j VW{x - y)dfiniy), for ah x G n. 

n\{x} 

We assume again that W is symmetric: 

l^(^) = W{-0, for ah C G M^. 

The corresponding entropy of the system in Q is 

5'n = y \v{x) + ^{W * iJ,n)ix)^d^inix), for ah n G 7^. 
n 

We again explicitly restrict ourselves to the situation that Xi^n £ O for all i G J7 and all 
n G AA; cf. Section |4.3.3 



Our aim is to derive the entropy inequality of Section 4.7.1 for the proposed microscopic 
measure. That is, up to 0(At^)-terms, we have that 

Sn+i>Sn, for allnG {0,l,...,iVT- 1}. 

Using the definition of /i„ (in terms of the sum of Dirac measures), the relation between Xi^n-\-i 



and Xin in (^4.7.27), we derive that 



Sn+1 = ^{V{xi^n+l) + ^ Yl W{Xi,n+l-Xj,n+l)] 

J , Ti + 1 7^ ^ i , n + 1 

|^(2;j,„ + At„i;„(xj,„)) -h ^ ^ W{xi^n + AtnVn{Xi.a) - Xj^n " AtnVn{x 



J , n -|- 1 7^ ^ i , n + 1 



(4.7. 
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Remark 4.T.6. Note that the statement Xj^n+i 7^ Xi^n+i is equivalent to Xj,„ ^ Xi^n, if Vn is 
a Lipschitz continuous function with the Lipschitz constant strictly smaller than 1/At„. The 
implication 

follows trivially from the unambiguity of the push forward, whereas the implication 



is a consequence of similar arguments as in the proof of Lemma 4.5.11 



We can thus proceed on (|4.7.28|): 

Sn+l 



nyXj^n, 



+ 1 E {W^(x.n-X,,„)+At„(.„(x.„, 



))-VW{Xi,n-Xj,n) + 0{Atl)y 



(4.7.29) 



Here we have used Taylor series expansions around (in V) and around Xi^n — xj^n (in W) , 
respectively. Note also that VW{xi^n — Xj^n) = —^W{xj^n — Xi^n)- We can thus write: 



Sn+l = 'Y{y{Xi,n) + ^ E W{Xi^n - Xj^n)^ 



+Atn E Vn{Xi^„) ■ VV{xi^ri) 



+ ^Atn E Vn{Xi^n) -^Wixi^n - Xj^n) 

+ ^Atn Yl Mxj,n) ■ ^W{Xj,n - Xi,n) + 0{Atl 



Sn + Atn E Vn{xi,n) ' VV{xi^n) 



) • VW{xi,n - Xj,n) + OiAti 



Sn + Atn Y 

Vn {Xi^n )-^W{xi^n)+ Y ^W{Xi^n-Xj^n)^ + 0{Atl) 

(4.7.30) 
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Omitting 0{Atl)-teims, ( 4.7.30| ) we arrive at 

= Sn + Atn j |w„(x)pd/i„(x) 



(4.7.31) 



The inequaUty in (4.7.31) is an analogon of the statement of Theorem 4.7.1 for a discrete 
mass measure. We do not give details on the multi-component case. 



4.8 Two-scale phenomena 

At this stage, we want to point out which specific types of measures we have in mind to 
use in our model. Throughout the previous sections we have already given clues about the 



choices we intend to make, for example in Sections 1.3 2.4 3.4 and |4.1[ As anticipated, we 
are mainly interested in two sorts of measures. 

The first option is a discrete measure (or: microscopic mass measure), as introduced already 



in Section 2.4.1 In general, a discrete measure may be the (weighted) sum of a countable yet 



infinite number of Dirac measures (cf. Lemma 2.1.6). We restrict ourselves here to a finite 



but arbitrary - number of point masses (read: people). This approach enables us to trace the 
exact motion of a particular individual. 

Secondly, if we work with an absolutely continuous mass measure (w.r.t. the Lebesgue mea- 
sure), then we are provided with a density by the Radon-Nikodym Theorem. This option 



(a macroscopic mass measure) has already been introduced in Section 2.4.2 In fact, this 
perspective corresponds to considering the crowd as a fiuid (or a large 'cloud'). We decide to 
do so if the number of people is large, and if we, moreover, are not interested in what happens 
exactly at the individual's level. Such local fluctuations have been averaged out. 

We only allow these two options, which means that from now on we explicitly exclude singular 
continuous measures from our area of interest. This is because it is not self-evident what in- 



terpretation we should give to such measure. This decision has been anticipated in Section 3.4 



In Section [T| we have already indicated that we are especially interested in the interplay 
between microscopic and macroscopic mass measures. Piccoli et al. [17] do so by including 
a discrete and an absolutely continuous part in one measure. Furthermore a tuning parame- 
ter (taking values from to 1) is used to enable a transition from fully microscopic to fully 
macroscopic. Speaking in terms of mixture theory (cf. Section [s]), this means that the dis- 
crete and absolutely continuous part together describe a single constituent. They can not 



act independently. We already remarked in Section 4.1 that we consider this to be not very 
useful. In the numerical scheme that is described in Section [5j we therefore do not include 
this superposition of a micro and a macro part in one measure. If desired, it can however be 
incorporated without too much difficulty. 
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Our intention is a bit different: We namely want to distinguish between subpopulations and 
assign to them a measure that corresponds to their 'character'. More specifically, we want to 
use a discrete measure only if this subpopulation consists of a limited number of individuals 
which are of special interest. This is why we only consider discrete measures that are consti- 
tuted of a finite number of Dirac measures. Large crowds are in our framework automatically 
represented by an absolutely continuous mass measure. Such subpopulations consist of people 
that are not special (otherwise they should have been included in an independent discrete 
measure), and therefore it suffices to be able to observe only global behaviour. It is of no use 
to 'do effort' to capture information from the microscopic level in a discrete measure, if we 
are not interested in that information any way. The use of both microscopic and macroscopic 
mass measures in one unified framework, makes that we work in a two-scale setting. Our 
perspective is twofold: on one hand we observe from a macroscopic point of view, while on 
the other hand we can detect behaviour at an individual's level if we have identified this 
individual as 'special'. 

Remark 4.8.1. In the future, we possibly want to come back to modelling the cloud as a 
particle system, consisting of a large number of individuals. We then aim at comparing the 
particle system's results to the results of this thesis. Note that these two approaches are 
incorporated in the measure-theoretical framework we are dealing with. 

By introducing separate mass measures, we allow the subpopulations to evolve differ- 
ently. For example, we assign to each subpopulation its own desired velocity field. That 
is, for each subpopulation there is a separate goal that it wants to achieve. Moreover, it is 
possible to have asymmetric interactions. We have already shortly referred to what we called 



"predator-prey relations" in Remark 4.3.6 In such asymmetric situation preys are only re- 



pelled from predators, but predators are attracted to preys. 



Nature turns out to provide very clear illustrations of systems that are of a two-scale character 
and at the same time display asymmetric interactions of the "predator-prey" type. We point 
out two examples. Fascinating interaction takes place between fiocks of small birds (starlings, 
say) and one or more larger predator birds (e.g. hawks), as is illustrated by Figure 4.4 The 
huge number of starlings makes it nearly impossible to distinguish between individuals in the 
fiock. One clearly perceives a continuum-like cloud of birds; modelling this cloud by means 
of an absolutely continuous measure thus seems natural. As the hawk attacks the starlings, 
the flock reacts to the approaching enemy as if there were some macroscopic coordination. 
Increase and decrease of the density in the flock can be observed. It is striking to see that the 
group does not fall apartj^ Similar effects can be seen in shoals of flsh being attacked. The 
interaction between the predator flsh and the shoal triggers complex structures to appear and 
enables sudden transitions between macroscopic patterns. The phenomena we refer to are 
described e.g. in |36l [5T]. 



Two-scale effects also occur in human crowds, when there is special interaction between a 
speciflc individual (or a limited number of them) and the rest of the crowd. These phenom- 
ena might be harder to visualize (because the individual and the group members have the 
same size), but this does not mean that they are not there. Phenomena of the "predator- 
prey" type occur, for example when a crowd is attacked by some criminal or terrorist. More 

■'Movies of such phenomena are available on the internet, see e.g. 
www . youtube . com/watch?v=b8eZJnbDHIg&f eature=related . 
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'friendly' interplay is present when the special individuals have the role of tourist guides, 
leaders, firemen, safety guards et cetera. Our aim in Section [5] and further is to capture these 
two-scale phenomena with our model. 




Figure 4.4: Illustration of two-scale, asymmetric interaction between a single predator bird 
and a flock of smaller individuals. 
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Chapter 5 

Numerical scheme 



In this section, we propose a numerical scheme for finding solutions of the time-discrete model 



of Section 4.4 See Definition 4.4.6 The scheme was originally developed in [17]. 



We consider the set J] C such that J7 := [0, C] x [0, W], where £, W € (0, oo). We restrict 
ourselves to the situation in which each component of the crowd is either discrete (concen- 
trated in a finite number of points) or absolutely continuous (w.r.t. A"^), as was motivated in 
Section 4.8 Note that we have seen in the proof of Theorem |4. 5.1 that the push forward of a 
discrete measure is again discrete (and similar for the push forward of an absolutely continu- 
ous measure). The mass measure of each constituent is thus of the same type throughout time. 



5.1 Types of measures 

5.1.1 Discrete measure 

Suppose that we have a discrete mass measure for constituent a, consisting of A^" distinct 
Dirac distributions. We choose the particular form 

1=1 

for the mass measure at time slice n. Here, M" £ M"*" is a proportionality constant that 
makes it possible to compare a sum of Dirac measures (they measure the number of people) 
to an absolutely continuous measure (that measures kilograms of people). By {xf^}fj^^ C O 
we denote the set of (time-dependent) centres at time t = tn- 

5.1.2 Absolutely continuous measure: spatial approximation of density 

If the constituent a has a corresponding density p^, then we need to approximate in space. 
Therefore, we fix Nc,Nyy G N and we subdivide VL into Nc ■ Ny^; cells. Define he ■= C/Nc 
and /lyv := W/A'w, the horizontal and vertical grid size, respectively. The size of a cell is thus 
hchy\!. Each of the cells Ej^k is given an index (j, /c) G /C := {l, 2, . . . , Nc] x {l, 2, . . . , A^w}, 
such that 

Ej,k ■■= [{j - l)hc,jhc] X [{k - l)hn>,khw]- 
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Note that the boundaries of these cehs overlap each other. This is, however, not a serious 
problem, since the boundaries are null sets w.r.t. A'^. Up to a null set the cells are mutually 
disjoint. 



We sketch such grid in Figure 5.1 



El,Nv 



Enc .Nv 



E. 



hy\! 



ENr 



he 



C 



Figure 5.1: The spatial discretization of the domain by subdivision into Nc ■ -/Vyw cells. 

We approximate the density by a piecewise constant function p^. This means that for 
each (j, k) £ IC there is a ,^ G such that 

for all X in the interior of Ej j^.. Since the boundaries of the cells form a null set, it is not 
important how we define there. 



5.2 Calculating velocities 

In Section [46] we have proposed time-discrete velocity fields v" defined by 



<es(^) + E / f^{\y-^Mo:y)j^M{y), 
a ^ J \y 



foraUxGl^. (5.2.1) 



We approximate by v'^. Where and how an approximation is needed, depends on the types 



of measures associated to a and /3. We discuss this aspect in more detail in Sections 5.2.1 
and [5X1 



5.2.1 Discrete /z" 

If /i" is a discrete measure, then the evaluation of f " is required in each of the points xf „ G 0. 
The desired velocity v"^^^ can be evaluated in these points without any difficulty. However, 
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attention has to be paid regarding the integral terms with respect to /in that arise in the 
social velocity part of (5.2.1). We distinguish between the cases: 

(i) /in is discrete; 

(ii) /in is absolutely continuous. 
Firstly, let for each n be given as 

1=1 

for some G M+, N^^ € N and a set of centres {x^^f^-^ C 0. The corresponding interaction 
term can be calculated exactly as follows: 

/ /"(|y-:r|)5(e^,)^^d/.^(y)=M^ /^"(|xf„-rr|Mr , -. (5.2.2) 

J \y-X\ ■^■^^,n x 

n\{x} 

Any of the points can be taken as a choice of x. 

Remark 5.2.1. Note that the exclusion of {x} from the domain of integration is intended 
to avoid interaction of a point mass with itself for the case a = (3. However, in the right- 
hand side of (5.2.2) the contribution of interactions are also excluded if the positions of two 
distinct point masses coincide. That is, if we have xf ,^ = Xj^ for a = 13 but i ^ j (two distinct 

points from the same subpopulation), or xf^ = Xj^ where a f3 (two coinciding points from 
different subpopulations). We require our velocity field and initial conditions to be such that 



neither of these two situations occur; see also Remark 4.5.15 



Secondly, we consider the interaction integral for the case that /in is absolutely continuous 
w.r.t. A'^. The approximation pn of the density is used to obtain 

f^{\y-xM9-y)f^dp^^{y) = [ f;^{\y-x\)g{e^^)f^p^MdX'{y) 
'" x\ J \y x\ 

f^{\y-x\)g{6^y)f^/MdX'{y) 

\y x\ 



Q\{x} 



(5.2.3) 



The integrals over Ej^k are approximated by a two-dimensional form of the trapezoid rule, 
using the four vertices of the rectangle (Newton-Cotes). Let these vertices be called y*^, for 
i G {1,2,3,4}. For example, if yjj^ is the top right vertex, then := {jhc,khy\;). 
We obtain 



/ 



^ " \y — x\ 4 ."^^ ^ -^yj.k \y- , — x\ 
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However, if x = ^ for some i, we need to adapt this approximation, since typically fo has a 



/3 

singularity at 0. We then choose to use the interaction with the midpoint of the cell instead. 

yj"k •~ (0 ~ \)hc, {k — |)^w) denote the midpoint of Ej^k- Consequently, we use the 
approximation 



E 



(m) 



X 



(m) 



(5.2.5) 



for any x coinciding with a vertex of the cell Ej^k. 



Remark 5.2.2. If the position x coincides with a vertex of a cell, this causes a problem solely 
from the numerical point of view. The general scheme (trapezoid rule using the four vertices) 
is in that case no longer applicable for approximating the value of the integral over that cell. 
The concerning singularity does not cause a problem from the perspective of mathematical 



analysis. Namely, if we demand /^(s) 
bounded, then 

\y x\ 

n\{x} 



1/s, and assume that the density is uniformly 

f'^{\y-x\)g{e^^y 
f^{\y-x\) g{e% 




\P'n\ 



f1{\y - x\) d\\y) 
rp(\y-x\)A\\y). 



pi{y) d^y) 
(5.2.6) 



B{x,R) 



Here, B[x^ R) is the unit ball in M? around x with radius R := + W-^. Due to the choice 
of this specific radius it is guaranteed that C B(x,R) for each x £ Let us now only 
take the repulsive part around x into consideration. The contribution of the attraction-part 



(if present, cf. Section 4.2 and Figure 4.1) is finite any way. 



Let Rr denote the radius of repulsion. We have 



/II / 

Pn Woo 1 


fl{\y-A) 


B{x,Rr) 





Rr 



Mpi\ 



27r||/ 



nlloo 



Rr 

r 



rdr 



1 \ rdr 



< 



T^KWPnl 
OO. 



(5.2.7) 



We have taken the repulsive part of in its most simple form, but yet satisfying /^(s) ~ l/'^j^ 
A multiplicative constant can of course be added without loosing finiteness of the integral. 

^In fact, it suffices to impose the weaker demand that /^(s) ~ s ' for any 7 < 2 
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Remark 5.2.3. Contrary to what is stated in Remark 5.2.1, here we do not wish to forbid 



the situation that a point mass coincides with the vertex of a ceU. If two distinct point masses 



are located in the same position (the case considered in Remark 5.2.1), this corresponds to a 
physically impossible situation. Although probably only rarely a point mass will be located 
exactly on a cell's vertex, this is not physically undesirable, and we do want to allow for it. 



For those situations, we have therefore introduced an alternative approximation (5.2.5) to 
circumvent the occurring problems in the numerical scheme. 



5.2.2 Absolutely continuous fx" 

If //" is absolutely continuous with respect to A'^, then is approximated by a piecewise 
constant function. Let ^"(yj™^) be (an approximation of) the velocity of subpopulation a in 
the midpoint of cell Ej^f.. Then for each x in the interior of Ej^i^, the approximated velocity 
field Vni^) is given by 



In each midpoint yj™^ evaluating the desired velocity v2q<, does not cause any problem. In 
order to deal with the integral terms in the social velocity, we again take into consideration 



which type of measure /in is 



Firstly, we treat the case when /i^ is a discrete measure. We specifically suppose that /i^ 
has the form 

Ar/3 

' l.TL 

1=1 



If, for a given cell Ej^k, the midpoint yj^ does not coincide with any of the centre points xf^^, 
then 

(m) NP „(m) 



i,n »j 



^^•"^ \y-yj,k \ i=i K,n-yj,k \ 

(5.2. 



If by chance one of the centres, say Xmn, of a Dirac mass coincides with the midpoint yfl! , 
then we have to adapt our scheme. We replace the direct interaction between these two 
points, by an average over the interaction between Xm,n and each of the four vertices ^ of 
the cell. That is, we use 

-j:f^{\<n-yiMo:. ) ;r"ti - 

i=l ^' \Xm,n Vj^kl 
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Using this expression, we obtain the fohowing approximation: 

(m) 

V' ,'V'\ (™ 



I 



(m) 



i=l 



13 _ (m) I 



yj,k ''^,r^ \x'] „ - W ^ ^ 



Xr 



yj,k 



(5.2.9) 



i=l \Xm,n lljkl 

We have merely overcome problems that might occur in particular if the spatial grid is too 



coarse. For a sufficiently fine grid, the condition in Remark 4.5.15 makes sure that there is 
a neighbourhood of zero density around the position of each Dirac mass. It is not important 
whether the velocity is properly defined in regions of zero density, or not. 



We look now to the case when is an absolutely continuous measure. We follow (|5.2.3) 



and (5.2.4), and derive 



_ (m) 

J yj,k y \y — Vj ' 



4 i (™) 

Pu,k),n 2^ Jp y\yj,k - yj,k \)9(>^jn.) 



ij,k)eJC i=l 
,{m) 



'^'^'^^''\ylk-y^\ 



(5.2.10) 



Keep in mind that the midpoint ^ can never coincide with any of the vertices of a cell in 

n. 



5.2.3 Summary: The approximation t;^ 



We now shortly summarize the way we have defined in Sections 5.2.1 and 5.2.2 



If /i" is discrete, we need the value of in the corresponding centres xf^ of the Dirac 
masses. 

Calculation of the value of v^^^ is straightforward. Those terms in the social velocity 
for which the measure Hn is discrete are calculated according to (5.2.2). If, on the other 



hand, the measure is absolutely continuous, then the approximation of (5.2.3) is 
used. Regarding the approximation of the integrals over all cells Ej^^, we distinguish 
between two cases: 



does not coincide with one of the four vertices of cell 



Use the two- 



dimensional trapezoid rule as given by (5.2.4) 



— xf^ coincides with one of the four vertices of cell Ej^k- Use the midpoint of the cell 
to approximate the integral; see (5.2.5). 
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If is absolutely continuous, the velocity in each interior point of a cell Ej^k is approx- 
imated by the velocity in the midpoint of that cell. Evaluation of v^^.^ the midpoint 
yj™^ is again straightforward. If the measure in the social velocity is discrete, then 
we distinguish between two cases: 



None of the Dirac-centres coincides with the midpoint yj™^- Then use (5.2.8) 
to evaluate the integral. 

One of the centres xf^ is located exactly in the midpoint Vj^^- Use the approxi- 
mation given in (5.2.9). 



If is absolutely continuous, then the approximation as given by (|5. 2.101). 



5.3 Push forward of the mass measures 

The only thing that is still to be included in our numerical scheme, is the push forward of the 



measure /u" to obtain We use (see Section 5.2, in particular Section 5.2.3) to obtain 

a modified version of the motion mapping (cf. Definition 4.4.1): 



If /i" is a discrete measure, the push forward can be found in a natural way, as was already 



suggested in Corollary 4.5.8 If /i" is given by 



1=1 



then the numerical approximation of the push forward is 



N° 



i=l 



Determining the push forward of this discrete mass measure boils down to updating the 
centres of the Dirac masses in the following way: 



X, 



i,n+l • 



If /i" is absolutely continuous, then a little more effort is required. Since is defined such 
that it is constant within a cell, the push forward of a cell Ej^^ is a translation by the vector 

/S.tnVniv'j^)- Indeed we have 



■j,k) 



{x'^:ix)\xGEj^k} 

{x + AtnVn{x)\x £ Ej^k} 

{x + Atnv:{y\f)\xGE,,k} 



The push forward of Ej ^ is indicated in Figure |5.2[ We want to keep the spatial grid fixed 
in time. Also, we require our approximated density p'^ to be constant within a cell, for every 
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Figure 5.2: The push forward of a ceh Ej^^ governed by the approximated velocity field. The 
'image' Xn{Ej^k) clearly lies within four cells of the spatial domain. 



n. In general, the push forward of Ej^k will not exactly coincide with a cell of our grid: cf. 
Figure 5.2, where XniEj,k) lies in four cells. We thus propose the following update of the 
density: 



l(^) = P'{j,k),n+l 



XHE. 



(p,<?)e/c 



(5.3.1) 



for all X in the interior of -Ej^fc- 

A positive contribution is given by each cell Ep^q that is (partially) mapped into Ej^k, since 
only then Ej^^ H Xn(-^p,g) is non-empty, and thus has positive Lebesgue measure. Note that 
P'ipq) n^'^i-^j^i^ ^ Xn{Ep,q)) is exactly the mass that is transferred from Ep^q into Ej^^- The 
density of Ej^^ then follows from adding the mass contributions of all cells Ep^q, and dividing 
by the area of Ej^^- This makes (5.3.1) a natural way of updating densities. 



The numerical approximation of the push forward of is fully determined by the itera- 
tive scheme for its density p° as given by (5.3.1). 



Remark 5.3.1. Using the update of the density in (5.3.1), Pn+i is by definition non-negative, 



as it is a sum of non-negative terms. Furthermore, the total mass contained in Q is conserved 
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in time. We namely have that 



/ 



(i,fc)e/c 

{j,k)eK. ^ {p,q)eK. 

(p,g)e/c (i,fc)e/c 

(p,9)e/c 



We have used in the fifth step that Xn is a translation, due to which X'^(^Xn{^p,q)) — ^'^{^p,q)- 
Conservation of total mass in $7 is only guaranteed under the assumption that no cell (of non- 



zero density) is mapped outside 0; cf. Part (m) of Assumption 4.5.13 
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Chapter 6 



Numerical illustration: simulation 
results 



In this section we present the results of our simulation experiments. Except for some simple 
test cases, we only consider two-scale situations. That is, in each experiment there are two 
subpopulations, one of which is discrete and the other is absolutely continuous. The compo- 
nent indexed a = 1 is a collection of discrete individuals. In our simulations, component 1 
consists mostly of one individual only, and exceptionally of two. The component with index 
a = 2 is a macroscopic crowd. 

To avoid effects at the boundaries (which we have not specified so far), the domain Q is 
'sufficiently large'. That is, the size of is such that no mass reaches the boundary. No 
problems occur as long as the cells of the spatial grid on the periphery of have zero density, 
and the positions of the individuals remain in fij^ 

In the following sections, the results of the simulation are presented by giving a graphical 
representation of the crowd in the domain. The positions of individuals are marked by red 
bullets, whereas the density of the macroscopic crowd is indicated by a grey shading. The 
darker the colour, the higher the density; a colour bar at the right-hand side of the graph 
shows what shading corresponds to a certain density. 

Note that we essentially use the scheme of |17[ 145) . The convergence and stability of the 
numerical solution is proven in these papers. We expect similar properties to hold for our 
setting. We, however, omit to give further details in this direction and focus directly on 
simulation results. 

The results we present will be interpreted only on a qualitative basis. At a later stage we 
will try to make some of these results quantitative by recovering experimental data by S. 
Hoogendoorn, W. Daamen and M. Campanella (Delft University of Technology, see Section 



1.4 of the Introduction) regarding experiments in a corridor. 



We admit that this is a rather pragmatic solution, but for the moment it is the best we can do. 
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6.1 Reference setting 



For attraction-repulsion interactions, we use the following expression in the interaction inte- 
gral: 



pAR 



R 



AR 



.i?AR 



i?AR(i^AR _ ^AR) 
0, 



if < s ^ R^^ ; 



{s-R^^){s-R^^), if i?AR < ^ ^ ^AR. (6.1.1) 



if s > R^^. 



The radii R^^ and R^^ are such that < R^^ < R^ . The factor F^^ G [0, oo) is a positive 
scaling constant. Note that this is exactly the function plotted in Figure 4.1 (left). Note that 
is differentiable in s = i?^^. 

Similarly, we define for purely repulsive interaction 

R?- 



if s > i??. 



(6.1.2) 







Again, we take R^ > and the scaling constant S [0,oo). This function is plotted in 
Figure 4.1 (right). 

Unless indicated otherwise, we use a set of 'standard' parameters in our simulations. As 
dimensions of the domain Q we take 

C = W = 50. 

The proportionality constant is only present for a = 1 and is assigned the value 

= 60, for a = 1. 

The desired velocity is independent of the space variable, and only the direction is different 
for the two subpopulations. We take 



^des(^) 



— 1.34ei, for all x € il, 
1.34 ei, for all x G $7. 



Here ei is the unit vector in the direction that corresponds to the horizontal axis in our 
graphical representation. 

Regarding the interaction functions, we choose the following reference parameters: 



i?AR 


= 0.03, 




= 0.03, 


pAR 
rLj. 


= 1.5, 


pAR 


= 3, 


pR 


= 4, 


a 


= 0.5. 
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Attraction-repulsion interactions take place within one subpopulation, whereas repulsion-only 
takes place in the interactions between distinct subpopulations. 

In the sequel we indicate clearly where we deviate from these standard parameters. 

6.2 Two basic critical examples: one micro, one macro 

We investigate here whether the simulation results are mathematically and physically accept- 
able. 

Consider two individuals both having desired velocity in the direction of ei. The magni- 
tude of the desired velocities is equal and constant, but they are directed oppositely. The 
initial positions of the two individuals differ only in the ei direction, and they are located 
such that they initially want to move towards each other due to their desired velocities. Their 
interaction is purely of repulsive nature, and has the same parameter values for any of the 
two. We expect these individuals to approach each other, up to a certain distance. At this 
distance the desired velocity in one direction is in balance with the (oppositely directed) re- 
pulsive effect in the social velocity. The total velocity is thus zero for both individuals. 
Let Xi and Xj denote the positions of the two pedestrians. For simplicity we take a = 1 and 
M* = = 1. The velocity of pedestrian i is given by 



v{Xi) = -y^es + f^{\Xj - Xi 



(6.2.1) 



with as in (6.1.2). The velocity is zero if 



-r{\xj-x. 



A necessary condition is 



/ {\Xj Xi 



(6.2.2) 



Note that {xj — Xi)/\ unit vector. Considering (|6.1.2l), condition (|6.2.2l) reads 



1 =Ke 



if \xj — Xi\ ^ B^. This yields that 



at the point where the velocity is zero. 



^des 



+ i?R 



(6.2.3) 



In Figure 6.1 we show the outcome of this simulation. The individual placed initially on 
the left has desired velocity v^cs = 1.34 ei, while the one on the right has desired velocity 
v^cs = — 1.34ei. Furthermore we take = 1. In Figure O the mutual distance is plotted 
against the time. Indeed the equilibrium distance is the one predicted by (6.2.3). 
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Figure 6.1: Two individuals approaching each other until a certain minimal distance is 
reached. The images are taken at t = (left), t = 15 (right). 




Time 



Figure 6.2: The mutual distance between two approaching individuals (blue). The predicted 
minimal distance is indicated in red. 



Note that the setting of Figure 6.1 is unrealistic and thus undesirable. In everyday life 
these two pedestrians would namely both move a bit aside and then walk straight ahead 



without any constraints. In Section 1.3 of the Introduction, we have indicated already that 
a small amount of random noise can be used to avoid these deadlocks. We expect that the 
equilibrium configuration is instable, and that deadlocks do neither occur if we perturb the 
initial data (in the direction perpendicular to the desired velocities). 



In the experiment presented in Figure 6.3 we observe that a macroscopic crowd tends to 
form a circular configuration. The crowd has no desired velocity: f^es = 0, and therefore 
(T = 1 is taken. This typical behaviour is well-known from molecular dynamics. The mutual 
interactions favour the minimization of the ratio circumference to area. A mathematical proof 
is given in [21] , for the same macroscopic (continuous-in-time) equation of mass conservation. 
The interaction potential in [23] is not the same as the one used here, but we expect that 
similar results can be derived along comparable lines of argument. 



6.3 Two-scale interactions of repulsive nature 



In Figure 6.4, we show the interaction between a macroscopic crowd and an individual that 



wishes to approach it, and eventually forces itself a way through. This is a situation in which 



we try to mimic the two-scale ('predator-prey') behaviour described in Section 4.8 



A circular area of (nearly) zero density is formed around the individual. The size of this 
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Figure 6.3: The simulation of a macroscopic crowd's motion. Initially, the crowd forms a 
square of uniform density p = 2. Its configuration evolves into a circular shape. The images 
are taken at i = (top left), t = 30 (top right), t = 60 (bottom left). 



region typically depends on B^. At a distance shorter than mass is driven away from 
the individual. If we decrease the radius B^, also the 'empty zone' around the individual 
decreases; see Figure 6^ for this effect. We estimate the distance in front of the individual 
that is empty, and use the ideas of Section 6^ to do so. Let us consider the point x in 
the crowd that is right in front of the individual. Disregarding the effect of the rest of the 
macroscopic crowd, the velocity of x is similar to the one given in (6.3.1): 



v{x)=vl,, + f'^{\z-x\) 



Here, the variable z is used to denote the position of the individual. The desired velocity 
of the individual v\^^ is coupled to v\^^ via the relation: v\^^ = — f^es- Assume that the 
individual moves at this speec^and that x remains right in front of z. Then the mass located 

^This is generally not the case. 
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Figure 6.4: The interaction between a discrete individual and a macroscopic crowd, 
images are taken at t = (top left), t = 5 (top right), t = 10 (bottom left). 



The 



in X is forced to move also with velocity v^^^. This leads to the following equation: 

v{x) = vl, + f^^ilz - x|)^— ^ = vl, = -vl,. (6.3.1) 



Following the lines of Section 6.2 we derive that a necessary condition is that 



X 



2\vlJ+MF^ 



(6.3.2) 



According to this estimate the size of the empty area around the individual scales linearly 
with (note that |z — x| ^ R^). In Figure 6.5 the configuration is given for two distinct 



values of Rf-. Indeed the size of the low density zone decreases as Rf^ decreases. The distance 



estimate in (6.3.2) is indicated by a circle in blue. Right in front of the individual there is a 



spot of high density. This is the spot x we were considering. We observe that x is located just 
outside the blue circle, by which the estimate turns out to be quite good. We also observe 
that the estimate loses its value more to the sides of the individual. After the individual 



passes by, we observe in Figure 6.4 that the crowd splits in two and that the two halves do 
not come together again. This effect depends on the choice of parameters. By adjusting the 
parameters we can achieve the crowd to 'envelope' around the individual. To this end, we 
decrease the radius of repulsion of the interaction between the individual and the crowd. Also 
we increase the radius of attraction of the internal interactions in the crowd. To do so, we 
take R^ = 2, R^^ = 5. This modification makes that the macroscopic crowd does enclose 



the individual. The result is shown in Figure 6.6 
We deduce from Figure 



6.5 



that the width of the empty zone is of the order of 2R^ (probably 

such that the 



a little bit narrower, as also (6.3.2) suggests). Let us define a critical radius Rc 



width of the empty zone is 2i?cr- We expect the crowd to come together again, if -R^ ^ ^ 2i?c 
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Figure 6.5: The size of the empty zone around the individual depends on i?^. On the left 
= 4 (this is the top right situation in Figure 6.4), on the right = 2. In blue a circle is 



drawn with radius equal to the distance predicted in (6.3.2). The images are both taken at 
t = 5 starting from identical initial configurations. 



15 20 25 30 35 40 45 50 



25 30 35 40 



Figure 6.6: The interaction between a discrete individual and a macroscopic crowd. Modi- 
fication of the radii of interaction makes that the crowd comes together after the individual 
has passed by. The images are taken at t = (top left), i = 4 (top right), t = 5 (bottom left). 



this is the situation in which the crowd in one half can also 'feel' the part on the other side 
of the empty zone. 



Remark 6.3.1. We have tried to achieve the effect of Figure 6.6 also by adding macroscopic 
mass above and below the initial square. That is, initially the crowd now has a rectangular 
shape twice as wide (in vertical direction in the graph). We hoped that the repulsive effects 
within the crowd would force the two halves to move towards each other again. However, this 
effect was not obtained. 

Remark 6.3.2. The desired velocity of the individual does not need to be necessarily in the 
direction of ei. If the desired velocities of the individual and the macroscopic crowd are at 
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Figure 6.7: The interaction between two discrete individuals and a macroscopic crowd. The 
individuals are initially a distance of apart. The images are taken at t = (top left), 
t = 5 (top right), t = 10 (bottom left). 



6.4 Interaction with two individuals 

Two-scale simulation experiments can involve two individuals instead of one. The outcome of 



such a simulation can be seen in Figure 6.7 The two individuals initially are at a relatively 
short distance from one another. Their mutual distance remains small and macroscopic mass 
is not able to separate the two. The macroscopic crowd only moves around them. In fact, 
the impact of the two individuals is not fundamentally different from the impact a single 
individual of double mass would have. 

If we increase the initial distance between the two individuals, then the macroscopic mass 
does find a way to fill the space between them. It even forces the individuals more apart. 



This can be seen in Figure 6.8 Note that it still requires some 'effort' to separate the two 



individuals. In the top right image in Figure 6.8 a region of high density is present at the 



left-hand side of the opening between the individuals. We expect that it depends on the 
initial distance between the two individuals compared to R^, whether macroscopic mass can 
pass through the opening between the two individuals or not. 
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Figure 6.8: The interaction between two discrete individuals and a macroscopic crowd. Ini- 
tially the distance between the individuals is slightly more than Ra^, and thus they do not 
'feel' each other. The images are taken at t = (top left), t = 5 (top right), t = 10 (bottom 
left). 



6.5 Modelling leadership 

Can we capture 'leadership' effects via our two-scale model? More precisely, can an individ- 
ual be enabled to drag macroscopic mass along? We interpret such behaviour as a tendency 
to follow the individual, even though the macroscopic crowd in itself has no special desire 
to go in a certain direction. The individual then acts as a leader (or guide) for the group. 
To achieve this effect, in this section, we choose attraction-repulsion as the influence of the 
individual on the macroscopic crowd, instead of solely repulsive interactions. 



In Figure 6.9, the individual is only driven by its desired velocity, which is directed to the 



left. The crowd does not affect its motion. Yet, the dynamics of the macroscopic crowd 
is influenced by the 'target particle' due to attraction-repulsion. The crowd has no desired 
velocity, thus a = 1 is taken. All other conditions are as in the reference setting. Figure 



6.9 shows that the individual is able to create a short slipstream of macroscopic mass, when 



moving through the crowd, but we can not convincingly call it a leader. 



The striking feature appearing in Figure 6.9 are the vertical lines, alternatingly of high 



and low density. These appear in the macroscopic crowd. Why does this happen? As Figure 
6.10| shows, these oscillations are apparently only a result of the numerical approximation. 
The density is smoothened if we decrease the grid size. 



We wish to enhance the effect of leadership (compared to Figure 6.9) by adjusting parameters. 
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Figure 6.9: The interaction between a discrete individual and a macroscopic crowd. The 
attraction-repulsion influence of the individual on the macroscopic crowd makes that the 
crowd is slightly dragged along. A small amount of mass evades on the left-hand side of the 
crowd. The images are taken at t = (top left), t = 7 (top right), f = 14 (bottom left). 



Firstly, we increase the magnitude of the attraction-repulsion influence of the individual on 
the crowd. More specifically, we take = 0.3 instead of = 0.03; for interactions 



within the crowd F^^ = 0.03 is maintained. All other conditions are as in Figure 6.9 In 



Figure 6.11 we see the effect of the performed parameter modification. 

Note that within the environment of the individual, the density attains very high values, 
which is physically not acceptable. Let us concentrate on qualitative behaviour only. The 
individual successfully attracts a large part of the macroscopic crowd, at least initially. As 
time elapses, it loses control over the crowd. Maybe, this fact is due to his too high velocity. 
At first, the mass in its direct environment is distributed more or less uniformly over a circle 
around it; this mass shifts more and more to its rear end. Also, we observe a trace of mass 
left behind by the individual. Once this mass is at a distance greater than from the 
individual, it remains near the place where it loses contact 

Effects of a similar kind can be obtained not only by modifying the interaction strength, 
but also by increasing the radius of attraction. In the attraction-repulsion infiuence of the 
individual on the crowd we restore the value F^^ = 0.03, and take = 6. Indeed this 



enables the individual to act as a leader, cf. Figure 6.12 



The macroscopic crowd is initially attracted to the individual, which is comparable to Figure 



6.11[ This situation differs from Figure 6.11 in the sense that the macroscopic crowd is much 
less compressed. A possible explanation is twofold. The area of attraction around the leader 
is much larger, and moreover, the attraction towards the centre of this region is less strong. 



''The scaling of the grey shading has been altered in order to be able to distinguish lower densities from 
zero density. Very high densities occurring in small regions dominate the scaling. The aforementioned trace 
of mass is therefore hardly visible, in spite of the fact the density is around p ~ 2. 
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Figure 6.10: Decreasing mesh size for the setting of Figure 6.9 at time t = 14. Top left: the 



original grid size of Figure 6.9 Top right: the number of cells has been increased by a factor 
1.5 in both horizontal and vertical direction. Bottom left: increase of the number of cells by 
a factor 2 compared to the top left situation. Oscillations in the density are diminished by 
decreasing grid size. 



Note that due to the weaker attraction-repulsion (compared to the previous case = 0.3), 
the leader is less successful in keeping his followers with him. Finally, the macroscopic mass 
has moved on average a smaller distance to the left. However, a positive aspect of this setting 



is that the maximal density is nearly a factor 4 smaller than in Figure 6.11 
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Figure 6.11: The interaction between a discrete individual and a macroscopic crowd. The 
influence of the individual on the crowd has been increased, and as a result the individual is 
able to drag the crowd along. The images are taken at t = (top left), t = 7 (top right), 
t = 14 (bottom left). 
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Figure 6.12: The interaction between a discrete individual and a macroscopic crowd. The 
radius of attraction has been increased, which enables the individual (temporarily) to take 
macroscopic crowd along. The images are taken at t = (top left), t = 7 (top right), t = 14 
(bottom left). 
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Chapter 7 



Discussion 



At this point, we wish to look back and hst a few open issues with respect to modelhng, 
analysis and simulation. We will investigate some of them in the near future. 

• Via (mainly formal) calculations we arrived at a fully continuous-in-time model. Sub- 
sequently, we deduced a time-discrete version of this model, which we analyzed math- 
ematically. We proved global existence and uniqueness of a time-discrete solution (in 



Theorems 4.5.1-4.5.2), and showed some properties of this solution (e.g. positivity and 
conservation of mass). A number of conditions on the (one-step) motion mappings and 
velocity fields were inevitably needed. Afterwards, we proposed a specific form of the 
velocity field, in which a dependence on the instantaneous configuration of the mass 
appears. Up to now, we have not proven that the chosen velocity field actually satisfies 
the imposed conditions. This is still to be done. We expect this step to be difficult, 
mainly due to the functional dependence of the velocity on the mass measures, and due 
to the fact that the dynamics of all subpopulations are coupled via their velocity fields. 

Formulating and proving results like existence and uniqueness of solutions to the continuous- 
in-time problem, is our next goal. However, this will be a much more difficult task than 
treating the discrete-in-time model. For instance, the one-step motion mapping of the 



mass measures is linear in r;" for the time-discrete model, cf. Definition 4.4.1 This 
is because we linearized by using Taylor series approximations in order to obtain the 
time-discrete model. The motion mapping will in general no longer be linear in the 
velocity if we consider the continuous case. Moreover, it is all but clear whether we 
can even speak so easily about a one-step motion mapping and the accompanying push 
forward of the mass measure. These concepts are likely to translate into continuous-in- 
time counterparts, relating the configurations at time instances that are not an a priori 
fixed time step apart. 

We would like to enlarge the class of interaction potentials for which an entropy in- 
equality can be derived. For the moment the potential is only admissible if 
Wj^(^) = W^(— .^) for all ^ G M'^, and interactions are symmetric. The latter con- 
dition means that = Wa- See Sections 4.3 and 

These restrictions on are simply to strong for many interesting settings. An idea 
at least to circumvent the condition W?(^) = Wg{—^) is inspired by [lOj (especially 



4.7 



and especially Assumption 
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Section 3.1 therein). We write Wj^ as the superposition of a symmetric part and a drift 
part: 

yy f), symm ' drift; 

where symm(0 = 3y^„,(-e) is satisfied for all ^ G W^. The second term W° ^^-^^ 
contains the deviations from symmetry. We are interested to see which restrictions on 
Wp are needed in order still to be able to derive an entropy inequality. 
However, some problems arise in the interpretation of such interactions. Consider the 
vector ^ in its representation in polar coordinates. Without the condition that VF^(C) = 
Wp{—^,), the interaction potential certainly depends not only on the length of ^ 
but also on its angle (the azimuth). Since the velocity field involves VW^, it inevitably 
has a non-zero azimuthal component. We should ask ourselves the question whether 
this is physically acceptable: are the interactions allowed to influence the velocity in a 
direction that is not parallel to the connection vector between two points? 
If we do not want to restrict ourselves to a gradient (in this case VW^), we could replace 
it by a vector field of the form VVF + V x (Helmholtz decomposition). Under which 
conditions on the vector field uj can we still derive an entropy inequality? 

• A wider range of domains can be considered, if we know how to incorporate walls (outer 
boundaries) and internal obstacles in our domain. We have not treated this aspect so 
far, but it surely is very important to do so. In [H], Piccoli proposes to project the 
obtained velocity on the space of admissible velocities. That is, those velocities that 
satisfy v ■ n ^ at the boundaries. In practice, this means that outward pointing 
normal components of the velocity are to be set zero. Perhaps, we would like to have 
an alternative way, which fits better to our measure-theoretical framework. Probably 
one needs to define boundary measures and Cauchy fiuxes, while also describing their 
throughput. 

• More fundamental mathematical challenges are related to discrete-to-continuum limits. 
We would like to investigate in what sense we can approximate the absolutely continu- 
ous part of our mass measures by particle systems (discrete mass measures in the case 
of crowd dynamics). Especially, the limiting process N ^ oo (where N is the number 
of Dirac measures in the particle system) is relevant. Such approximation, in combina- 
tion with the Wasserstein metric, might turn out to be useful for proving existence and 
uniqueness properties for the continuous-in-time model, and also to derive alternative 
numerical schemes for computing the density based on particle simulations. 

The question is which scaling one needs to obtain our macroscopic perspective in the 
limit. Moreover, we are interested in whether other types of scaling would lead to differ- 
ent limits, e.g. a mesoscopic Boltzmann-type description ('intermediate level' between 
micro and macro, see e.g. [15]) or the twofold micro-macro setting of [17]. Perhaps 
there even exists a way to exploit the fact that we have multiple subpopulations. That 
is, we might apply different scalings to distinct subpopulations, such that e.g. for one 
subpopulation the limit is of Boltzmann-type, and for another the limit is a two-scale 
micro-macro description. What is the actual interpretation of such distinction between 
the two subpopulations? 

• Regarding our simulation results shown in Section [6) We observe that we have success- 
fully circumvented problems that might occur at the boundary, by considering simula- 
tion runs as long as all mass is contained in the interior of 0. Such ad hoc solution can 
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be avoided in the future if we find appropriate tools for including boundary effects (see 
above). When outputting the configuration, our simulation tool is equipped to check 
whether the total mass is actually conserved. This turns out to be indeed the case for 
every instance of our simulation. 

When examining the results of our simulation, several questions arise. All of them are 
the subject for further research: 



Figure 6.3 indicates that the equilibrium configuration for attraction-repulsion in- 
teractions is a ball. If instead, we test the same situation for purely attractive 
interactions, does the crowd shrink to a point? Is there a mathematical way to 
show that the limit is a Dirac mass? Note that in this case the right-hand inequality 



in (4.5.1) can no longer be satisfied. 



— In Section 6.3, we estimate the radius of the (nearly) empty area in front of the 
individual. It might be possible to characterize the shape and dimensions of the 
rest of the empty zone induced by the individual. In the same section, we show 
that a modification of parameters makes the crowd clog together again after the 
individual has gone past. What are the typical length and time scales at which 
this reunion takes place? 

— The interaction between two individuals and a macroscopic crowd was simulated in 



Section 6.4 If we zoom out, under what conditions can we deduce from the macro- 
crowd's behaviour that more than one individual was present (cf. Figure 6.8)? 



When do we obtain merely the same results as in the presence of one individual of 



double mass (cf. Figure 6.7), and how does this effect depend on the initial mutual 



distance between the two individuals? 



In Section 6.5 we tried to create leaders. We observed (especially in Figures 6.11 
and 6.12) that the individual at first is able to take a part of the macroscopic 
crowd along, but that he loses control as time goes by. Regarding the long-time 
behaviour of the crowd, will eventually all mass be lost by the leader? When 



do leadership effects end? In Figure 6.11, bottom left, the 'tail' of mass behind 



the individual consists of two regions. Closest to the individual the density is 
relatively high, whereas after quite sudden transition the density is lower. What 
are the characteristics of both regions (shape, dimensions), and why is there an 
apparent strict separation between them? 

In all of the cases above: How exactly do the observed phenomena depend on the 
parameters? What can we say about the (in)stability of equilibrium configurations? 
Can we tune parameters in such a way that physically unacceptable densities do 
not occur? Can we identify a set of dimensionless quantities, that characterizes the 
dynamics in the model? Also the role of the initial conditions needs to be taken 
into consideration. 
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Appendix A 



Proof of Theorem 2.1.4 



This proof follows the lines indicated by [22], p. 42. 

Proof. 1. Define the set £ := {Q' € B{Q) | X{Q \ il') = 0}. Note that £ is non-empty since 
at least 17 G Choose a sequence {i?fc}fceN+ ^ such that 

fi(Bk) ^ inf n(n') + 1, for aU /c G N+. 

Such sequence exists by definition of the infimum. 

Define B := flfc^i-^fei ^ is an element of B{0,), since any cj-algebra is closed w.r.t. 
countable intersections. By one of De Morgan's laws we have that Q \ (HfcLi ^fc) = 
UfcLi \ Bk) , for which we should note that this is not a disjoint union. Therefore we 
have 

oo 

oi^ x{n\B) i^Y^ \{n \ Bk) = 0, 

k=l 

since \ Bk) = for all k G N"*". We have thus proven that B £ £. 
In fact we even have that niB) = inf^j/g^- To see this, note that 

oo ^ 

^{B) = //(^pl ^ fi{Bj) ^ inf + -, for all j G N+. 

k=l ^ 

Since this statement is true for any j G N"*', it follows that n{B) ^ inf^/g^: Due to 

the fact that B £ £, it is clear that also fJ-{B) ^ inf^/g^- thus = infj^/g^- /i(Q'). 

2. Let ;Uac : IB(J7) M+ be defined by /iac(J^') = /^(J^' n B), and similarly : IB(J]) M+ 
by fJ,s{^') = /^(f^' n (il \ -B)) . Note that it follows from this definition that /i = /iac + /^s- 

Let A G ^(17), be such that A C B and A(A) = 0. Assume that fj,{A) > 0. Note 
that B\Ae B{n), and 

x{n\{B\ A)) = \B)uA) = x{n \b) + x{A) = o. 

Here we use that Q \ B and A are disjoint (since A C B) and that A(il \B) = (since 
B e£). We conclude that B\Ae£. 

As a result of the fact that yl C -B, it holds that fj.{B) = fi{B \ A) + ^{B n A) = 
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fi{B\A)+ fj,{A), and because we assumed that ij-{A) > it fohows that fi{B\A) < n{B). 
As B \A £ £, this contradicts fi{B) = inf^/gf /x(O'). Thus fi{A) = must hold. 
Let n' e B{n) now be such that A(il') = 0. Then //ac(f^') = fJ-i^' n B) = fohows if we 
set A := ^l' n B C B in the lines of arguments above. Hence, //ac ^ 

Furthermore fis{B) = PI (il \ B)) = ^(0) = 0, and, due to the fact that B € £, 
X{Q \B) = 0. This implies /^s -L A. 

3. To prove uniqueness of this decomposition, assume that we have fi = fi}^^ + /ig and 
/X = /x^j, + fig, where ^ _L A for i G {1, 2}. For all Q' G B{Q,) we thus have 

Neither in the left-hand side, nor in the right-hand side we necessarily have a positive 
measure. However it can quite easily be seen that (/i^^, — /i^^) <^ A, because X{Q') = 
implies = for i = 1, 2 (as fi\^ <^ A), and thus also (/u^c — /"ac)(^') = 0- 

Furthermore there exist Bi,B2 G I3{fl), such that nK^l \ Bi) = \{Bi) = for all 
i G {1, 2} and for ah Q' G B{^). Define B := BiU B2, and note that n\B C n\Bi 
for each i G {1, 2}. It follows that ^ \ -B) ^ fil{n \ Bi) = 0, for i G {1, 2}, thus 
=0 and thus 

{fil-f,l){n\B) = 0. 
Also, ^ X{B) = X{Bi) + X{B2) = 0, and thus 

X{B) = 0, 

from which we conclude that — fil) ± X. 



We now have jl <^ X and /i _L A, and will show that this implies jl = 0. Use statement 
^ from Lemma 2.1.3 to characterize singular measures, by which we have disjoint 
Ai, G B{n) such that i2{n') = n Ai) and X{n') = A(0' n A2) for ah J^' G B{n). 
For any Jl' G B{n), we have A(J7') = A(f]' \ A2) + A(l^' n A2) = X{n' \ A2) + X{n'), 
which implies A(r2' \ A2) = 0. Because /i is absolutely continuous w.r.t. A, we also have 

li{n' \ A2) = 0. Thus: = ji{n' \ A2) = \ A2) n ^1) = n ^1) = 

Since fl{Q') = for all Q,' G B{Q,), we have that ft = and hence //^c = fJ-'ic ^^'^ A's = Ms- 

□ 
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Appendix B 



Proof of Theorem 2.1.9 



The inspiration for this proof comes from [37], pp. 45-46, although the theorem presented 
here is more general than the one stated in [37] . 



Proof. 1. Due to the fact that /Xg is singular w.r.t. A, we are provided two disjoint sets 
Ai,A2 G B{n) such that = fisi^'nAi) and X{n') = A(0'n^2) for ah Cl' G B{n). 

Now define the sequence of sets {Bn}n(^n+ given by 

Bi := {x G Ai |/xs(x) ^ 1}, 

Bn := {x e Ai\- Hs{x) < ^—}, for n G {2,3, ...}. 
' n n — 1 

Each Bn contains only finitely many elements, because is a finite measure. It follows 

that B := IJ^i is a countable set. Note that B = {x ^ Ai\ ij,s{x) > 0}. Also note 
that X ^ B implies /Us(x) = 0. Indeed, ii x £ Ai\ B this is trivial. Furthermore, if 
X ^ Ai then /Xs(a;) = fis{{x} Ai) = figi^) = 0. 

2. Let /id : M+ be defined by Hd{^') = ^s(f^' n B), and similarly fi^c : B(fi) M+ 

by ^sc{^') = Ais(^^'n(r2\i?)). Note that it follows from this definition that fis = fid + fJ-sc- 

It is readily seen that /id is discrete w.r.t. A, because i? is a countable set of points in 
0,, for which we have 

/id(J^ \ B) = /id((0 \B)nB)= fxam = 0, 

and (since B C Ai) 

X{B) = X{B n A2) ^ X{Ai n A2) = A(0) = 0. 

We now show that /isc is singular continuous w.r.t. A. Let x G be arbitrary. fisc{x) = 
fMs{{x} n (17 \ B)). Thus, if X ^ Q\B then fisc{x) = /isc(0) = 0. On the other hand, if 
X £ Q\B then /isc(x) = ij,s{x) = because x ^ B (the last step has been shown above). 
Thus, fJ-sc{x) = for all x G il. 

Consider the set Ai. First of ah, A(^i) = X{Ai n A2) = A(0) = 0. Moreover, fi^d^ \ 
Ai) = fis{{^ \ Ai) n {i^ \ B)) = /is(r2 \ j4i), because B C Ai. Since /is is singular, we 
have fis{^ \ Ai) = /is((0 \ Ai) n Ai) = /is(0) = 0. 
We now have that /isc is singular continuous w.r.t. A. 
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3. To prove uniqueness of this decomposition, assume that we have Hs = fJ-\ + Msc 
)Us = /x^ + //gc) where is discrete w.r.t. A and fil^ is singular continuous w.r.t. A for 
i e {1, 2}. For all n' G B{n) we thus have 

{Hl - f^lm = {nl - nlm') =: A(O'). 

Assume that Ai and A2 are the countable sets, such that Md(^ \ ^i) — K^i) = 
for i G {1,2}. Define A := AiU A2, which is obviously countable. We have for each 
i e {1,2}, and all Q e B{n) 

n A) < i4{n') = ni{n' nA) + \ A) < ^4^{n' n A) + ni{n\ Ai) = ni^in' n A), 

and thus l4^{^') = Mdl^' ^ This implies /i(O') = n A). Note that Q' n A is 
a subset of A and is thus a countable collection of points in Q, say {yi,y2, • ■ •}• We 
can also write i}' n A = ^'^^liuk}, which is a disjoint union. We thus have (for each 
i G {1, 2}): /Ugj,($7' n ^) = J2T=i AciVk) = 0, because each term of the sum is zero by 
definition of singular continuous measures. As a result 

A(o') = c^A) = fiHn' n A) - nA) = o, 

for all W G B{Q,), by which we have uniqueness of the decomposition. 

□ 
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Appendix C 



Proof of Lemma 2.3.1 



C . 1 Proof of Part [T] of Lemma [2X1 

This proof is based on [26j, p. 133. 

Proof. By assumption u is a positive measure. Note that the Radon-Nikodym Theorem pro- 
vides the existence of di'/dfi. For simpHcity of notation, let us write: du/dfj, = h. We now 
first prove that the positivity of i' imphes that h ^ almost everywhere w.r.t. 

The proof goes by contradiction. Assume there exists an O G B{Q) satisfying /u(r2) > 
and h < /i-almost everywhere on ^l. We introduce the following sets: 

Fo := {x £n\h{x) <0}, 

Fn := {xenl h(x) ^ --}, for all n G N+. 

' n 

Consider the disjoint union Q = (Q\Fq)[J {^ICi Fq) . From the assumption that h < /i-almost 
everywhere in 0, it follows that /i(il\i<o) = 0. Note that since has strictly positive measure 
fj,{0,), we get 

< = \ Fo) + n Fq) = n Fq). 

Furthermore, we have for all n G 

o^u(nnFn)= I hdfx^ I --di^ = --^i(nnFn) i^o, 

J J n n 

nr\Fn nr]Fn 

which implies that ^{Vt n Fn) = for any n G N^. 

Note that Fo = Ur=i holdsQ Since is a positive measure, this leads to the follow- 
ing inequality: 



^i{^r\Fo) = ^l\^r\i^]Fn) ^ J]] ^(fi n f„) = o. 



n=l / n=l 



^If X £ UJ^i then there exists a specific h G such that x G Fn. The following holds: h[x) ^ 
— 1/n < 0, so i G -Fo- If i G Fo then h{x) is strictly negative, so h{x) ^ — 1/n < for h :— \l/h{x)\. Thus: 
£ G Fft C Ur=i-F'"- 
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where the last step follows from the fact that n Fn) = for any n G N"*". We now have 
a contradiction with the statement H Fq) > 0. Thus, there is no such set Q of positive 
measure, on which h < almost everywhere w.r.t. fi. This implies that h ^ /^-almost 
everywhere in 0. 

Now let h satisfy h{x) ^ for all x G 57 (note that it is possible to choose such a repre- 
sentant). Every nonnegative measurable function is the (pointwise) limit of an increasing 
sequence of nonnegative simple functions; assume that {/ifcj^i is such sequence. Simple 
functions are always measurable, and therefore the following is true: 

lim / hj^dfi = / hdfi, 
fc^oo J J 

n n 

see [26], p. 112. 

This limit is also valid if the domain of integration is any arbitrary Q' G B{^}), as we will show 
now. Let 1^2' denote the characteristic function of $7'. The sequence {h^lQ'} is increasing 
and consists of nonnegative, measurable functions (products of two measurable functions, are 
again measurable). Furthermore In'/i is (by definition of {/ifc}) the pointwise limit of the 
sequence {/ifclc}- Consequently (cf. [26], p. 112): limk^oo J^hklfi'dfi = J^hl^'dfi, which 
can also be written as 



lim / /ifcd/i = / hdfj., for all Q' G B{Q). 

k^oo J J 



For convenience, write dfi/dX = p (the existence of which is guaranteed by the Radon- 
Nikodym Theorem). Following the same arguments as for h, we can conclude that it is 
possible to take a representative of p, such that p{x) ^ for all x gQ. The Radon-Nikodym 
Theorem furthermore ensures the measurability of p. The sequence {h^p}'^^ consists thus 
of measurable functions, as each element is the product of two measurable functions. Both p 
and hk (for all k G N"*") are nonnegative, so the same holds for the elements of the sequence. 
The sequence is increasing, because {/ife} is increasing. As a result, we can conclude that 
linik^oo J^hkpdX = j^hpdX (cf. [26], p. 112). Multiplication with 1^/ and some further 
reasoning yields, like before, that 

lim / hkpdX = / hpdX, for ah Q' G B{Q). 
fc^oo J J 
n' w 

The functions are simple functions, so for each k G there exists a finite collection 
{A^^^}^^^ of mutually disjoint measurable subsets of Jl, where N^^'' G N"*", such that 

Ar(fe) 

hk{x) = ^ al''^l^(k){x), for ah x G 



i=l 



for some collection of coefficients {a-^^}^'^' C For any measurable set A, the following 
identity holds: 

J iAdfi = fi{n' nA) = J pdx = JiAPdX, ioi aim' £ B{n). 



n' fi'nA 
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Therefore 



y hkpdX, for all G B{n) and for all A; G N+. 



J^, hdfi and J^^, hpdX are thus limit values of the same sequence, which means they must be 
equal: ^{Q') = f^, hdfi = J^, hpdX. From this, it follows that the integrand hp is in fact the 
Radon- Nikodym derivative dv/dX, which finishes the proof. □ 



C.2 Proof of Part U of Lemma 2.3.1 



Proof. Define the positive part and the negative part g of the function g by 

5r+(x) := max{c/(x),0}, 
g-[x) := -mm{g{x),0}. 

These two functions are measurable, because g is measurable (see [iBl . p. 15). 

Assume for the moment that neither of g^ and is identically zero. 

If we define for all il' G B{Q) the following: 

7+/- (J]') := / g+/-df,, 



then 7+ and 7" are measures, e.g. due to |l6j (p. 23). Note that by the superscript '+/—', we 
indicate that the statement holds in both the '+' case and the '— ' case. Both g^ and g^ are 
nonnegative, so 7^ and 7" are positive measures. g~^^~{^l) > follows from the assumption 
that g~^ and g~ are not identically zero. Since g G L^(il), it follows that 7^/~(J7) < 00. 
The Radon-Nikodym Theorem can now be applied to the measures 7"^ and 7", and pro- 
vides the existence of d'y^/~ /dfi. The uniqueness of Radon-Nikodym derivatives implies that 
(i7+/-/(i/i = £/+/-. 

Now apply Part [l] of Lemma 2.3.1 by setting i> = 7^ and u = subsequently. Note that 
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7"^/ ^ /X is guaranteed by defining 7+/ as an integral with respect to /i (cf. [26], p. 104). 
The result is 

- ^^'-<"') - / - / - / 

for all rj' € where Part [T] has been used in the third equality. 



Note that 

h' n' 

is automatically satisfied if g^^^ = 0. We can thus proceed without distinguishing between 

g+/- = and 5+/" ^ 0. 

Finally, we conclude that for all 0,' G B{n) 



gdfi = J {g^ -g )dn 

g^dfx - g'dfi 



g~^-j^dX— [ g -j^dX 



dX / " 



by which we have the desired statement. □ 

C.3 Proof of Part [3] of Lemma [2XT 



Proof. Set V = X and apply Part [T| which states 

dX dX dfj, 
dX dfj. dX 

We now claim that dX/dX = 1, which can easily be shown. Obviously A <C A, thus the 
Radon-Nikodym Theorem yields the following: 



X(n') = [ '^dX, for ah n' G B(n). 
J dX 



However, also A(O') = In'dX = J^, IdX holds. The uniqueness of Radon-Nikodym deriva- 
tives implies that dX/dX = 1. We thus have 

^ dX dfi 

dfj, dX ' 
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and hence 

dX \djiJ 



□ 



C.4 Proof of Part H of Lemma [2:3J] 

Proof. Let Vt' G B{Vti x $^2) be arbitrary. Define the fohowing set: 

Oil := {xi € r^i I {xi,X2) G ri' for some X2 € 02}- 

Also define 

{^2 G 1 {xi,x2) G r^'}, 

for any choice of xi E Oi. The set fi^.^ is called the Vti-section of Vt' determined by xi (cf. 
[26] . p. 141). Note that is a subset of O2. In fact, xi can be regarded as a parameter 
here. It can be proved that every section of a measurable set is a measurable set (see |26j . 
p. 141). 

We first prove that vi ® V2 ^J-i ^2- 

Assume that $7' G B{^ll x VI2) is such that (/ii ® fi2)i^') = 0. By definition (see [26], p. 144) 
of the product measure ni ^2 we have 



(^1 fX2){^') = J ^l2{^'xM^'^i^^) = 0- 

As a consequence (cf. [26], p. 147) 

fj,2i^'xi) = 0, for fii-a.e. xi G Sl'^. 

Define 

■.= {x,eQ[\fi2{^',J^0}. 

Let xi be such that xi ^ JIJ^^j then of course /i2(f7^.^) = 0. From the hypothesis that 1^2 ^ 1^2, 
it follows that ^2(0^^) = 0. If we define 

nl^:={x,€n\\u2{n',^)^o}, 

then we have thus proved that Jl^^ ^/l2' 

Since 112(0,'^-^) = for ^ui-a.e. xi G 17']^, we clearly have that fJ-ii^^^) ~ ^- hypothesis of 
the lemma ui is absolutely continuous w.r.t. /ii, and thus fii{Q^^) = implies = 0. 

We now use ilj^^ ^^^2' conclude 

by which we find that = 0. This statement can also be written as 

^2(f^Li) = 0> fo'^ i^i-a.e. xi G fi'^. 
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It follows trivially that 

We have thus shown that z^i (SD 1^2 ^ /^i (E> /i2- 

We can now apply the Radon-Nikodym Theorem, which provides the existence of the unique 
Radon-Nikodym derivative 

d{vi (g) V2) 

Let Vl' G B{Qi X ^12) be arbitrary, and let 1^ denote the characteristic function of the 
set A, where yl is a set in B{Q-i)^ B{Q-2) or B{Vii x 0.2). By hypothesis of the lemma the 
Radon-Nikodym derivatives dvi/d/ii and du2/dii2 exist. We can now perform the following 
calculations: 

{vi®iy2){^') = i^2{K,)di^i{xi) 



ni 



n' 12' 



O' 12' 



di^2ix2) j dui^xi) 

'^dlJ,2{x2) I ^dfii{xi) 
dfi2 I dill 

dvi du2 
dfj.1 d^2 



dfl2{x2)dfJ,l{xi) 



dui dv2 / \\ 7 , \ T / \ 

f^rf^^^'' "^^"1 ld^2(2;2)d^i(xi) 

( ^^"'"'^'1*'^^'* """^^i '^^^^ ) '^^^^ ^ ^^2){X1,X2), 

where the last step follows from [26], p. 147. 

Let (xi,X2) be an element of VLi x O2. By definition of il^- , we have that (X2) = 1 
if and only if (xi,X2) G ^' ■ 

If (xi,X2) G ri', then also 1q^i^{xi) = 1, by definition of $7'^. We thus have 

In'i(Si) = 1- 

On the other hand, if (xi,X2) ^ f^' then (52) = 0, and thus 

^1 

In'i(ii) = 0. 
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We conclude that for all {xi,X2) G x the following identity holds: 

This means that we can write 

(fl (Si iy2){^') = j {'^^'^^'^Q.'{xi,X2)\d{lll®ll2){xi,X2) 

Uniqueness of the Radon-Nikodym derivative now makes sure that 

d{y\ (g) 1/2) _ dv\ dv2 
d{iJ,i (g) /X2) cZ/Lti diJ.2 ' 

which finishes the proof. □ 
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Derivation of the entropy density 
for an ideal gas 



The following derivation is due to [50j. More details can be found in |12) (pp. 16-19 and 
96-102). 

We consider a gas that is compressible, non-viscous and thermal (i.e. the temperature in- 
fluences its behaviour). The set of independent variables is Q := {p,T,j}, where p is the 
density, T the absolute temperature and j = VT, the temperature gradient. All other de- 
pendent variables are functions of Q. We assume the following conservation laws: 



1. Balance of mass: 



|^ + v.(H = o. (D.l) 



2. Balance of energy: 

Dt 

3. Entropy inequality (Clausius-Duhem): 



De 

p— = -pV -v + V ■q + pr. (D.2) 



.^-V.(i)-f,0. ,D.3) 

In the balance of energy, e is the internal energy density. Furthermore, q denotes the heat flux 
and r is the heat supply density. The entropy inequality is the local version of the Clausius- 
Duhem Inequality. It is precisely the inequality given by Postulate |3. 5.1 for one component. 



Moreover = q/T has been taken, as was also suggested by (1.6.9) in [Hj (the results on 

p. 29 have to be reduced to one component to see this). 

Note that we require T > for the entropy inequality to make sense. 



Since (D.l) can also be written as 



dp „ „ 

= -:f+pV-v + Vp-v 
ot 

= (D.4) 
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it fohows that V • v 
Note that 



1 Dp 



De fde ^ 

^ (dp ^ , . 

= p-+pVe.v + e[-+V-{pv) 



=0 

de dp „ , . „ , , 
p-g^ + e-^ + Ve-{pv)+eV -{pv) 



d_ 
dt 



{p£) + V -{pve), 



and a similar identity holds if we take r] instead of e. 

Let F denote the Helmholtz free energy, defined by F := e — rjT. Taking the derivative 
D/Dt in the definition of F and multiplicating afterwards by p, leads to 



De DF DT Dn 



-pV ■ V + 'V ■ q + pr. 



(D.5) 



The latter equality is due to the balance of energy (D.2). Equation (D.5) can also be written 
as 



r^--, DF 



prj— pV ■ V +\/ ■ q + pr. 



(D.6) 



A suitable combination of the entropy inequality (D.3) and (D.6), V • f = and of the 

I I I I p Dt 

fact that r > 0, yields 



1 Dp 



DF DT p Dp 
- P'^Wt ^ -pDi 



or 



pT— - V -q- pr 
Dri 

pT—l - TV 
^ Dt \T 

t(p^-v., 

V Dt \T 



pr 

q_ 

' T 
pr 



VT — pr 

q 



T 



■ VT 



>0 



T 



■VT, 



DF DT pDp 1 



Since F is a depends only on Q = {p, T,j}, we have 



(D.7) 
(D.8) 



DF_dFDp dF DT OF Dj 
15t ~ ~d^lDi^ dflyt ^ IfjWt' 
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by which (D.8) transforms into 
p dF\Dp 



p[r] + 



dF\ DT 
dfJlH 



dFDj 1 



1 



The coefficients of Dp/Dt, DT/Dt and Dj/Dt, and -q ■ VT are independent of Dp/Dt, 

DT/Dt and Dj/Dt. Furthermore, Dp/Dt, DT/Dt and Dj/Dt can be chosen arbitrarily 
(pointwise), and thus the following constitutive equations must hold 



p = p 



,dF 
dp ' 



V 



and since T > also 



dF 
q-VT^O. 



and 



dF 



0, 



(D.9) 



(D.IO) 



Remark D.l. It follows eventually from ( |D.10 ) that q = X{p,T, j)VT, where A ^ 0. If 
linear thermal conduction is assumed (that is, q is linear w.r.t. j), then A = X{p,T) and thus 
Fourier's Law is obtained. However, this result is not so important in the sequel. 



If we want a more explicit form of the constitutive equations in (D.9), we are required to 
give also an explicit choice of F. It is clear from the third constitutive equation in (D.9) that 

F = Fip,T). 

Definition D.2 (Ideal gas). A gas is called ideal gas, if 



F{p,T) := c{T - To) -cT log 



+ RT log 



Po 



where Tq and po are an arbitrary reference temperature and density, c is the specific heat, and 
R is the universal gas constant. 

In fact, c = cy is the specific heat at constant volume, which can only be a function of T 
(cf. [52], p. 114). We have made the assumption that cy is constant. In nature, this assump- 
tion is only valid over wide temperature ranges for monoatomic gases; see [52], pp. 114-115. 
For the sake of clarity, we will not go into details for non-constant specific heat cy = cy (T) . 



For an ideal gas, it follows from (D.9) that 
p = RpT, and 



r/ = c log 



i^log 



po 



In p = RpT we recognize the ideal gas law. Furthermore, note that e can be determined 
explicitly using e = F + rjT. The internal energy is a function of the temperature only: 
e = e{T) = c{T-To). 
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Appendix E 



Modifications in the proofs of 



Theorem 4.5.1 and Corollary 4.5.4 



We describe the details that need to be adapted in the proofs of Theorem 4.5. 1| (global 
existence of the time-discrete solution) and Corollary 4.5.4 (conservation of mass) to make 
them compatible with the new assumptions made in Section |4.5.3[ 



E.l The proof of Theorem 4.5.1 



In Part [2] of the proof of Theorem 4.5.1, the identity {Xn) ^i^) = is used. This identity 
holds if the motion mapping is invertible. Using the pre-image, we only have that 

We modify the proof that the absolutely continuous part fi^c n+i finite accordingly, by 
inserting an inequality: 



00. 



Consider the arguments to prove that ^J-^n+i is finite in Partpjof the proof of Theorem 4.5.1 
There, we used that {xi}i^j C implies {Xn{^i)}ieJ C 0, due to the fact that the motion 
mapping is a homeomorphism mapping from 0, to Q. Note that {xi}i£j C supp /i" holds. In 
Section 4.5.3 we require that {Xn|supp^'» maps from supp^" to O. Thus {Xn{xi)}i€j C is 
still true. 



Regarding Part [4] of the proof, note that 

no longer holds if pre-images are used. Instead we have that 

n\Br.c{x?iy'{^xU^\B, 

We need to prove that if 
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holds, then 



Since 



we have that 



It follows that 



thus 



Now use that 



and 



n\{xZ)-'{xZ{Bn))cn\B^. 
= fi^un \ Bn) ^ (n \ (x^)-^ ix^iBn))), 
^^Z,n{^\ix'^^r'{x'^iBn)))=o. 
ix?^r'm c n, 



to obtain 



/^°c,n(f^\(X°)-HXn(^n))) ^ l^l,n{{x'i)-\^) \ {x'nY^X'^XBn))) 

= f^sc,n{iX?^r'{^\X?iiB^))) 

= fi2,^^Jn\xUBn)). 



Because 



we also have that 



and thus we are done. 



f^?c,n{^\iX?i)-Hx?^{Bn)))=0. 
/^"c,n+lff^\xj^(i?n)) =0, 



We also need to show that X'^{Bn) = implies that A'^(x°(-B„)) = 0. This follows directly 
from the newly made assumption (cf. Section 4.5.3) that there is a constant Cn > such that 



Cn\'^{Xn{^')) ^ {^'), for ah n' e B{n). 

When proving finiteness of fJ-^cn+i again can only use that {Xn)~^i^) ^ instead of 
{xZr'^in) = n. It follows that ' 



oo. 
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E.2 The proof of Corollary |4.5.4 



We redo the proof of Corollary |4.5.4| completely, because it is slightly more involved. 

Proof. For each n € {0, 1, ... , Nt — 1} and a £ {1,2, ... consider the measure where 
CO G {ac, d, sc}. By definition of n+i (^^^ constructive proof of Theorem 
[2l[3]and[4l): 



4.5.1 



Parts 



.a 



For each n we thus have 



11°' 



Note that supp/z^n ^ supp /i" C (x°)~^(f^)- This implies that 



'u),n+l 



= <„((x°)~'(f^)nsupp/x5?)+<„((x^)-i(17)\supp/i^) 

= /^2,n(supp/i^) +/X°„(0) 

= /^2,n (SUPP IJ-n n SUpp /^2,n) + ^2,n (supp Mn \ SUpp //2,n) 



= /^2,n(supp/i2,n) 

By an inductive argument: /i2n(f^) = /"2o(^) ^^^^ each n. 

Since /i" = /^ac.n + A^d n + /^sc,n ^r all n G N, the above implies trivially that 



for all n G AA. 



□ 
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Paper 'Modeling micro-macro 
pedestrian counterflow in 
heterogeneous domains' 

The following pages contain the full content of [231 : our paper 'Modeling micro-macro pedes- 
trian counterflow in heterogeneous domains', published in Nonlinear Phenomena in Complex 
Systems. 
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Modeling micro-macro pedestrian counterflow in heterogeneous domains 



Joep Evers* 
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PO Box 513, 5600 MB Eindhoven, The Netherlands 

Adrian Miintean"'' 
Centre for Analysis, Scientific Computing and Applications ( CAS A ) 
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Institute for Complex Molecular Systems (ICMS) 
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We present a micro-macro strategy able to describe the dynamics of crowds in heterogeneous 
spatial domains. Herein we focus on the example of pedestrian counterflow. The main working 
tools include the use of mass and porosity measures together with their transport as well as suitable 
application of a version of Radon-Nikodym Theorem formulated for finite measures. Finally, we 
illustrate numerically our microscopic model and emphasize the effects produced by an implicitly 
defined social velocity. 

PACS numbers: 89.75.Fb; 02.30.Cj; 02.60.Cb; 47.10.ab; 45.50.Jf; 47.56.+r 
Keywords: Crowd dynamics; mass measures; porosity measure; social networks 



I. INTRODUCTION 

One of the most annoying examples of collective 
behavior[l] is pedestrian jams - people get clogged 
up together and cannot reach within the desired time 
the target destination. Such jams are the immediate 
consequence of the simple exclusion process [2, 3], 
which basically says that two individuals cannot occupy 
the same position a; € C R*^ at the same time 
t E S :=]0, T[, where T g]0,cx)[ is the final moment at 
which we are still observing our social network. 

Observational data (cf. e.g. [4]) clearly indicates 
that such jams typically take place in certain neigh- 
borhoods of bottlenecks [5] (narrow corridors, exits, 
corners, inner obstacles/pillars, ...). The effect of 
heterogeneities [6] on the overall dynamics of the crowd 
is what motivates our work. 

In this paper we start ofT with the assumption that inside 
a given room (e.g. a shopping mall), which we denote by 
f2, there are a priori known zones with restricted access 
for pedestrians (e.g. closed rooms, prohibited access 
areas, inner concrete structures) [7] . whose union we call 
fig. Let us also assume that the remaining region, say 
flp, which is defined by Op := 1] — 0^, is connected. 
Consequently, f2p is accessible to pedestrians. The exits 
of fi - target that each pedestrian wants to reach - are 
assumed to belong to the boundarj^ of ftp. The way we 
imagine the heterogeneity of O is sketched in Figure 1. 



* j.h.m.evers@student.tue.nl; Corresponding author, 
t a.muntean@tue.nl 




FIG. 1. Schematic representation of the heterogeneous 
medium fl. The little black discs represent the pedestrians, 
while the dark gray zones are the parts where the pedestri- 
ans cannot penetrate (i.e. subsets of Qs). The pedestrians 
are considered here to be the microscopic entities, while the 
light grayish shadow indicates a macroscopic crowd; see Sec- 
tion II A for the precise distinction between micro and macro 
made in terms of supports of micro and macro measures. 



In this framework, we choose for the following working 
plan: Firstly, we extend the multiscale approach devel- 
oped by Piccoli et al. [8] (see also the context described 
in [9] and [10]) to the case of coimterflow[ll] of pedestri- 
ans; then we allow the pedestrian dynamics to take place 
in the heterogeneous domain f2, and finally, we include 
an implicit velocity law for the pedestrians motion. The 
main reason why we choose the counterflow scenario [also 



Modelling Crowd Dynamics 115 



2 



called bidirectional flow [12]] out of the many other well- 
studied crowd dynamics scenarios is at least threefold: 

(i) Pedestrians counterflow is often encountered in the 
everyday life: at pedestrian traffic lights, or just 
observe next week-end, when you go shopping, the 
dynamics of people coming against your walking di- 
rection [especially if you are positioned inside nar- 
row corridors]. 



(ii) The walkers trying to move faster by avoiding local 
interactions with the oncoming pedestrians facili- 
tate the occurrence of a well-known self-organized 
macroscopic pattern - lane formation; see, for in- 
stance [13]. 



(iii) We expect the solution to microscopic models posed 
in narrow corridors to be computationally cheap. 
Consequently, extensive sensitivity analyses can be 
performed and the corresponding simulation results 
can be in principle tested against existing experi- 
mental observations [4, 14]. 

The presence of heterogeneities is quite natural. Pedes- 
trians typically follow existing streets, walking paths, 
they trust building maps, etc. They take into account 
the local environment of the place where they are 
located. If the number of pedestrians is relatively high 
compared to the available walking space, then the 
crowd-structure interaction becomes of vital importance; 
see e.g. [15] for preliminary results in this direction. 

As long term plan, we wish to understand what 
are the microscopic mechanisms behind the formation 
of lanes in heterogeneous environments. In other words, 
we aim at identifying links between social force-type 
microscopic models (see [16, 17], e.g.) and macroscopic 
models for lanes (see [13, 18], e.g.) in the presence of 
heterogeneities. Here we follow a measTU~e-theoretical 
approach to describe the dynamics of crowds[19]. Our 
working strategy is very much inspired by the works by 
M. Bohm [20] and Piccoli et al. [8]. 

The paper is organized as follows: In Section II we 

introduce basic modeling concepts defining the mass 
and porosity measures needed here, as well as a coupled 
system of transport equations for measures. In Section 

III we present our concept of social velocity. Section 

IV contains the main result of our paper - the weak 
formulation of a micro-macro system for pedestrians 
moving in heterogeneous domains. We close the paper 
with a numerical illustration of our microscopic model 
(Section V) exhibiting efiiects induced by an implicitly 
defined velocity. 



II. MODELING WITH MASS MEASURES. THE 
POROSITY MEASURE 

For basic concepts of measure theory and their inter- 
play with modeling in materials and life science, we re- 
fer the reader, for instance, to [21] and respectively to 
[9, 20, 22]. 

A. Mass measure 

Let fi C M*^ be a domain (read: object, body) 
with mass. Since we have in mind physically relevant 
situations only, we consider d € {1,2,3}. However, most 
of the considerations reported here do not depend on the 
choice of the space dimension d. Let Umi^') be defined 
as the mass in fi' C fl. Note that whenever we write 
f2' C fi, we actually mean that ii' is such that W (E B{i}), 
where B{Cl) the cr-algelira of the Borel subsets of fi. As a 
rule, we assume to be defined on all elements of B{Q,). 

In Sections II Al and II A 2, we consider two spe- 
cific interpretations of this mass measure that we need 
to describe the behavior of pedestrians at two separated 
spatial scales. 

1. Microscopic mass measure 

Suppose that Q contains a collection of A'" point masses 
(each of them of mass scaled to 1), and denote their po- 
sitions by {pk}k=i ^5 N £ N. We want jim to be 
a counting measure (see Section 1.2.4 in [23], e.g.) with 
respect to these point masses, i.e. for all i]' G B{il,): 

^lrn{n') = #{Pke^-^'}. (i) 

This can be achieved by representing /im as the sum of 
Dirac measures, with their singularities located at the pk, 
k e {1,2,. . . ,N}, namely: 

N 

Mm = X! ^P" ■ (2) 

k=l 

We refer to the measure fim defined by (2) as microscopic 
mass measure. 

2. Macroscopic mass measure 

Let us now consider another example of mass measure 
lim- To do this, we assume that the following postulate 
applies to fim- 

Postulate II. 1 (Assumptions on /im) (i) ^ 0. 
(ii) jim is (J-additive. 
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(Hi) Hm ^'^> where X'^ is the Lebesgue-measure inM''. 

By Postulate II. 1 (i) and (ii), we have that /i,,; is a posi- 
tive measure on 12, whereas (iii) implies that there is no 
mass present in a set that has no volume (w.r.t. A"^). 
A mass measure satisfying Postulate II. 1 is in this con- 
text referred to as a macroscopic mass measure. Radon- 
Nikodym Tlieorem[24] (see [21] for more details on this 
subject) guarantees the existence of a real, non-negative 
density p G L^a(Cl) such that: 

Atm(^^') = / pix)dX'^ix) for all Cl' £ S(0). (3) 

Jn' 

Similarly, we introduce time-dependent mass measures 
Ht, where the time slice t £ S enters as a parameter. 



B. Porosity measure 

Let f2 C M'^ be a heterogeneous domain composed of 
two distinct regions: free space for pedestrian motion 
and a matrix (obstacles) such that f2 = Sl^ U f2p (disjoint 
union), where Ctg is the matrix (solid part) of fl and flp 
is the free space (pores). This notation is very much 
inspired by the modeling of transport and chemical 
reactions in porous media; see [25], e.g. 

Let /Ltp(f2') be the volume of pores in f2' C fl. 

Postulate II. 2 (Assumptions on fip) (i) jip ^ 0. 

(ii) lip is a-additive. 

(iii) < A''. 

By Postulate II. 2 (i) and (ii), we have that /ip is a measure 
on f2. We refer to ^p as a porosity measure (cf. [20]). The 
absolute continuity statement in (iii) formulates mathe- 
matically that there cannot be a non-zero volume of pores 
included in a set that has zero volume (w.r.t. A'*). As- 
sume that £1 is such that A''(f2) < oo. Then the Radon- 
Nikodym Theorem ensures the existence of a function 
(p G L\{0) such that: 

IXp{VL')= I (t)d\'^ for all fi' G B(fi). (4) 
Jo,' 

Note that jip{Q') measures the volume of a subset of fl' 
(namely of f2' fl fip). So, we get that 

= x''{n' n Op) < x^{n') for aii n' e B{n). (5) 

We thus have J^, cjidX"^ < ^, dA'', or ^^,{\ - (P)dX'' > 0. 
Since the latter inequality holds for any choice of 0,' , it 
follows that < 1 almost everywhere in f2. 



C. Transport of a meeisure 

For the sequel, we wish to restrict the presentation to 
the case d = 2. For our time interval S and for each i G 
{1,2}, we denote the velocity field of the corresponding 
measure by v^{t,x) with {t,x) G 5 x il. Let also /ij, and 
/Lt| be two time-dependent mass measures. Note that for 
each choice of i, the dependence on t of is comprised 
in the functional dependence of on both measures 
and This is clearly indicated in (9). 
The fact that here we deal with two mass measures ii\ 
and transported with corresponding velocities and 
v^, translates into: 

'^+v-(My) = o, 

< for ah (t,a;) G 5 X fi. 

(6) 

These equations are accompanied by the following set of 
initial conditions: 

= Ato as t = for i G {1, 2}. (7) 

It is worth noting that (6) is the measure-theoretical 
counterpart of the Reynolds Theorem in continuum me- 
chanics. To be able to interpret what a partial differ- 
ential equation in terms of measures means, we give a 
weak formulation of (6). Essentially, for all test func- 
tions tp^jip"^ G 6*0(12) and for almost every t £ S, the 
following identity holds: 

^ / ^\x)dni{x)= I v\t,x)-Vr{x)diilix) (8) 
Jn Jn 

for all i G {1,2}. 

Definition II. 1 (Weak solution of (6)) The pair 
({/ij }t^Oi (m? }t^o) is called a weak solution of (6), if for 
oZH G {1, 2} the following properties hold: 

1. the mappings t H- ^^^^{x)dii\{x) are absolutely 
continuous for all tp'^ G CQ(f2); 

2. v'eL^[S;Ll,{Q)); 

3. Equation (8) is fulfilled. 

We refer the reader to [26] for an example where the exis- 
tence of weak solutions to a similar (but easier) transport 
equation for measures has been rigorously shown. 

III. SOCIAL VELOCITIES 

We follow very much the philosophy developed by Hel- 
bing, Vicsek and coauthors (see, e.g. [13] and references 
cited therein) which defends the idea that the pedes- 
trian's motion is driven by a social force. Is worth noting 
that similar thoughts were given in this direction (motion 
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of social masses/networks) much earlier, for instance, by 
Spiru Haret [27] and Antonio Portuondo y Barcelo [28]. 
Moreover, other authors (for instance, Hoogendoorn and 
Bovy [29] ) prefer to account also for the Zipfian principle 
of least effort for the human behavior. We do not attempt 
to capture the least effort principle in this study. 



A. Specification of the velocity fields 

Until now, we have not explicitly defined the velocity 
fields {i G {1,2}). Very much inspired by the social 
force model by Dirk Helbing et al. [16], the velocity of 
a pedestrian is modeled as a desired velocity Wj^^ per- 
turbed by a component i 21. The latter component 

iPt 'i^t J 

is due to the presence of other individuals, both from 

the pedestrian's own subpopulation and from the other 
subpopulation. The desired velocity is independent of 
the measures pij and and represents the velocity that 
an agent would have had in absence of other pedestrians. 

For each i (E {1,2}, the velocity is defined by 
superposing the two velocities ii^^g and as 
follows: 



for all t G (0, T) and x G fi. 



(9) 



For a counterflow scenario, the desired velocities of 
the two subpopulations follow opposite directions. We 
thus take 

fixed (for i G {1, 2}) and 



The component Wj^i ^2] models the effect of interactions 

with other pedestrians on the current velocity [30]. Since 
the interactions between members of the same subpopu- 
lation differ (in general) from the interactions between 
members of opposite subpopulations, we assume that 
V? 1 21 consists of two parts: 



+ / r^i\y-x\)g{al^)l^M{y), 



y-x\ 



(10) 



for i G {1,2}, where j = 1 if i = 2 and vice versa. In (10) 
we have used the following: 

• f°^^ and /°PP are continuous functions from M_|_ 
to M, describing the effect of the mutual distance 
between individuals on their interaction. Compare 



the concept of distance interactions defined in [22] . 
yown incorporates the influence by members of the 
same subpopulation, whereas /°pp accounts for the 
interaction between members of opposite subpopu- 
lations. is a composition of two effects: on the 
one hand individuals are repelled, since they want 
to avoid collisions and congestion, on the other 
hand they are attracted to other group mates, in 
order not to get separated from the group. /°pp 
only contains a repulsive part, since we assume that 
pedestrians do not want to stick to the other sub- 
population. 

• al.y denotes the angle between y — x and f^es(^)- 
the angle under which x sees y if it were moving in 
the direction of v'^^^{x). 

• 5 is a function from [— tt, tt] to [0, 1] that encodes 
the fact that an individual's vision is not equal in 
all directions. 

Regarding the specific choice of f°vp g^^d g we 

are very much inspired by [16] and [8], e.g. However we 
do not use exactly their way of modeling pedestrians' in- 
teraction forces. We list here the following forms for the 
functions /°™", /°pp and g that match the given charac- 
terization: 



(i?°PP)2)' ^ ^"""(U) 

if s > i?°PP, 



r""(s) 



ffs<i?°™; ^"^^^ 

1 0, iis> 

g{a) :=<T+(l-f7)^-t^^, for a G [-7r.7r].(13) 

Here F°pp and F°'"'^ are fixed, positive constants. The 
constants i?°PP, {radvi of repulsion) and 

(radius of attraction) are fixed and should be chosen 
such that < < and < i?°PP. Furthermore, 
the restriction max{ii°™°, R°pp} < L has to be fulfilled. 
The interaction /°pp is designed such that individuals 
"feel" repulsion (i.e. /°pp < 0) from another pedestrian 
if they are placed within a distance i?°PP from one 
another. The corresponding statement holds for if 
individuals are within the distance Additionally, 
an individual is attracted to a second individual if their 
mutual distance ranges between i?™" and 
The function g ensures that an individual experiences 
the strongest infiuence from someone straight ahead, 
since g{0) = 1 for any a G [0,1]. The constant a 
is a tuning parameter called potential of anisotropy. 
It determines how strongly a pedestrian is focussed 
on what happens in front of him, and how large the 
influence is of people at his sides or behind him. 



In the remainder of this section, we suggest four 



different alternatives for the definition of vf„i ,,2, by 



118 Modelling Crowd Dynamics 



5 



indicating various special choices of distance inter- 
actions and visibility angles (conceptually similar to 
Oi'^y) as they arise in (10). All of them boil down to 
including an implicit dependency of the actual velocity 
= vL= + 'y? 1 21- Note that this efTect increases 
the degree of realism of the model, but on the other 
hand it makes the mathematical justification of the 
corresponding models much harder to get. 



Modification of the angle a\. 



/opp(|y 

1 

Aim. 

and 
1 

respectively. 



L 
L 



by: 



rP{\{y + vi{t,y)i 



{x + v\t,x)T)\)dT, 

(14) 

{x + v'{t,x)T)\)dT, 

(15) 



We defined the angle al.y as the angle between the 

vector y — X and i'des(^)- However this is not a good 
definition if the pedestrian in position x is not moving in 
the direction of v\^^{x) (or, in a broader sense, if the ac- 
tual speed cannot be approximated sufficiently well by 
the desired velocity). Therefore we suggest to define 
a* = Of* (t) as the angle between y — x and v'^{t, x). 



2. Prediction of mutual distance in (near) future 

Up to now the functions Z"""" and /°pp depended on 
the actual distance between x and y at time t. However 
pedestrians are likely to anticipate on the distance they 
expect to have after a certain (small) time-step (say, 
some fixed At G M). In practice, this means that at 
a time t £ S & person will modify his velocity (either 
in direction, or in magnitude, or both) if he foresees a 
collision at time t + At £ S. 



4- Weighted prediction 

Since an individual probably attaches more value to 
his predictions for points in time that are nearer by than 
others, one additional modification comes to our mind. 
Let h : [t,t + Aimax] — ^ [0, 1] be a weight function. Then 
instead of (14) and (15), we propose 

^rr— / r-i\iy+v'it,y)T)-ix+v\t,x)r)\)hiT)dT, 

^i-max Jo 

(16) 

and 

-XT— / rn\iy+v'it,y)T)-ix+v%x)T)\)hiT)dT. 

(17) 

If h is decreasing, then the influence of ti is larger than 
the influence of tg, if < ^2 (which matches om: intu- 
ition). 



To predict the mutual distance between x and y 
at time t + At, the current velocities at x and y are 
used for extrapolation. The predicted distance is: 
\{y + v{y, t)At) — {x + v{x, t)At)\. Consequently, sticking 
to the notation in (10), the interaction potential f°'"'^ and 
/°PP should depend on \{y + v^{t,y)At) - (x + v'{t,x)At)\ 
and on \{y + v^{t,y)At) — {x + v'''{t,x)At)\ respectively 
(where j = 1 if i = 2 and vice versa). 



3. Prediction of mutual distance within a time interval in 
the (near) future 

The disadvantage of using \{y + v{y,t)At) — (x + 
v(x,t)At)\ is that At is fixed. A pedestrian can thus 
only predict the distance at an a priori specified point 
in time in the future. However, people are able to an- 
ticipate also if they expect a collision to occur at a time 
that is not equal to t + At. We assume now that we are 
given a fixed Atmax £ such that an individual can 
predict mutual distances by extrapolation for any time 
T G {t,t + Attnux)- Thus, Atmax imposcs a bound on how 
far can an individiial look ahead into the future. To cap- 
ture this effect, we suggest to replace /"""(jy — x\) and 



B. Two-scale measures 

We now consider the explicit decomposition of the 
measures /ij and /if. Let the pair (^1,^2) be in [0,1]^, 
and consider the following decomposition of fi^ and fM^: 

4 = diml + (1 - ei)Mi, i e {1, 2}. (18) 

Here, is a microscopic measure. We consider 

{p'k{t)}k=i ^ to be the positions at time t of iV' cho- 
sen pedestrians, that are members of subpopulation i. 
We want to be a counting measure with respect to 
these pedestrians, i.e. for all f2' G B{Q): 

mi{n') = #{pUt)en'}, ie{h2}. (19) 

We thus define mj as the sum of Dirac masses (cf. Section 
II A 1), centered at the pi, k = 1,2, . . . , N^: 

N' 

< = ^^lit)' »e{l,2}. (20) 

fc=i 

is the macroscopic part of the measure, which takes 
into account the part of the crowd that is considered 
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continuous. We consequently have <C A^, since a set 
of zero volume cannot contain any mass. Note that we 
are thus in the setting of Section II A 2. Now, Radon- 
Nikodym Theorem guarantees the existence of a real, 
non-negative density p^{t, •) G L]^2{^) such that: 



Jn' 



x)dX'^{x) 



for all O' G B{n) and all ie {1,2}. 



IV. MICRO-MACRO MODELING OF 
PEDESTRIANS' MOTION IN 
HETEROGENEOUS DOMAINS 



(21) 



porosity measure). A basic property of Radon-Nikodym 
derivatives now gives us: 



dMi _ dMl djjLp 
dX^ djjLp d\^ 

for each i G {1, 2} and for almost every i G 5. 



(22) 



We have already defined p^{t,-) := 



and ( 



dUp 



. If 



dA2 ^' dA2" 

we now denote by p^{t,-) the Radon-Nikodym derivative 



dMl 
d^ip 

for ah i G {1,2}. 



, the following relation holds: p^{t,-) = p^{t,-)(j){-) 



We have already made clear that we want to model the 
heterogeneity of the interior of the corridor. In practice 
this means that pedestrians cannot enter all parts of the 
domain. As described in Section II B, we have a measure 
Pp corresponding to the porosity of the domain (which 
is fixed in time). However, we note that the concept of 
porosity (cf. Section II B) is a macroscopic one. For this 
reason only the macroscopic part of the mass measure in 
(18) needs some modification with respect to the poros- 
ity. In this context, one should be aware of the analogy 
with mathematical homogenization. This technique dis- 
tinguishes between microscopic and macroscopic scales, 
where we also see that some (averaged) characteristics 
are only defined on the macroscopic scale. For more de- 
tails, the reader is referred to [25] or [31]. In M.'^, we 
have Pp <C A^. Furthermore 4; Pp for i G {1,2} 
and a.e. t £ S. This is obvious, since no pedestrians 
can be present in a set that has no pore space (i.e. zero 



A. Weak formulation for micro-macro mass 
measures 

We now have the following measure; 

pi = Oiml + (1 - ei)Mi, i G {1, 2}, (23) 

as was given in (18), where now: 

n' 

ml = Y.^Plit)' dMiix) = p\t,x)cj>{x)d\^{x). 

fc=i 

(24) 

This specific form of the measure will now be included 
in the weak formulation (8), with velocity field (9)-(10). 
The real positive numbers 6i {i G {1,2}) are intrinsic 
scaling parameters depending on A/"'. 

The transport equation (8) takes the following form: 



d 
dt 



(e, V V'(pi(t)) + (1 - Oi) [ ^\x)p\t,x)ct>ix)dX''ix)) = 
^ k=i ■'^ ' 

eiY^^'it^Pkit)) ■ v^A'HW) + (1 - / v\t,x) ■ V7p'{x)p\t,x)4>ix)dx^{x), 



(25) 



for alH G {1, 2}. Here we have used the sifting property 
of the Dirac distribution. In the same manner, we specify 
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fc=i ' ' 

/ r%\y-A)9{<y)^;^/{t.y)ct>{y)dX\v) 
Jn \y ~ ^1 

+0j E /°''''(bi(*)-^l)ffK 



+ / r^^^^{\y-x\)g{al^)j^i^{t,y)<p{y)dX\y), (26) 

12/ ~ X\ 



for i e {1,2}, and j as before (j = 1 if i = 2 and vice 
versa). We have omitted the exclusion of {a;} from the 
domain of integration (in the macroscopic part), since 
{a::} is a nuUset and thus neghgible w.r.t. A^. Note that 
the sums may be evaluated in any point a; G f2 (not 
necessarily x = p\.{t) for some i and fe); the integral parts 
may also be evaluated in all x, including x = p\.{t) for 
some i and k. 



V. NUMERICAL ILLUSTRATION 

We wish to illustrate now the microscale description 
of a counterflow scenario (i.e. for 9i = 62 = 1) by 
presenting plots of the configuration of all individuals 
situated in a given corridor at specific moments in time. 

We consider a specific instance in which there are 
in total 40 individuals (20 in each subpopulation) . The 
dimensions of the corridor are c? = 4 and L = 20. The 
velocity is taken as defined in (10)-(13). Furthermore, 
the following model parameters are used: v\^^ = 1.34ei, 
vl^^ = -1.34ei, F°PP = 0.3, F°™ = 0.3, i?°PP = 2, 
^own ^ 2, R°^^ = 0.5, = 0.5, = 0.5, CT = 0.5. 

In Figure 2, we show the configuration in the cor- 
ridor at times t = Q,t = 7.5, and t = 15. The individuals 
of the subpopulation 1 are colored blue, while the 
individuals of the subpopulation 2 are colored red. 
Clearly, self-organization can be observed in the system: 
Pedestrians that desire to move in the same direction 
form lanes (in this case, three of them). This feature is 
observed and described extensively in literature, cf. e.g. 
[13]. 

Another feature, pointed out by Figure 2, is the 
following: Within the three already formed lanes, small 

clusters of people are formed. This flocking is a result of 
the typical choice for in (12). Members of the same 



subpopulation are repelled if their mutual distance is in 
the range (0,i?°™); they are attracted if their mutual 
distance is in the range (J?"""", i?°^"). No interaction 
takes place if individuals are more than a distance R'^^ 
apart. 

The attraction part of the interaction causes indi- 
viduals that are already relatively close to get even 
closer, until they are at a distance i?"""". For distances 
around i?™", there is an interplay between repulsion 
and attraction, eventually leading to some equilibrium 
in the mutual distances between neighboring individuals 
in one cluster. In Figure 2, we observe self-organized 
patterns even within clusters. 

We should mention here that further simulation is 
needed, in particular on multiply-connected domains. 
Note however, that as the model is now, moving individ- 
uals are not a priori hindered to penetrate any obstacle. 
Moreover we wish to point out that the number of 
particles should be sufficiently large, in order to be 
statistically acceptable. 
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